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SUMMARY 

This is the first of two volumes devoted to a computer program 
for predicting the performance of an annular combustor. In this volume 
the correlations and calculation methods used in the program are presented, 
and the assumptions involved in applying them are discussed. The results 
of a number of test cases are also described. 

The analys's falls naturally into three parts: 

1. Performance of the diffuser. 

2. Air-flow, pressure, and temperature distributions in the 
flame tube and annuli, including calculation of the combustor 
total -pressure loss. 

3. Heat transfer and the calculation of the f lame-tube-wal 1 
temperature distribution. 

Each of these parts is discussed in detail in Volume I of this report. 

Volume 11 gives a complete description of the computer program, 
forming a self-contained operating manual that permits the program to be 
operated without reference to Volume I. 
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INTRODUCTIO N 

The cut-and-t ry methods generally used to design combustion 
chambers have been giving way in recent years to a more systematic approach 
based on analysis and correlation of experimental data. This approach 
is potentially capable of substantially reducing development time and 
expense while improving combustor performance. While the analytical 
methods and correlations presently available are crude in many respects, 
they are nonetheless sufficiently complex that computers can effectively 
be used in their solution. 

The aim of the work described in this report has been to develop 
a computer program for the analysis of fluid flow, combustion, and heat 
transfer in annular combustors with diffusers, making use of currently 
available analytical methods and correlations. It is expected that the 
resulting program, used carefully in conjunction with experimental data, 
will prove to be a useful tool for the design of high-performance annular 
combustors . 

The program has been extended to include analysis of rectangular 
test combustors. An approximate analysis of tubular geometries may be ob- 
tained from the program for rectangular combustors. 

Background 

A number of attempts have previously been made to calculate the 
air flow and heat transfer in gas-turbine combustors. Graves (Ref 1) and 
Grobman (Ref 2) analyzed the pressure loss and air-flow distribution in 
tubular combustors with uniformly tapered (or constant cross-secti or ) flame 
tubes. This work included the effects of momentum transfer between the 
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gas streams in annulus and flame tube, annulus wall friction, heat release, 
hole discharge coefficients, and compressibility. Flow in the diffuser 
was not considered, and instantaneous mixing between cold-air jets and 
the gas in the flame tube was assumed. No heat-transfer analysis was 
undertaken. The results of these calculations exhibited many of the effects 
observable in practical combustors. 

Samuel (Ref 3) used the same general approach as Grobman, with 
instantaneous mixing between cold-air jets and the main gas stream. A 
simple diffuser analysis was included, and ; oth tubular and annular 
geometries could be treated. In the case of annular combustors, the 
calculation followed the three parallel streams in the inner and outer 
annuli and through the flame tube, and iterated on the initial mass-flow 
split between these streams until boundary conditions at the end of the 
combustor were satisfied. A similar approach is used in the present work. 

The output of Samuel's program was used as input for a heat-transfer program. 

The most comprehensive heat-transfer analysis of aircraft-type 
combustors that has been published to date is that of Lefebvre and 
Herbert (Ref 4). This involved the solution of a heat-balance equation 
for each element of the wall, taking into account radiation from flame to 
wall and wall to casing, and forced convection on the inside and outside 
of the wall. Net heat transfer in the longitudinal direction was assumed 
to be negligible compared with radial heat transfer. Reasonable agreement 
with limited experiment was obtained. 

Another heat-transfer analysis was carried out by Tipler 
(Ref 5). whose work dealt with an industrial combustor and included the 
effects of radiation transfer in a longitudinal direction. 

Over the last few years, a number of new correlations and 
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techniques have been produced which provide the basis for a more 
sophisticated approach to the analysis of flew and heat transfer in 
combustors. 

Sovran and Kiomp (Ref 6) have provided a way of generalizing straight- 
walled two-dimensional and annular diffuser performance. An almost identical 
approach was used by Reneau, Johnston, and Kline (Ref 7)* who extended the 
work to cover different inlet boundary- layer thicknesses; this aspect is 
particularly important for gas-turbine combustor diffusers. For the first 
time, therefore, there now exists a set of maps which enables the perfor- 
mance of straight-walled diffusers to be predicted with fair certainty for 
any length and area ratio and inlet boundary- layer thickness, at least 
within the range of interest. 

Although many sets of data on the behavior and mixing of jets 
have been accumulated over the last 20 years, it is only recently that 
such data have been brought together, correlated, and presented in a 
generalized fashion. For example, a review of 16 experiments on jets 
penetrating a gas stream at an angle has recently been prepared by Northern 
Research. Spalding (Ref 8) refers to 18 experiments on wall jets in 
developing a new theory for this type of flow. 

In the field of heat transfer, Spalding (Ref 9) has correlated 
the results of 9 experimenters to produce a new correlation for film- 
cooling heat transfer, and Sturgess (Ref 10) has shown that this correla- 
tion works moderately well for the practical film-cooling slots used in 
actual combustors. Schirmer and Q,uigg (Ref 11) have produced the most 
useful work yet published on the effect of pressure on radiation from 
luminous flames. Finally, Gruber (Ref 12) and Sparrow (Ref 13) have shown 
how view-factor data can be used to make a more comprehensive analysis of 
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radiative heat transfer in combustors than has hitherto been possible. 

The works mentioned above, and others, provide useful new tools 
for combustor analysis, and open the way to removing some restrictions 
that have previously limited the utility and validity of such analyses. 

Overall Objectives of Present Program 

A computer program was therefore envisaged to achieve the 
following objectives: 

1. Analyze the air mass-flow and pressure distributions in a 
combustor whose geometry is given, for specified inlet 
cond i t i ons . 

2. Compute the temperature distribution on the flame-tube wall. 

3. Predict whether separation occurs in the diffuser. 

The rest of this report discusses thesci objectives further, and shows how 
they are achieved. 

Arrangement of Report 

The section following this introduction contains a summary of 
the overall approach used in this work, including the major assumptions. 

In the three subsequent sections are described the calculation methods and 
organization of the diffuser, air-flow, and heat-transfer analyses. Program 
limitations are donsidered in the next section. 

The computer program has been applied to a number of sample cases 
designed to test individual subprograms as well as the overall program. 

The results of these cases are discussed in the next section. A final 
section contains conclusions derived from this work. 

Detailed derivations, where not germane to the main discussion, 
are relegated to appendices. 
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OVERALL APPROACH 


I introduction 

This section is a general introduction to the computer program 
and the methods and assumptions upon which it is based. It is meant to 
provide the reader with a framework for the detailed discussions that follow 
in succeeding sections. Topics covered are; 

1. Major assumptions* 

2. Combustor geometry and nomenclature. 

3. Input to the computer program. 

4. Calculation options. 

5. Functions of major subprograms; control, diffuser, air-flow, 
and heat-transfer. 

6. Interaction between major subprograms; major interative 
procedures. 

7. Output of the computer program. 


Major Assumptions 

Assumptions affecting the program as a whole are listed below. 
Assumptions affecting individual parts of the program are introduced in 
the applicable detailed sections. 

1. Conditions are steady; no quantity varies with time. This 
assumption is reasonable but it precludes the treatment of 
transient phenomena such as combustion oscillation. 

2. Combustors are of annular or rectangular cross section. 

This precludes the direct simulation of tubo-annular (cannular) 
and tubular combustors. While flow conditions for these 
types may be approximated by appropriate selection of 
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dimensions, other parts of the simulation (for example, 
radiation heat transfer) will be of reduced accuracy. 

3. Conditions are uniform around the circumference of the 
combustor; i.e. there is no circumferential variation in 
any quantity. This assumption is a major simplification; 
in practice, cyclic variations around the circumference 
occur, corresponding to the fuel-injection points. The 
results of this analysis will be useful, however, in 
indicating output parameters averaged around the circum- 
ference. In the case of combustors of rectangular cross 
section, this assumption implies that end effects 

due to the finite width of the combustor are ignored. This 
does limit the extent to which calculated results can be 
expected to agree with experiment. 

4. The flow can be treated as quasi -one-dimensional. To 

represent the complex flow pattern existing in a real 
combustor is beyond Uie power of present analytical methods. 
Instead, the problem is broken down into two parts: flow 

in the main gas streams (one-dimensional) and flow and 
mixing of jets (predicted from correlated experimental data). 

5* Radial -flow or reverse-flow combustors cannot be 
analyzed. 

6. Flow conditions in the combustor are within the range of 
validity of the correlations used. The correlation are 
based on experimental data presently available. Outside 
their range of established validity, their accuracy is 
likely to be reduced. 
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Combustor Geometry and Nomenclature 

The geometry of an annular combustor and the major terms 
used in describing it are shown in the sketch below: 

Secondary Cooling 



by the snout into three streams, which flow into the snout and the outer 
and inner annuli of the diffuser. After diffusing, the three streams enter 
the flame tube: the snout air through the swirler and the holes in the 
dome, the inner- and outer-annulus air through holes and cooling slots 
in the flame- tube walls. 

Swirler 

A swirler is essentially a bladed passage, similar to a fan 
stator, which imparts swirl to the air flowing through it. Individual 
fuel injectors (of which there are many in an annular combustor) are often 
surrounded by swirlers to promote rapid mixing of air with the fuel spray. 
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The computer program includes a rather crude method for selecting swirler 
dimensions. A method for calculating swirler pressure drop is also 
i ncluded. 

Holes and Slots 

Apertures in the walls of the flame tube and in the dome are 
classified as penetration holes or cooling slots, according to their 
function. Penetration holes direct the entry of air jets into the flame 
tube for mixing and dilution of the hot combustion gases. Cooling slots 
introduce films of cool air just inside the flame-tube walls to act as 
thermal barriers. The cooling slots are assumed to be of the circumferentially 
continuous, dynamic-head type shown in the sketch below: 

Cooling Air __ j 



Hot Gases 


Correlations for the flow and mixing of penetration and wall jets are part 
of the computer program, as are film-cooling correlations. 

Secondary Holes and Primary Zone 
In this report the holes whose upstream edge marks the down- 
stream boundary of the primary zone of the flame tube are defined as 
secondary holes. (Holes farther downstream are referred to as dilution 
holes.) Part of the air from the secondary holes recirculates upstream 
into the primary zone; the remainder passes downstream and is treated 
as ordinary pane t rat ion-jet air. It should be noted that recirculation 
is considered to occur only at the secondary holes, and that all secondary 
holes have their center-lines at one axial location in the combustor. The 
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secondary holes n^ed not to be the first holes in the flame-tube wall. 
Input to the Computer Program 

Input quantities required by the computer program are of three 

kinds : 

1 . Library data. 

2. Combustor geometry. 

3. Case data (inlet flow conditions, calculation options, and 
so forth). 

These quantities will now be described in general terms. A detailed 
description of the input format, units, and so forth, is given in Volume 
I I . 


Library Data 

The library data are tables of quantities available for use by 
the program: 

1. Discharge coefficients and jet angles for 100 hole types. 

2. Performance data for straight-walled two i imens iona 1 diffusers 
and straight-walled annular diffusers. 

3. FI ame-emi ss i vi ty data (none initially supplied). 

The library data may be altered or supplemented by the user; appropriate 
instructions for so doing are given in Volume II. 

Combustor Geometry 

The following information must be supplied by the user to specify 
the combustor geometry. 

1. Selection of annular or rectangular type. 

2. Dimensions of walls of diffuser, snout, dome, flame tube, 
and casing. 
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3. Hole and cooling-slot data: 

a. Hole type (selected from library data). 

b. Number of holes in row. 

c. Axial location; inner or outer wall. 

d. Area (or cooling-slot height). 

4. Identification of secondary-hole row (marks end of primary 
zone) and first hole row or slot in the flame-tube wall, 
as distinct from the dome (marks end of diffuser and start 
of combustor annulus). 

5. Specification of swirler (optional): 

a. No swi rler. 

b. Swirler designed within program. 

c. Swirler dimensions supplied by user. 

Case Data 

The following information must be supplied for each case considered 
by the computer: 

1. Flow conditions at inlet (compressor exit): 

a. Total temperature. 

b. Total pressure (weight-mean average). 

c. Air mass flow rate. 

d. Boundary- layer blockage and shape factors. 

e. Veloci ty prof i le. 

2. Other flow conditions: 

a. Initial estimate of flow split among the snout and the 
inner and outer diffusing passages. (Optional. If not 
specified, program uses internal estimate.) 
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b. Initial estimate of downstream blockage in diffuser. 

c. Fraction of secondary-hole air recirculating upstream 
into primary zone. (Optional. If not specified, program 
calculates flow recirculation.) 

d. Bleed air flowing from combustor annuli (optional): 

i. Location (up to three locations). 

ii. Fraction of inlet-air mass flow bled. 

3. Information on combustor walls: 

a. Emissivity and absorptivity of flame-tube walls. 

b. Thermal conductivity and thickness of flame- tube walls. 

c. Emissivity of casing wails. 

d. Temperature distribution of casing walls (optional; if 
not specified, assumed equal to compressor-discharge 
temperature). 

4. Fuel characteristics and distribution: 

a. Overall fuel-air ratio. 

b. Lower calorific value. 

c. Hydrogen-to-carbon ratio. 

d. Distribution of fuel burning rate, specified as an 
upper limit to the fraction of fuel burned up to each 
hole row. 

5. Choice of optional calculation methods to be used within the 
program, and specification of data peculiar to the correla- 
tions employed in these methods. These options are discussed 
below, in connection with the major subprograms in which they 


appear. 
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6. Miscellaneous data: 

a. I terati on 1 imi ts. 

b. Calculation-step sizes. 

Overall Structure of Computer Program 

The computer program comprises 35 subroutines; these are conveniently 
grouped by function into four major subprograms: 

1. Control subprogram. 

2. Diffuser subprogram. 

3. Air-flow subprogram. 

4. Heat-transfer subprogram. 

The general functions of these subprograms and the flow of information among 
them are summarized in flow-chart form in Figure 1. Their operation is 
discussed in the following four subsections. 

Control Subprogram 

The control subprogram performs the preliminary operations 
necessary to prepare the computer for the main calculations. Its functions 
may be described as follows: 

1. Read in and print out the case data and the combustor 
geometry. 

2. Read in the library data; assemble and print out a short 
list of data (such as hole dimensions) actually needed for 
the cases being considered, 

3. Calculate the combustor reference area (the maximum cross- 
sectional area bounded by the inner and outer casings). 

4. Set up the system of geometrical indexing used in the 
diffuser, air-flow, and heat-transfer subprograms. 


As familiarity with the geometrical indexing systems is basic to under- 
standing the details of the various calculation methods, this subject is 
discussed in detail below. 


Geometrical Indexing System 

Four interrelated geometrical indexing systems are used within 
the program: 

1. Input stations. 

?. Hole rows. 

3. Diffuser stations. 

4. Flame-tube calculation stations. 

These are illustrated by the sample numbering system on the sketch below: 


Geometr i c 
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Input Stations 

The shapes of the walls of the diffuser, snout, dome, flame 
tube, and casing are described by their axial and radial dimensions at 
geometric input stations. Axial dimensions are measured from the diffuser 
inlet (compressor exit); radial dimensions from the engine center-line. 

For rectangular combustors widths are measured from an arbitrary datum 
point. Station numbers begin with Station 1 at the compressor exit. 

Hole Rows 

Holes are indexed in groups referred to as hole rows. A hole 
row comprises holes of identical size and shape, centered at the same axial 
location, and spaced at equal intervals around a single wall of the flame 
tube or dome. Up to six hole rows may be specified at a given axial 
location: hole rows may be interspersed on the same wall or paired on 

opposite walls of the flame tube. Indexing starts with the first hole 
row downstream of the compressor exit. 

The indexing system for slots is similar to that for holes, 
except that slots are circumferentially continuous (as mentioned above) 
and hence only one slot may be placed on a given wall at any axial location. 
Noncont i nuous siots may be approximated by continuous slots of equal area. 
Slots and holes are indexed together. 

Dummy hole rows (holes of zero diameter) may be specified as 
input, for purposes of arranging the location of calculation stations (see 
the subsequent section on calculation stations). 

Diffuser Stations 

The calculations in the diffuser make use of the input-station 
indexing system. In addition, the diffuser is divided into major regions 
by a separate indexing system, labeled as follows in the above sketch: 
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1 Entrance to the diffuser (compressor exit). 

2 Last geometric input station upstream of the start of the 
snout. 

2' First geometric input station on the snout. 

3 Last geometric station where diffusion takes place. 

4 Entrance to the annulus (diffuser exit). This station 
coincides with the upstream edge of the first hole row in 
the flame tube, as distinct from the dome. 

Flame-Tube Calculation Stations 

The primary zone of the flame tube is bounded by the dome and 
swirler, the flame-tube walls, and (at its downstream end) by a plane 
passing through the upstream edges of the secondary holes. This zone is 
treated as a homogeneous stirred reactor; no calculation stations are 
used . 

Calculations in the annulus and in the flame tube downstream 
of the primary zone are related to a system of calculation stations 
located as follows: 

1. At the upstream edge of each hole row that is centered at 
a distinct axial location. Thus, for hole rows paired on 
opposite flame-tube walls only one calculation station is 
ass i gned . 

2. At each dummy hole row. 

3. At the axial location of each cooling <Mot having a distinct 
axial location. 

4. At up to five axial locations spaced at a specified interval 
downstream of each cooling-slot location. 

5- One station at the very end of the combustor. 

Indexing of calculation stations starts with the first hole row in the 
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flame tube (as distinct from the dome) downstream of the compressor exit. 

The first calculation station thus coincides with the end of the diffuser 
annuli (Diffuser Station 4) and the start of the combustor annuli. 

Diffuser Subprogram 

The diffuser subprogram receives from the control subprogram 
inputs of geometry, inlet flow conditions, and certain quantities peculiar 
to the calculation option being used; it also receives a first estimate 
of the mass-flow split among the snout and the two annuli from the control 
subprogram, and subsequent estimates (during iterations on the flow split) 
from the air-flow subprogram. 

Using these inputs and the calculation options chosen by the 
user, the diffuser subprogram performs the following functions: 

1. Determines diffuser performance parameters: 

a. Ideal pressure-recovery coefficient. 

b. Actual pressure-recovery coefficient. 

c. Effectiveness. 

2. Evaluates flow mis-match at the snout. 

3. Determines whether or not flow separation occurs and, if so, 
i ts locat i on. 

4. Calculates flow conditions on the dome and at the ends of 
the d i f fuser annul i : 

a. Static pressure. 

b. Static temperature. 

These conditions are supplied as input to the air-flow 
subprogram. 

For convenience in performing calculations the diffuser is divided 
into three diffusing regions (between Stations 1-2, and 2' -3 in the two 
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annuli, as shown in the sketch on page 14), and two mixing regions (between 
Stations 3 and 4). Within the diffusing t sgions the three diffuser calculation 
options (empi ricai -data method, streamtube method, and mixing-equation 
method) may be used in any of the following combinations: 


Inlet Reqion (1 to 2) 

Annular Reqions (?.' to 3) 

Streamtube 

Streamtube 

Streamtube 

Empirical Data 

Streamtube 

Mixing Equation 

Empirical Data 

Empirical Data 

Empirical Data 

Mixing Equation 


For diffusers without snouts, diffusion in the annular regions 
2 to 3 is treated by assuming that the flow expands isentropical ly from 
the dome stagnation pressure. 

The mixing-equation method is always used in the mixing regions 
(Stations 3 to 4). In addition, area mismatch at the snout (Stations 2 

to 2') is accounted for in all methods by calculating the total -pressure loss 

due to sudden expansion or contraction. 

* 

The three calculation options will now be briefly described. 

Empi ricai -Data Method 

This method is based on the direct evaluation of the effectiveness 
of a particular diffusing passage. There are two ways this can be done, 
at the user's option: 

1. Diffuser effectiveness may be supplied directly as a program 
input. 

2. Diffuser effectiveness may be calculated by the program from 



19 


correlations of experimentally measured effectiveness versus 
diffuser geometry that are contained in the library data. 

The results are corrected for deviation of the inlet 
boundary- layer blockage from that present in the experiments 
on which the correlations are based. 

The remaining diffuser-performance parameters (ideal and actual 
pressure- recovery coefficients) and the outlet conditions are calculated 
directly from the effectiveness, the geometry, and the assumption that the 
total pressure is constant. 


Streamtube Method 

In this method the flow passage is divided into N streamtubes 
within which flow is assumed to be uniform and isentropic. In setting up 
the streamtubes, the distribution of boundary- layer displacement thickne. 3 
supplied as a program input is used to adjust the flow area. The develop- 
ment of the flow (static pressure and velocity) is computed for each 
streamtube at each station along the diffusing passage. Boundary- layer 
calculations are performed to obtain a revised estimate of boundary- layer 
displacement thickness, which is then used in a revised streamtube analysis. 
This process is continued to convergence of the displacement thickness. 

The result of this method is the velocity distribution, static pressure, 
and blockage at the outlet from the diffusing passage. 

Mixing-Equation Method 

This is a rather crude method that can only be used between 
Stations 2 and 4. It takes into account the following effects occur ing 
in these passages: 
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1. Pressure loss due to sudden expansion or contraction at the 
snout (Stations 2 to 2'). 

2. Mixing in the diffusing passage from Station 2' to 4. 

3. Pressure loss due to curvature of the flow passage from 2' 
to 4. 

A flow diagram of the diffuser subprogram is given in Figure 2. 
Air-Flow Subprogram 

The air-flow subprogram receives the following inputs from the 
control subprogram or program input: 

1. Geometry of walls, holes, and swirler (if specified), at 
appropriate indexing stations. 

2. Fuel data. 

3. Jet-mixing model and entrainment constant. 

4. Initial estimate of mass- flow split among snout and annuli. 
From the diffuser subprogram it receives the static pressures and static 
temperatures on the dome and at the ends of the diffuser annuli. 

Using these inputs, the air-flow subprogram performs the following 

functions: 

1. Calculates the flow conditions in the primary zone (static 
pressure and temperature, mass-flow rate) and at calculation 
stations in the flame tube and annuli (static and total 
pressure and temperature, mass-flow rate, and velocity). 

2. Calculates the combustor total -pressure loss. 

3. Directs the iterative process by which the diffuser and air- 
flow subprograms together arrive at the correct flow split 
among the snout and the two annul i . 
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Opti ons 

The following options are available in the air-flow subprogram 
and must be specified by the user: 

1. Jet-mixing models. The air-flow subprogram contains four 
models to represent the mixing of penetration and wall jets 
into the main gas stream. The user must choose a model, 
which is then used for both types of jet. Associated with 
each model is a constant that determines the actual rate of 
jet mixing. An individual value of this constant must be 
specified for each jet type. 

2. Heat transfer from the flame-tube walls to the annulus air- 
may be calculated or ignored at the user's option. If this 
quantity is calculated, an iteration in the overall calculation 
is required, because of the influence of heat input on the 
flow conditions in the annulus. 

Calculation Procedure 

The main steps in the calculation procedure within the air-flow 
subprogram are as follows: 

1. Calculate the primary-zone pressure from the pressure on the 
dome, the snout air- flow rate, and the relations between 
pressure drop and flow rate for the swirler and the dome 
holes. 

2. Calculate the air-flow rate into the primary zone through 
holes up to and including the secondary holes (recirculating 
part) . 

3. Calculate the primary-zone exit static pressure and tempera- 
ture from the total flow rate, inlet conditions, and fuel 
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burned (up to stoichiometric). The primary zone is 
treated as a stirred reactor. 

4. Proceed from calculation point to calculation point down 

the annuli and flame tube, calculating the mass-flow rates, 
velocities, and static and total pressures and temperatures. 
At the end of this calculation procedure the flows at the ends 
of the annuli are inspected. If these flows are not correct (equal to 
zero, or the correct bleed flow, if any), the mass-flow split within the 
diffuser is altered according to the following scheme: 

1. If the flows in the two annuli are o* opposite sign, the 
annulus flows are adjusted, while the dome flow is held 

f i xed . 

2. If the flows in the annuli are of the same sign, the dome 
flow is altered and proportionate changes in annulus flows 
are made. 

The entire sequence of diffuser and air-flow calculations is then repeated. 
This process is continued until the residual air flows at the ends of 
the annuli approach zero within a specified tolerance, or until the 
specified interation limit is reached. 

A flow diagram for the air-flow subprogram is given in Figure 3* 

Heat-Transfer Subprogram 

The heat-transfer subprogram receives as input the geometry of 
<-he fiame tube and casing from the control subprogram and the axial 
distributions of velocity and temperature of the flame-tube gases from the 
air-flow subprogram. Additional data that may be specified as input at 
the option of the user are: 
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1. Temperature distribution of the outer casing. 

2. Choice of calculation options. 

3. Data needed for the calculation option chosen: 

a. A constant characterizing cooling-film effectiveness. 

b. The permeability coefficient for porous walls. 

c. Flame-emissivi ty data (optional). 

d. Thermal conductivity and thickness of flame- tube walls. 

4. Fuel hydrogen/carbon ratio. 

5. Flame-tube and casing absorpt i vi t ies and emissivi ties. 

The primary objective of the heat-transfer subprogram is to 

establish the axial distribution of temperature along the flame-tube walls 
for various program options. Under one program option, a further objective 
is to provide the air-flow subprogram with the axial variation of the heat- 
transfer rate from the flame-tube wall to the annulus air. 

In operation, the subprogram evaluates various heat-flux components 
at a point on the flame-tube wall in terms of the wall temperature; the 
heat-flux components which are considered may include several or all of 
the following, depending upon the program options specified: 

1. Convection from the flame-tube gases. 

2. Convection to the annulus air, 

3. Radiation from the flame. 

4. Radiation to the outer casing. 

5. Radiation interchange between the flame- tube walls. 

6. Longitudinal conduction along the flame-tube wall. 

7< Heat transfer to transpiring air for porous wal ls. 

The heat-balance equation is then solved to determine the unknown tempera- 
ture. These operations are performed at calculation stations along the 
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flame tube; no provision is made, however, for the inclusion of the flame- 
tube dome in the heat-transfer analysis. 

Program Options 

Several of the program options concern the heat-transfer subprogram; 
these options are listed below. 

1. The flame- tube wall may be cooled by: 

a. Convection to the annulus air and radiation to the 
outer casing. 

b. Combined effects of (a) and film cooling. 

c. Combined effects of (a) and transpiration cooling. 

2. The emissivity of the flame may be computed from: 

a. One of two equations for nonl urn i nous flames. 

b. One of three equations for luminous flames. 

c. Emissivity data supplied as special input by the user. 

3. The radiation from the flame to the flame- tube wall may be 
calculated on the basis of one of two assumptions: 

a. Radiation is purely radial (one-dimensional). 

b. Radiation to an element of the wall occurs from all 
(radial and axial) parts of the flame (two-dimensional). 

4. The heat-balance equation may: 

a. Include a term for heat conduction along the flame-tube 
walls. 

b. Include a term for radiation interchange between the two 
flame-tube walls. 

c. Include the terms in both (a) and (b). 

d. Neglect the terms in (a) and (b). 

5. An option is available to exclude the heat-transfer subprogram 
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from the program routing for cases in which heat-transfer 
results are not required. 

A flow diagram of the heat-transfer subprogram is given in 

Figure 4. 

Output of the Computer Program 

A general description of the information printed out by the 
computer program follows. A sample printout is given in Volume 11, Appendix || | . 

input and Geometrical Data 

The following information is printed out by the control subprogram. 
Part of the data are taken directly from input; the remainder are generated 
by geometrical manipulation of the input data. 

1. Combustor geometry at input stations. 

2. Swirler design data. 

3. Details of holes at each hole row, separately for penetration 
holes and cooling slots. 

4. At each hole row, ratios of cross-sectional areas of flame 
tube and annuli to reference area. 

5. Total hole area in flame-tube wall, for penetration holes 
and cooling slots. 

6. Ratio of total hole area, including swirler, dome holes, 
penetration holes, and cooling slots, to reference area. 

7. Identification of various key locations in the combustor: 

a. Last input station before the snout (diffuser Station 2). 

b. Last input station in diffusing part of diffuser (diffuser 
Station 3). 
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c. First hole row in the flame tube, as distinct from thf; 
dome (upstream edge is diffuser Station 4 and the start 
of the combustor annulus). 

d. Secondary hole row (upstream edge is end of primary zone). 

8. Fraction of secondary air recirculating into primary zone 
(if specified as input). 

9. Fraction of inlet air bled, and location. 

10. Fuel data: 

a. Overall fuel-air ratio. 

b. Lower calorific value. 

c. Hydrogen-to-carbon ratio. 

d. Axial distribution of fuel burning rate. 

11. Heat-transfer input data: 

a. Thickness, thermal conduct ivety, emissivity, and 
absorptivity of flame-tube walls. 

b. Emissivity of outer casing. 

12. Inlet flow conditions (at compressor exit): 

a. Total temperature. 

b. Total pressure. 

c. Mass- flow rate. 

d. Boundary- layer blockage and shape factors. 

e. Veloci ty prof i le. 

13. First estimate of downstream blockage in diffuser. 

14. Overall fuel-air ratio. 

15. Fuel burning rate at calculation points. 

16. Casing temperature at calculation points. 

17. Details of program options: 
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a. Diffuser options for various diffusing passages. 

b. Jet-mixing model and entrainment constants for penetration 
and wall jets. 

c. Heat transfer to annulus air (considered, or not considered). 

d. Flame luminosity correlation. 

e. Heat-transfer models for radiation and wall cooling. 

Output of Diffuser Subprogram 

The following quantities are printed out by the diffuser sub- 
program: 

1 . Ref e rence cond i t i ons : 

a. Reference area. 

b. Reference velocity. 

c. Inlet Mach number. 

d. Reference Mach number. 

e. Reference dynamic pressure. 

2. Diffuser parameters in the three diffusing passages, (1) 
compressor exit to lip of snout), (2) inner passage between 
snout and casing, and (3) outer passage between snout and 
casing: 

a. Ideal pressure-recovery coefficient. 

b. Actual pres sure- recovery coefficient. 

c. Diffuser effectiveness. 

d. Fractional total -pressure loss. 

e. Velocity profile at exit of the first (1) diffusing 
passage (diffuser Station 2). 

f. Boundary- layer displacement thickness and shape factor 
at inlet to the passage. 
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g. Input point at which separation occurs, if any. 

h. Mismatch area ratio at the snout. 

Output of Air-Flow Subprogram 

The following results of air-flow calculations are printed out 
for each calculation station in the flame tube and the inner and outer 
annul i : 

1. Total temperature. 

2. Total and static pressures. 

3. Bulk velocities. 

4. Mach numbers. 

5. Accumulated pressure loss due to friction and heat addition 
(annul i on ly) . 

6. Accumulated pressure loss due to expansion (annuli only). 

7. Rate of fuel burning (flame tube). 

8. Friction factor (annuli only). 

In addition, the following total -pressure-loss factors are printed out 
for each annulus: 

1. Combustor total -pressure loss relative to reference dynamic 
pressure. 

2. Combustor total -pressure loss relative to compressor delivery 
pressure. 

3. Expansion total -pressure loss relative to compressor delivery 
pressure . 

4. Total -pressure-loss factor due to friction and heat addition, 
relative to compressor delivery pressure. 

Also, the overall total -pressure loss for the combined diffuser and 
combustor, relative to the compressor delivery pressure, is printed out. 
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The following flow quantities for holes, calculated in the air* 
flow program, are printed out for each hole-row station: 

1. Pressure-loss factor. 

2. Discharge coefficient. 

3. Effective hole area (area x discharge coefficient). 

4. initial jet angle. 

5. initial jet velocity. 

6 . Fraction of current annulus air flowing through current hole 
row. 

7. Accumulated fraction of inlet air flow in flame tube. 
Finally, the air mass-flow split is printed out* 


Output of Heat-Transfer Subprogram 

The heat-transfer calculations are performed for each calculation 
point; heat-transfer quantities are, therefore, available as functions of 
axial distance along the flame tube. In addition to the flame-tube wall 
temperatures, which are the quantities of principal interest, the following 
intermediate quantities are printed out in the computer results: 

1. Adiabatic-wall temperature. 

2. Fi lm-cool ing effectiveness. 

3. Flame intensity. 

4. Flame emissivity. 

5. Heat transferred to transpiration air in the wall, where 
applicable, 

6. Radiation interchange between each wall, where applicable. 

7. Heat-transfer rates to and from the flame-tube wal Is by 



radiation and convection. 

Reynolds, Prandtl, and Nusseit numbers of the flow in the 
flame- tube and in the annuli* 
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DIFFUSER SUBPROGRAM - ANALYTICAL METHODS AND PROGRAM DEVELOPMENT 
Introduction 

The portion of a combustor treated by the diffuser subprogram 
is indicated in the sketches on the following page, for the two configur- 
ations which can be treated by the subprogram. The objectives of the 
diffuser subprogram are: 

1. To provide the necessary inlet flow conditions to the annuli 
(Stations 4A and 4B); static pressure and temperature are 
the specific properties calculated, although all others can 
be derived from these two and the mass flow. 

2. To determine the diffuser performance parameters of effective- 
ness, ideal pressure- recovery coefficient, and actual press*. ' H- 

recovery coefficient for the sections 1-2, 2-3A, 2-3B, 2-4A, and 
2-4B . 

3. To determine whether separation occurs in the diffuser and, 
if so, at what points. 

4. To evaluate the magnitude of the flow mismatch between 
Stations 2 and 2‘ (for diffusers with a snout). 

The specific input quantities required by the diffuser subprogram and the 

o 

output quantities provided have been previously listed in the OVERALL 
APPROACH sections, and will not be repeated here. 

This chapter presents the options available in the subprogram 
to accomplish these objectives, and the assumptions made and equations 
used in the analyses. 

Major Assumptions 

The assumptions employed in the diffuser analysis which are 
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common to all calculation optipns are: 

1. The flow in the upstream part of the diffuser (1-2) is 
unaffected by the flow downstream of Station 2. An analysis 
based on any other assumption is beyond the present state 

of the art. This assumption is expected to be valid in well- 
designed diffusers (see 2 below). 

2. If there is a snout, it is reasonably well matched to the 
flow; that is, the streamlines in the portion 1-2 which 
bound the flow required by the annular passages should have 
radial locations at Station 2 nearly equal to those of the 



flow in the upstream part of the diffuser will be affected 
by flow in the downstream part, in which case the results 
of the diffuser subprogram are highly questionable. 

3. If there is a snout, its performance as a diffusing passage 
is specified by an input value of the total -pressure-loss 
coefficient (i.e., ratio of total -pressure loss between 
Station 21' and the dome to the dynamic head at 2s ' ) • This 
is based on the presumption that the flow in the snout will 
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be dominated by either the colander (or orifice plate) or 
the downstream resistances supplied by the swirler and 
dome holes; in either case, it is not expected that the total - 
pressure loss due only to diffusion in the snout will be 
of sufficient magnitude to warrant a separate analysis. 

4. If there is no snout, the performance of the passages 2-4A and 
2-4B is specified by a single input value of total -pressure- 
loss coefficient between Stations 2 and 4 (i.e., ratio of 
total -pressure loss between Stations 2 and 4 to the dynamic 
head at Station 2). This is based on the presumption that 

in such cases the resemblance of these passages to a diffuser 
will be remote. 

5. The static pressure across any diffusing passage is constant 
across any section normal to the passage center-line. This 
assumes that the flow will be subjected to mild curvatures, 
and it seems to yield reasonable results for we 1 1 -des i gned 
diffusers. 

6. The air in the diffuser may be treated as a perfect gas, with 
y = 1.4 and c ®0.24 Btu per lbm-deg R. 

r 

Other assumptions which relate only to specific calculational 
options in the diffuser subprogram will be' discussed subsequently. 

Definition of Terms 

To facilitate the subsequent discussion of analysis methods and 
calculational procedures employed in the diffuser subprogram, it is convenient 
to define some terms which will be frequently used. These terms fall 
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generally into two categories: flow parameters and diffuser-performance 

parameters. 


Flow Parameters 

The flow parameters which are most frequently used are the area- 
average velocity, the mass-average velocity, the ratio of equivalent free- 
flow area to total flow area, or blockage, and the dynamic head. All 
of these parameters refer to the flow at a given axial location, where 
the static pressure is assumed constant across any section normal to the 
passage center- line. 

The area -average velocity, t7 a , is defined by 


- _JL r A 
u a “ A J 


udA 


( 1 ) 


where u is the velocity component in the direction of the passage center- 
line (normal components are neglected) and A is the passage cross-sectional 
area normal to the passage center-line. For incompressible flow, u is 

d 

the velocity the fiow would have if the profile were uniform across the 
entire passage. 


The mass-average velocity, u^, is defined by 
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For incompressible flow, it is the velocity which when multiplied by the 
mass flow yields the total momentum flux. 

The blockage is 1-E, where E is defined by 


E « 
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m 


IL 

£; 


-2 
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For incompressible flow, E is the ratio of the area which would be 
occupied by a uniform flow, with a total pressure equal to the mass- 
averagtd total pressure of the actual flow, to the actual area of the 
passage. A related term which is frequently used is the profile 
parameter, (3 , which is merely the reciprocal of E. A derivative of 
the general inlet blockage, 1-E, is the inlet boundary-layer blockage, 
l-E^p where E^ defined by 


(4) 


where is the area occupied by the boundary-layer displacement thickness. 
If the flow outside the boundary layer is uniform (i.e., u = constant across 
the passage) then E^ = E. 

The dynamic head based on mass-average velocity, q, is defined 


by 
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where p is the average ~dens i ty across the passage 


m 
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m 

For incompressible flow, q is the actual mass-average dynamic pressure; 
for compressible flow, it is treated as a reference quantity. 


P 


(5) 
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Diffuser Performance Parameters 
The diffuser performance parameters used herein are the ideal 
incompressible-flow pressure-recovery coefficient, the actual pressure-recovery 
coefficient, and the diffuser effectiveness. 

The ideal pressure-recovery coefficient is defined by 
Maximum obtainable value of (p^ - Pj) 

^p , i dea 1 — 


( 7 ) 
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where the subscripts 1 and 2 refer to inlet and exit, respectively. Two 
specific forms of this relation are used. First, the ideal incompressible- 


flow pressure-recovery coefficient, C . , defined by 


C . = 1 - 
P' 


1 

AR' 


( 8 ) 


where AR is the passage area ratio, A^/A^ . This expression in fact represents 
the maximum static-pressure rise which can be achieved in incompressible 
flow if the velocity profile at the diffuser inlet is uniform. If, however, 
the inlet velocity profile is non-uniform an additional increase in static 
pressure may be expected due to the mixing process. It is shown in Appendix 
I that the maximum achievable pressure rise is obtained by instantaneous 
mixing of the profile followed by diffusion, with the result that the ideal 
pressure- recovery coefficient in the presence of mixing can be expressed 


as 
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Although both Equations 8 and 9 are strictly meaningful physically only 
for incompressible flow, they are used here as reference quantities for 
compressible flow. 

The actual pressure-recovery coefficient, C , is defined by 
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The diffuser effectiveness,^ , is defined as the ratio of the 

actual pressure-recovery coefficient to the ideal pressure-recovery 

coefficient. Two definitions are used here, corresponding to the 

pressure-recovery coefficients defined by Equations 8 and 9: 

C 

C . 
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Alternative Analysis Procedures Available in the Subprogram 

For the purposes of analysis, the diffuser is considered as 
three separate diffusing passages: 1-2, 2-3A, and 2-3B; and two mixing 

passages: 3A-4A and 3B-4B (see sketch on page 32). The essential element 

of the diffuser analysis is then the evaluation of the performance of a 
single diffusing passage, including the determination of both the performance 
parameters and the outlet conditions. The diffuser subprogram contains 
the following three options for determining the performance of a single 
diffusing passage (the specific passages for which they can be used are 
indicated in parentheses): 


Option 1 
Option 2 
Option 3 


Empi r 5 cal -Data Method (1-2, 2-3A, 2-3B) 
Streamtube Method (1-2, 2-3A, 2-3B) 
Mixing-Equation Method (2-4A, 2-4 b) 


Option 1 has two variations: the empirical performance data can be obtained 

from empirical tabulations included within the program or it may be specified 
directly as input. Option 3 is always used in mixing passages 3A-4A and 
3B-4B, regardless of the option(s) used up to Stations 3A and 3B. 

The following three sections present the analyses used in these 
options; subsequent sections deal with the overall calculational procedure 
used in the diffuser subprogram to integrate these individual analyses. 


Empi r leal -Data Method 


I 


In addition to the general assumptions listed previously* the 
following assumptions are made in evaluating diffuser performance with 
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empirical data: 

1. The velocity profile at any section can be characterized 
as a "top-hat" profile; specifically, for a particular 
value of blockage (l-E), the profile is given by 


u = 0 for -4p->Y and Y>-^y^ 

('3) 

ti = u for -1~-<Y£-^ 
m2 2 

(14) 


where Y is the nondimens i onal annulus height. A "top-hat" 
profile is shown in the sketch below. 



2. The total pressure is conserved in the flow. This results 
in the diffuser effectiveness being reflected solely by 

a change in blockage from inlet to outlet; the tota l -pressure 
loss is assumed to be associated with the sub ^ "mixing 
out" of the prof i 1 e . 

3. The effect of the inlet Mach number on effective 

ness is negligible. This assumption is supported by the work 
of Henry, Wood, and Wilbur (Ref 14), Van Dewoestine and Fox 
(Ref 15), and Young and Green (Ref 16). Typical results from 
Reference 16 are shown in Figure 5; it is evident that Mach- 
number effects are not appreciable for M^ 0.6. 



Input Information Required 


The input information required to assess the performance of 
a diffusing passage is as follows: 

1 . The inlet condi tions defined by m, p^ , , and Ej (the 

subscript 1 refers to the inlet of the particular diffusing 
passage and is not necessarily coincident with diffuser Station 1). 

2. The nond i mens i ona 1 diffuser geometry defined by the ratio 
of the length along the passage center-line to the passage 
height at inlet (L/A r^ for annular diffusers, L/W for two- 
dimensional diffusers), and the area ratio (A^/Aj). 

3- If the option to specify directly the diffuser effectiveness 
is used, then this value is of course required in place of 
the nond i mens i ona 1 length. 

Emp i r i ca 1 Correlations Used for Diffuser Performance 

Two empirical correlations are included in the subprogram, in 
the form of tables, which yield values of the diffuser effectiveness, ^ , 
as a function of area ratio and nond i mens i ona 1 length, for a fixed value 
of inlet blockage ( I -E ^ ) . One correlation is based on data obtained for 
straight-walled, two-dimensional diffusers by Reneau , Johnston, and Kline 
(Ref 7), the data on which the tabulation is based are shown graphically 
in Figure 6. The second correlation is reproduced from work on straight- 
walled annular diffusers by Sovran and Klomp (Ref 6); this is shown graphically 
in Figure 7* The latter work was performed with diffusers having a hub- 
shroud ratio in the range of 0.55-0.7, and should only be used for diffusers 
in this range. 

For the specified geometry, the effectiveness of the passage 
is found from the tabular forms of the data represented by Figures 6 an^l 
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7 , depending upon whether a two-dimensional or an annular configuration 
is specified. The resulting uncorrected value of effectiveness is applicable 
only to the inlet blockages for which the data apply [(1-Ej) = 0.015 for 
two-dimensional diffusers and (1-E^) - 0.02 for annular diffusers^j. 

To account for different inlet blockages, the following empirical 
correlation has been developed; 

£ - I = i -f [ (1 “ E 1 ) ’ cos [ ( ?" O 5 ) 

where | = value of § for Ej = = 0.85 

This correlation applies in the ranqe l > E 1 > 0.7 and 0.5 < % 0.9 

but it is mathematically unique for all feasible values of E ^ and \ . 

This correlation is compared with some experimentally determined inlet 
blockage effects (from Refs 7 and 17) in Figure 8. 

The value of diffuser effectiveness corrected for inlet blockage 
is then obtained by first substituting the uncorrected value of 4 and the 
inlet blockage Ej for which the tabulation is valid into Equation 15 to 
determine and then solving Equation 15 again for the corrected value of 
^with the known values of i, , E ^ , and the actual inlet blockage, Ej. 

if the effectiveness of the diffusing passage is a specified 
input value, then of course the empirical correlations are not used. 

Determination of Diffuser Performance Parameters 
and Outlet Conditions 

The diffuser performance parameters, in addition to the effectiveness, 
which are determined when the emp i r i ca 1 -data method is used are the ideal 
and actual pressure-recovery coefficients. The ideal pressure-recovery 
coefficient is determined from Equation 8 and the actual pressure-recovery 



42 


coefficient is determined from Equation 11. 

The outlet conditions determined in this method are the static 
pressure, the mass-average velocity, and the outlet blockage, The 

static pressure is determined from the actual pressure coefficient. 

p 2 = p, + v, 


( 16 ) 


The outlet blockage is computed on the basis that the flow is incompressible 
and that the stagnation pressure is constant; the definition of stagnation 


pressure then yields 


2 C 1 

£ = I 

2 % + Aft (1-&) 


07 ) 


The remaining properties u m 2 ^ at the diffuser exit are determined 
from a simultaneous solution of the relations: 

m = Pi Hn2 E 2 A 2 

_ 2 
u o 

T = T = T + 
o2 1 01 '2 2g J C 

o p 
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Pi R T 2 


( 20 ) 


Limitations of the Method 

The major limitations of the empi ri cal -data method are due to 


two factors 


1. The empirical data used. 

2. The one-dimensional treatment of the flow. 

The empirical data incorporated within the program are valid only 
for straight-walled diffusers with uniform inlet velocity profiles and low 
inlet blockage. Although some data are available on curved-wall diffusers 
(Refs 18 to 21), they are either insufficient to produce a correlation, or 
they have been obtained for diffusers in which the exit profile has been 
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allowed to mix to substantially uniform conditions. No systematic data 
is available for diffuser configurations of the actual split-passage type 
found in combustors. Hence, the extrapolation of th<* present data to 
other conditions is at best uncertain. 

The one-dimensional nature of the treatment precludes accurate 
analysis of diffusers with non-uniform inlet profiles; the effect of inlet 
nonuniformities is known to be quite important (in Reference 22 for example, 
a 1 per cent change in the annulus total pressure due to a change in 
diffuser velocity profile was calculated to produce an 18 per cent change 
in the mass flow through the dilution holes). Also, the velocity profile 
at the exit of the first portion of the diffuser (Station 2) has a large 
influence on the resulting mass-flow split in the annuli; this cannot be 
treated accurately in the present method. 

Streamtube Method 

This method involves a streamtube analysis to calculate the 
development of the velocity profile within the diffuser, as well as an 
analysis of the wall boundary layers to provide the location of separation 
points, if any. The model is indicated schematically below: 

Passage 
Center-1 i ne 

St ream 1 i ne 
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In the following subsections, the subscripts U and L will refer to the 
upper and lower walls of the passage, respectively. 

Assumptions 

The following assumptions specific to the streamtube method 

are made: 

1. The flow can be represented by a number of streamtubes in 
each of which the flow is uniform and isentropic. 

2. The boundary layer is identical to that on a flat plate 
with the same pressure gradient that exists in the diffuser, 
and its effect on the inviscid flow is to produce an effective 
displacement of the diffuser wall equal to the displacement 
thickness of the boundary layer. 

3. The streamline slope varies linearly from one wall to the 
other at any axial station. 

4. Separation of the boundary layer is assumed to occur when 

the shape factor H exceeds a critical value, H se p. The latter 

value is an input to the program. For conical diffusers, 

Carmichael and Pustintev (Ref 23) show that a value of H se p = 3.0 

is typical. For airfoil sections, values of H se p = 1.8-2. 2 

are more usual (Ref. 24). A value of H = 1,9 has been 

sep 

assumed for this program. Figure 9 shows the results of 
the present calculat iona S procedure, compared with the experi- 
mental results of Reference 7* 

5. Subsequent to boundary- layer separation, no further static- 
pressure rise occurs in the section of the diffuser being 
considered. 
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Input information Required 

The input information required for a streamtube analysis is: 

1. The geometry of the diffuser in the form of radial position 
of each wall as a function of axial distance. 

2. The following inlet conditions: 

a. Weight mean total pressure, Pj. 

b. Uniform stagnation temperature, T q j. 

c. Mass flow rate, m. 

d. A velocity profile given by: 

u ■ f < Y > 

where U = an arbitrary normalization velocity 
f(Y) = a tabulated function of jj against Y 
Y = nond i mens, i ona 1 annulus height 

e. The boundary- layer displacement thickness, § p 

f. The boundary- layer shape factor, Hp 

3. The number of streamtubes into which the flow is to be divided, 

N. 

4. An initial estimate of the boundary- layer displacement thick- 
ness on each wall at every axial location of a geometric input 
point. 

Outline of Method of Solution 
The method of solution proceeds in the following steps: 

1. The determination of the static pressure at the inlet, and 
the mass flow, stagnation pressure, and velocity in each 
streamtube. 

The calculation of the static pressure at each downstream 
location (specified by the location of the geometric input 


2. 
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points), and the resulting radial locations of, and velocities 
in, each streamtube. The first time this calculation is 
performed, the estimated values of boundary- layer blockage 
at all axial locations (supplied as input) are used; subsequent 
calculations are based on revised estimates of this blockage. 

3. The determination of the boundary- layer displacement thick- 
ness on each wall, using the static pressures and wall- 
streamtube velocities obtained in Step 2. 

4. Repeat Steps 2 and 3 until the displacement thicknesses 
determi ned .i n Step 3 are within a specified tolerance of 
those used in Step 2. 

The details of these steps are presented in the following subsections 
where, for convenience, the geometrical relations are given only for 
annular configurations. 


Inlet Condi ti ons 

The inlet conditions for each of N streamtubes are determined 
as follows (the subscript j refers to a particular streamtube): 

1. Calculate the dimensionless coordinate Yj of the midpoint 
of each (equal -area) streamtube from the relations 






( 21 ) 

( 22 ) 
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At the compressor exit (Diffuser Station 1) it is assumed that 

, /u . . 

cos ( 2 ) “ i 

2. Find u^./U for each streamtube by interpolation from the 
velocity profile (tabulated as a function of V ) provided 
as input. 

3. Calculate p c ^» a characteristic density to provide a 
first guess for p. . from 

P 


1 


P ch “ RT 


ol 


where = weight mean total pressure at inlet. 

4. If the static pressure p^ at inlet is not given, make a 

first guess at this pressure from 

-.2 


Pch 


P 1 = P 1 2g 


m 


Pch- A l 


(23) 


where A. = NA. 

1 J 


5. Calculate the first guess at U from 


U = 


m 


u . 


r u 


p , A. 
H ch j 


(24) 


J 


Using this value of U, calculate from the results of 
step 2 and the temperature, T. , from 


u j 

T i o' ' 2 C p 9 0 J 


( 25 ) 


7. Calculate the density, p. , from 


P 1 

P j = RT. 

J 


( 26 ) 
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8. Calculate U from 


U 


m 


p. .A. 
J J 



(27) 


9. If this value of U is significantly different from the 
previous value of U return to 6. If not continue. 

10. Calculate trie total pressure P. in each streamtube from 



where Y = ratio of specific heats (l.k). 

11. Compute the mass flow in each streamtube from 

th . = O. u . A . 

J r J J J 

J. 

and the critical area of each streamtube A." from the usual 

J 

compressible-flow relations (p/P is sufficient to determine 
A /A*) . 


12. If the static pressure at inlet is calculated in Step 4, 
calculate the new value of weight mean total pressure 
from 

(p.) w -4 y m. .p. • 

I new m i i 

J J J 


where m. = mass flow in the j th streamtube 
J J 

13- If this new value of total pressure is significantly 

different from the given value of total pressure, calculate 
a new value of p^ from 


(P,) 


l'new P 1 TpT) 


1 


1 7 nevi 


and return to Step 7. 


( 28 ) 


ks 


Hence, the quantities Aj , p . , u j , Aj , are known at the inlet to each 
streamtube. The inlet dynamic head is calculated from 



( 29 ) 


Determination of Streamtube Properties 
at Downstream Axial Locations 

The conditions in each streamtube at each downstream axial location 
are determined in an iterative manner by estimating the static pressure 
at the axial location, computing the resulting total area of the streamtubes, 
and repeating the procedure until the computed total flow area equals the 
actual flow area normal to the passage center- line. The actual flow area is 
determined from the effective wall geometry obtained from the actual wall 
position and the estimate of boundary- layer displacement thicknesses. 

The details of the procedure are as follows (where the subscript 
j refers to a streamtube and the subscript i refers to an axial location 
of a geometric input point): 

1. As a first estimate, assume that the relative mid-streamtubc 
streamline slope at i+1 is equal to that at i: 


where 


< 7 rj > 

rj i+i 


<7r j> 


r _ f 7 l + 7u 

‘'rj -'j 2 


( 30 ) 


■ mid-streamtube streamline slope 
For the first axial location downstream of the inlet, assume 

<V i+1 -[/rL^^rU-yrL)],,, 

(since the relative streamline slope at the inlet is always 
assumed to be zero). 
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2. Estimate the pressure at the next axial location (i+1) 
from 

4 

J I 

and P i + I * P; + Ap 


TA. 


i + 1 


N 

z 

J-l 


V 2 Pj 9p 


n»j/2 Pj 9 q A? - AP 


(3D 


where TA. 

i 


= effective total area at axial Station i 


= TK 


cos Y .. 

r u' 4 » 'Cosy^+y^ ■ ( VA 


cos y r L 

cos y L + /u 





3. Calculate the stat ic-to-total pressure ratio for each stream- 

tube, p..,/P., and determine T./T ,, u. , and A. at i+1 from 
r i +1 i j ol j j 

the usual compressible flow relations. 

4. Compute total flow area based on initial estimate of stream- 
1 i ne slope from 

N 

(TA)' -(V A. cos Y .) (32) 

1+1 , J rJ i+1 

J=1 


5. 


Reestimate the relative slope of each mid-streamtube streamline 
from 


where 



and k is a dummy i ndex . 

6. Compute revised total flow area based on new estimate of 
streamline slope from Equation 33* Repeat Step 5 until 

i 

two consecutive values of (TA). + j are equal. 
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7. If (TA). + j does not equal TA. + j obtained from Equation 31. 
then repeat entire process from Step 2 with p., Pj j . A j . , 
f- j 'n Equation 30 being replaced by the current values of 

p i+r.^j,i+r A j,i+r^fj,i+r 

8. Continue the preceding process until p^ u j2 

are obtained at the exit station of the diffuser. 


Determination of Boundary-Layer 
Displacement Thickness 

With the velocity, u. f and density,^., in the bounding stream- 
tubes (j=i and j=N) known, a new estimate of the boundary- layer displace- 
ment thickness along each wall is made in the following way: 

1. The momentum thickness at each axial station i is computed 
from the following relation, derived in Appendix It from 


6 


j\i+l 


the momentum integral equation for boundary layers: 
~ 7/6 25/6 . 7/6 1/6 

e, , <" 


u j_j j ( Ai - ■ 0 - 0 °7&v 


6/7 


! , I U. . .. 

j , 1+1 


•) 


/j,i+l „ 25 / 6 . 7/6 

u j . i+1 “j » 1 


( A 

\ +] A 

X X J 


7/6 


u? dx* 
J 


(34) 


where j=l or M (referring to L and U walls, respectively) 
and 


x : 


i+1 


- x l - < x i+ i - x i ) 


cos 


X 


2 ^. .<* 


ii *71 , I 

4 'i+1 


(35) 


2. The shape factor, Hi of the boundary layer along each wall 

is computed from the empirical relationship of Dussord (Ref 25): 


H. . .. = H. . + 70(H. . - 1.05) 
j , i+l j , i J , i 


ill 

dx 


(36) 


Hj is never permitted to be less than 1.1, since this relation- 
ship is not valid for highly accelerating flows, and since 



the characteristics .of the boundary layer when subjected 
to highly favorable gradients are of little interest ijn 
the present application. Hj is also never permitted to 
exceed 3*5; this value exceeds any reasonable separation 
value. 

3* For all axial locations at which H does not exceed H on 

sep 

either wall, calculate the resulting displacement thickness 
along each wall from 

6'j.i ’ h j.i d j . i 


For axial locations where H>H _ on either wall, proceed to 

sep 

Step 5* 

4. Since this value of displacement thickness tends to interact 
substantially with the bounding streamtube velocity, a new 
estimate of the displacement thickness is determined from 
the following relation: 

where the subscript o refers to the initial estimates supplied 
as input to the boundary -layer calculation. & J j an< * S ^ j 
are determined from the solution to the following equations 
(derived in Appendix III): 

a" ( | + f|,|fN,| ’fl.l] ’•’fl.l « , 0|,| +< *‘oN,l ) 
4>n,i <5i\. ♦<%.! + »3 li.i ■ 6n.i p-tB.O 


where 




& i.i 
V r A 


240 


H 


J.t 


*1 * *M 
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M • * 1 • 05 A Q • • 1 

- 70 7 1 

’ 1 H.. x . - x . , 

j , 1 1 i-l 


II 


II 


Afi = .3 ^_1A-1 tjLizi ^oU-l JoMJ-l 

Aa j,i-l 3 *' d j,i-l 

U , i - 1 ” L, i-l 

If no separation occurs along either wall at any axial location, 
proceed to Step 7 - 

5. If H>H on either wall at an axial location, set H = 3-5 
sep 

on the separated wall for reference purposes, and calculate 
the displacement thickness along the separation wall from, 
for example: 

V I (TA) cosy>7~^ 

r cos ,-4 cosf „) S S B — L, .» 

U J c 1 Jr U . Jf 


(r L c°sy c| +rf,co5 ; frL ) 


where 


and 


(TA) = effective flow area at the separation point 
sep r 


=7T 


(r u cosf cr & N cosf rU )\ (r L cosf cl +S l cosf rL) $ 

cosf cl cos (f rL +/ r(J ) > ep 


J 


y U + JT L 


C 1 


This assumes that the static pressure and hence the flow 

area remains constant after separation. 

In the separated flow region, estimate the new values of 

displacement thickness from 

I 


/A. • = (6 • • + rf • •) 

°j,i 2 "-'oj,i °j,r 


If the newly estimated values of S. . are not equal to the 

J » 1 

originally estimated ones, <S . . , return to the streamtube 
calculation with the newly estimated values. Repeat both 
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streamtube and boundary- tayer calculations until agreement 
is obtained. 


Determination of Diffuser Performance Parameters 
and Outlet Conditions 


The diffuser performance parameters of ideal pressure-recovery 
coefficient, actual pressure-recovery coefficient, and effectiveness are 
calculated in the following way: 


C . = 1 
P 1 


-ftf 


C = 
P 


(38) 

(39) 


P 


(40) 


where q, is the inlet dynamic head obtained from Equation 29. 

The outlet conditions of are directly available 

from the procedure, as well asp., u. t m. , A. for each streamtube. The 

' J J J J 

blockage due to displacement thickness is also computed from 


r I r i i + <r,, • cosf rU 


E b. 1 1 - 


. , cos ^ rU 2\ 2 

U2 U2 cos y c i 2 j “ l 42 ° L2 cos J* 


cl2 


2 2 
r U2 " r L2 


(41) 


In the computer program, only u^r), P^, p(x . ) , S (x . ) , and H(x.) are printed 
as output information; the remainder of the quantities are used internally 
in transferring from one diffusing passage to another. 


Limitations of the Method 

The major limitations of the method are as follows: 

!• Diffusers which are highly curved or in which streamline 

curvature effects are otherwise important (such as in mildly 
stalled diffusers) cannot be treated. In these cases, the 
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static-pressure gradients across the passage are significant. 

2. The results of the method tend to be quite sensitive to the 
inlet boundary-layer properties; an accurate estimate Of 
these properties at the diffuser inlet is difficult in 
practical situations. (A method for estimating the inlet 
shape factor is outlined in Appendix IV.) 

3. The determination of the location of separation is subject 
to error. This is an age-old fluid-mechanics problem which 
remains beyond the reach of the present technological era. 

The method does, however, represent a substantial improvement over the 
empi r i cal -data method, and should be capable of treating a large number 
of diffusers accurately, particularly as experience with the use of the 
method increases. 

Mixing-Equation Method 

The mixing-equation method can only be applied to the diffusing 
passages 2-4A and 2-4B, and represents a very crude model of the flow in 
these diffusing passages. 

Assumptions 

The assumptions employed in the method are: 

1. The losses in the diffusing process can be represented as 
the sum of mixing losses, losses due to curvature of the 
passage, and expansion or contraction losses at the snout 
lip. 

2. The mixing losses can be evaluated from the determination 
of the pressure-recovery coefficient for incompressible 
flow, wherein the variation of the profile parameter, yS, 


with passage area is of the form 


= constant 


fi - K 

where K > /3^ . This is merely a qualitative indication 
that the mixing process is delayed by diffusion. 


Input Information Required 

The information required by this method is: 

1. The geometry of the passage, in terms of inlet area, exit 
area, the total angle difference between the slope of the 
passage center-line at inlet and that at exit, and the 
mean radius of curvature of the passage center- line. 

2. The inlet conditions defined by m, p^, q^, and^; the 
subscript 1 here refers to conditions determined at diffuser 
Station 2. 

3. The ratio of the area at diffuser Station 2 occupied by the 
mass flow in the annulus under consideration to the annulus 
area norma! to the engine center-line at diffuser Station 2 * , 


A^/A^,. That is, for example, 


(A thL - (r„) 2 


2A 


th- 


2A 


2A 


(42) 


where (r .) is determined from 
th 2A 


■<r. u ) 


th' 


2A 


m A 


2 7T rp udr 


(43) 


2A 


where m^ = mass flow in inner annulus 

= radius of inner casing at Station 2A 
p , u = properties from diffuser solution (Station 2) 
A similar calculation may be performed for the outer annulus. 
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2. Calculate the diffuser performance between Stations 1 and 
2, and the outlet conditions at 2, by either the empirical- 
data or streamtube method. 

3. Determine inlet conditions for the second stage of diffusion 
in a form suitable for the diffuser analysis method to be 
used between Stations 2-3A and 2-3B, and the inlet conditions 
to the snout (if any). 

4. Determine the stagnation pressure on the dome. 

5. If the diffuser has no snout, determine the inlet conditions 
to the annuli by assuming an isentropic expansion around the 
dome. 

6. If the diffuser has a snout, calculate the diffuser performance 
between Stations 2-4A and 2-4B and the outlet conditions at 

4A and 4B, by using the empi ri cal -data method or the stream- 
tube method between Stations 2-3A and 2-3B and a mixing 
analysis between Stations 3A-4A and 3B-4B, or by using the 
mixing equation method between Stations 2-4A and 2-4B. 

The following sections indicate the significant details of these 
steps, for the various combinations of analysis methods which can be used. 
These are five such conbi nations for diffusers with snouts: 


Stations 1-2 

Stations 2-4 

Streamtube 

Streamtube 

Streamtube 

Empirical data 

Streamtube 

Mixing equation 

Empirical data 

Empirical data 

Empirical data 

Mixing equation 
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It is pointed out that when the streamtube or empi ri cal -data methods are 
used between Stations 2-4, it is always implied that the methods are used 
between Stations 2-3 followed by a mixing analysis between Stations 3 and 4. 

For diffusers without snouts, the analysis of the r 1ow between 
Stations 2-4 is independent of the option employed between Stations 1-2; 
this is indicated subsequently. 

Determination of Inlet Conditions, Station 1 

The inlet conditions to the diffuser which are available from 
the program input are: 

1 . Mass flow, m. 

2. Weight-mean total pressure, P j. 

3. Total temperature, T p 

4. A nondimensional velocity profile, u/U. 

5- The shape factor of the boundary layer, . 

6. The inlet blockage, 1-Ej, the empirical-data method is 

to be used, or the fractional area occupied by the boundary- 
layer displacement thickness, l-E^, if the streamtube method 
is to be used . 

7. The fraction of 1-E, , which is on the inner wall. 

b I 

If the streamtube method is used between Stations 1 and 2, this information 
is in the proper form to enter the method directly. 

If the empi ri cal -data method is used between Stations 1 and 2, 
it is necessary to determine the static pressure and the inlet dynamic head 
from the above information. The static pressure is determined by an 
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iterative solution of the following relations: 

2 


T * T + 
ol 'l 


“ml 


2g jc 
3 o p 


P, -J5| R T, 


* Vl“ml E l A l 


!i - flUM 


i 


T ol 


which also yields anc * hence q^ 


Diffuser Performance Between Stations 1 and 2 
This is accomplished by either the empi ri cal -data method or 
the streamtube method, as previously described. A summary of the input 
and calculated quantities is given below: 


Streamtube Method Empi ri cal -Data Method 


[Input 

Calculated 

Input 

Calculatec 

m 

P2,p(xj)* 

m 

p 2 

_p . 

p 2 (r) 

pi 

E 2 

. 

mmm 


T ol 

u 2 (r) 

"l 

fit 

Uj(r)/U 
E b» 1 


E . 

~ * 
u m 2 

\ 

H i 

E bl2»^ x i^ H ^i^ 
* 

c p. 

C * 

P 

4* 


c * 

s* 



/ 


62 

The asterisks in the above table denote quantities which are printed as 
output data by the program. The subscript i refers to axial locations 
of geometric input points. 


Inlet Flow Conditions for Second 
Stage of Diffusion 

For diffusers without snouts, the only additional information 
required is the mass flow split, m fl , m B , m, . This is supplied as input 
from th= air-flow subprogram. 

For diffusers with snouts, the information required in addition 
to the mass flow split, consists of the flow properties A^, u. a ,£(or E) 
to be associated with the flows into the inner annulus, the outer annulus, 
and the snout. These are obtained somewhat differently, depending upon 
whether the empirical-data method or the streamtube method is used between 
Stations 1 and 2. 

If the empi ri cal -data method is used between Stations 1 and 2, 
these quantities are determined as follows: 

1. The nond imens i ona 1 annulus heights at Station 2 which bound 
the flow into the two annuli are calculated from (see sketch): 

(50 


y a - 


1-E 2 rn 

£_ + -J2 c 

2 m L 2 


1 -E, 


m 


y b • 1 - < 


E 2> 


(52) 



2. The areas of the passage at Station 2 which are occupied by 
the mass flows m^, m^j and m^ are calculated from; 
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A th,2A * A 2 Y A 

A th,2B “ A 2 ^ ” V 

A th ,2s = A 2 A th ,2A " A th,2B 


(53) 

(54) 

(55) 


3. The inlet blockages associated with the flows into the two 
annuli are calculated from: 

E /- S 

& 2a ” m Y a 

i _ m B E 2 


’2B m 1 -Y. 


(56) 

(57) 


t i 


t j 


and hence 


-1 


4. 


^2A,2B = E 2A,2B 

It is assumed that yS 2s = E 2s = 1. 

i 

The area-average velocities are determined from; 

= E, 


u* 


a,2A 2A u m2 


u a,2B = E 2B u m2 


(58) 


(59) 

(60) 


u a,2s =u m2 < 6l > 

If the streamtube method is used between Stations i and 2, 

these same quant i t ies are determined as follows: 

1. The area of che passage at Station 2, A^ 2A and A^ which 

are occupied by the mass flows and m^ are determinded from 


Equations 42 and 43 applied to both mass flows. Then A 
is determined from Equation 55» 

2. The inlet area average velocities are determined from 

'A 


th , 2s 


th,2A 

U = 7 \ udA 


1 


a, 2A A 


th,2A Jo 


(62) 


and a similar equation for tr 


a ,2B* 


• The inlet velocity profile parameters are calculated from 


3 



( 63 ) 
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A 


2A 


- f tH * 

Jo u 


2A 

2 


~r 


A th,2A u a, 2 A 


4 . 


and a similar equation for $2B* A 9 a,n » ' s assumed that 
*2s “ '• 

if the empi r leal-data method is to be used between Stations 

2 and 3. the inlet blockages associated with the flows in 

* 

the two annuli are calculated from: 

& 2k 


E 2A " 


^2A +< ^2B 


*b| ,2 


t . <hi E 

' 2B *2A + <5zb 


bl ,2 


( 64 ) 


( 65 ) 


For the purposes of the diffuser analysis between Stations 2 and 
3, the areas of the diffusing passages at Station 2' are assumed to be 


A th,2A'* A th,2A COS «if C,A ^ 

and 

A th,2B'“ A th,2B C0S ^C,B ^ 

where is the angle of the slope of the passage center-line relative 
to the axis of the combustor. This assumes in effect that the inlet 
areas to the annular diffusing passages are the areas occupied by the 
respective flows at Station 2, projected normal to the passage center- 
line. It is obvious that this assumption is not valid when the areas 
A th 2A' an< * A th 2B* substantially from the actual geometrical areas 

of the diffusing passages at Station 2' (A ' ^ cos^). The output of the 
computer program includes the ratios k^ 2B ,7A 2B* A th 2A'^ A 2A' ( < * #note< * 
as "area mismatch at the snout" in the program). 

The inlet dynamic heads at the diffusing passages are calculated 
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from: 


'2A 2g 


J_ 02A u a.2A 


( 68 ) 


o A 


th,2A' 


and a similar equation for q 2 g. The inlet dynamic head at the snout is 


determined from: 


i A u „ 
— _ 1 s m2 


(69) 


2s 2g 


o A 


th,2s 


Conditions on the Dome 

If the diffuser has no snout, the stagnation pressure on the 
dome is assumed to be constant and is determined by 


where 


P d - P 2 + 0 - * d > q 2 
^2 “ Tg Pi “m2 


(70) 

(71) 


The coefficient is an input quantity representing the number of dynamic 
heads lost in the mixing and flow processes occurring around the dome. 

The magnitude of will depend upon the shape of the dome and the 
effectiveness of the diffuser between Stations 1 and 2, and is difficult 
to determine a priori. It is expected that its value would be in the 
range 0.2 - 0.5. 

For diffusers with snouts, the static pressure just inside the 
snout, P 2s , , is determined by assuming a sudden contraction or expansion 
process between the areas A^ 2j and A 2s , (the actual geometric area of 
the snout entrance); this is accomplished by the application of Equation 
44 or 45. The dynamic head just inside the snout is calculated from 

q 2s' ( A 2s , ' q 2s 

The total pressure on the dome is then calculated from 
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P d " «>2s‘ + (1 ‘ K d> ^2s' (72) 

where Is an Input quantity representing the number of dynamic heads 

lost in the flow through the snout. 

Diffuser Exit Conditions If 
No Sr.out is Present 

if no snout is present, the diffuser exit conditions at Stations 
4A and 4B (the entrances to the annuli) are calculated by assuming that 
the flow expands isentropical ly from the dome stagnation pressure determined 
by Equation 70. For either annulus, the conditions at Station 4 are 
completely determined by (or A B ) , T Q , P rf , and A^ A (or A^ 0 ). 

Diffuser Performance Between Stations 2 
and 4 for Diffusers with Snouts 

The procedure for determining the performance of the diffusing 
passages between Stations 2 and 4 depends upon the diffuser-analysis 
method employed. The procedures are indicated below. 

Mixing-Equation Method 

\ 

In this case, the procedure is identical do that described 
previously. The inputs required for the analysis are i& A , p 2> q 2 ^, ^A* \h 2A' 
and a similar set for the outer annulus; these are determined in the manner 
indicated in the previous paragraphs. The calculated quantities are 

C pmA» C pA' $mA* P 4A’ T 4A» P 4A and a sirn,,ar set for the outer annu!us * 

All of these quantities are printed out by the computer program. 

Streamtube Method 

In this method, three steps are Involved since the basic stream- 
tube analysis is only applicable between Stations 2 and 3» end does not 
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consider any losses due to expansion or contraction of the flow at the 
snout lip. The steps are as follows: 

1. A basic streamtube analysis, as previously described, is 
performed between Stations 2A-3A and 2B-3B. The inputs 
required for this analysis are * A , A th 2A , , p 2 , ^(O.c^a* 
H 2A and a similar set for the outer annulus; in addition 

it is assumed that the boundary- layer properties on the 
snout lip are S- 0 and H = 1.4. The only calculated 
quantities which are used subsequently are the static 
pressures, p, A and p,_, and the profile parameters, and 

jA }u 3A 

which are computed from Equation 63 applied to Stations 
3A and 3B. The diffuser performance parameters C^, C^. , and 
as obtained from this analysis are printed as program 
output, however. 

2. It is assumed that the profile mixes to a uniform one 
between Stations 3 and 4. The static pressure at Station 
4 in the absence of area mismatch at the snout lip is 
accordingly obtained from Equation 5*7 in Appendix V: 


4A 


- P 


3A _ 1 


*3A 


( , + ifi . 

fi 2 1 A 4A 3 A A 4A 
3A 


) 


(7*f 


3. The effects of any area mismatch at the snout are rather 
crudely approximated by assuming that the total-pressure 
loss between Stations 2 and 4 due to this process is given 
by the sudden expansion or contraction relations (Equations 
46 and 47), and that this loss may be applied directly at 
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Station 4; that is, in the inner annulus for example: 

A. 

i i _ 

2 


P 4A P 4A 


1 / | 2A' 

“ 7 I' - T 


-)( *) 2 


( 74 ) 


zA 


th,2A' 


2A 1 


if A^,./A . OA /l. P', A is the total pressure corresponding to 
2A' th,2A' 4A 

p, 4A- V *A’ and A 4A- 

The calculated quantities are the diffuser performance parameters 

C , C and y. between Stations 2A-4A and 2B-4B, and the flow conditions 

pm p , 61 m 

P4A' p 4a' t 4A’ p 4b ’ P 4B ’ and T 4b’ A " ° f tneSe entities are part of 
the printed output of the computer program. 

Emp i r i ca 1 -Data Method 

As in the streamtube method, three steps are involved in the 


procedure: 


1. The basic empirical-data method, as previously described, 
is performed betweea Stat i ons 2A-3A and 2B-3B. The inputs 
required for this calculation are m A> A th 2A , , p 2> q 2A> E 2A * 
and a similar set for the outer annulus. The only calculated 
quantities which are used subsequentlv are the static pressures 
p^ A and p^, and the profile parameters ( = E^ A and 

( = E^ a ^ ) . The diffuser performance parameters C^, C^. , 
and ^ as obtained from the basic method are printed as 
program output. 

2. A mixing analysis is applied to the flow between Stations 
3A-4A and 3B-4B. This is indent ical to the procedure used 
for the streamtube method. 

3. The affects of any area mismatch at the snout are included 
in the manner used for the streamtube method. 


The calculated quantities and the program output are identical 
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to those of the streamtube method. 

Structure of the Subprogram 

The diffuser subprogram consists of: 

1. A subroutine which sets up starting conditions and directs 
the diffuser calculation through the other subroutines in 
the appropriate sequence. 

2. A group of subroutines which perform the streamtube analysis. 

3. A group of subroutines which perform the empirical-data 
analysi s. 

4. A group of library subroutines which provide gas properties 
and interpolation procedures. 

Detailed flow charts and the program listings for the diffuser subprogram 


are presented in Volume II. 
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AIR-FLOW SUBPROGRAM - ANALYTICAL METHODS AMD PROGRAM DEVELOPMENT 
Introduction 

The main objectives of the air-flow subprogram are: 

1. To calculate the air mass flow and pressure distribution 
in the flame tube and annuli and, hence, the overall 
pressure loss. 

2. To provide a starting point for the heat -transfer 
calculation by working out the temperature and velocity 

distributions of the gases in the flame tube and annuli. 

To facilitate analysis, flow in the combustor is broken down into 

distinct but interacting streams: 

1. The annulus air (inner and outer annuli). 

2. The air flowing into the flame tube through: 

a. Swirler. 

b. Dome holes. 

c. Penetration holes in the flame tube. 

d. Wall -coo ling slots *n the flame tube. 

e. Porous flame-tube walls. 

3. The mixture of fuel, air, and combustion products flowing 
through the primary zone (the flame tube upstream of the 
secondary holes). 

4. Parallel streams flowing through the flame tube downstream 
of the primary zone: 

a. The hot stream, a mixture of combustion products, 
unburned fuel, and air. 

b. The cold streams, those portions of the penetration 
jets and wall-cooling films, originating upstream, 
that remain unmixed at a given axial position. 
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This section of the report is devoted to a detailed description of 
the calculation methods and general organization of the air-flow subprogram. 

Assumptions 

The air-flow calculation rests on the following main assumptions, 
in addition to those discussed under OVERALL APPROACH: 

1. The primary zone can be treated as a stirred reactor: 

a. Static pressure, temperature, and fuel-air ratio 
are uniform throughout the zone. 

b. Mixing and burning occur instantaneously. 

This is one of the major simplifications in this work, and 
it precludes the possibility of calculating, for example, 
pressure distributions in this zone, it will be particularly 
invalid at low pressures, when burning in the primary zone 
is reaction-rate limited, unless an accurate estimate of 
the fuel-burning-rate distribution (see assumption 3) has 
been made. 

2. The flow can be considered to be one- d? mens i ona 1 . This 
represents a considerable simplification, particularly 
In the case of the gases in the flame tube. However, 
reasonably simple methods are not available at the present 
time for treating flow and heat transfer in a gas stream 
having non-uniform temperature and velocity profiles. 

3. The rate of fuel burning as a function of axial length 
can be specified as input. This assumption makes the 
calculations easier but a realistic estimate of the fuel- 
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burning rate will be difficult to make, particularly 
(as mentioned above) at low pressures. If the fuel -air 
ratio at any axial position, which will not be known in 
advance, exceeds the stoichiometric value, the excess 
fuel is assumed to be available for burning at the next 
downstream position. 

Calculation Procedure - Outline 

In calculating the flow conditions throughout the combustor, 
the air-flow subprogram goes through the following routine for each 
iteration on the mass-flow split: 

1. The mass-flow rate through the dome is obtained from the 
flow split, and the total pressure on the dome is obtained 
from the diffuser calculation. 

2. The pressure drop across the dome to the primary zone is 
found by equating the combined swirler and dome-hole flows 
(functions of pressure drop) to the total dome flow. 

3. With the primary-zone static pressure established, flow 
conditions in the annulus and flow though flame-tube holes 
into the primary zone are computed for calculation stations 
up to the secondary holes (the end of the primary zone). 

4. The portion of the secondary-hole flow that recirculates 
into the primary zone is calculated, either directly by the 
program or from a fraction specified as input. 

All of the contributions to primary-zone air flow (through the swirler, 
dome holes, and flame-tube holes up to the secondary holes) are now known. 
In addition, the fuel flow into the primary zone is known from input. 
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5. The temperature rise due to combustion is calculated for 
the primary zone treated as a stirred reactor. (Excess 
fuel above stoichiometric is carried downstream.) Flow 
conditions at the end of the primary zone are then found. 

For the remainder of the calculation the combustor is divided into axial 
control volumes bounded by adjacent calculation stations k and k + 1. 
Calculations proceed downstream from station to station. In general, the 
flow conditions at the upstream station, k, are known from the results of 
the previous calculation step. Conditions at the downstream station are 
found as follows: 

6. From the annulus-to-f lame-tube pressure drop at station k 
the mass flow rate through the holes just downstream of 

k is found. The mass flow remaining in the annulus is then 
known. 

7. For the annulus, the equations of momentum, energy, continuity, 
and state are solved for the pressure, temperature, velocity, 
and density at station k + I . 

8. Correlations for flow and mixing of wall and penetration 
jets are used to compute the fluxes of mass, axial momentum, 
and enthalphy contained in the jets entering the control volume 
(through holes at station k and in the form of residual jets 
originating upstream of k) and leaving the control volume 
(residual jets at k + 1). The net transfer of mass, 
momentum, and enthalpy from the jets Into the main gas stream 
in the flame tube may thus be found. 

9. The equations of momentum, energy, continuity, and state are 
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written for the hot stream in the flame tube (excluding 
residual jets) between stations k and k + 1 . Heat addition 
is included for fuel burned, up to stoichiometric; any 
excess fuel is carried downstream. The solution of these 
equations yields the temperature, pressure, density, and 
mass flow at station k + 1. 

The sections that follow describe the detailed calculations 
that are performed in each of the above steps. 


Flow Through Holes and Swirler; Swirler ?. ,siqn; Recirculation 

Flow Through Holes 

The mass flow rates through the holes in the dome and in the 
flame-tube walls are calculated from the discharge equation: 


m. - C . A, ( 2g p (P - p £ ) 
h d h 3 o K an an r ft 

where P * total pressure in annulus 
an 



an 


'an 


p u 2 /2g 
r an an 3 o 

discharge coefficient 


(75) 


Discharge Coefficients 

For dome holes a value of 0.6 is used for the discharge coefficient. 

For flame-tube holes a table of experimentally measured discharge 
coefficients for 100 hole types is supplied as library data. An index to 
these hole types is given in Table I of Volume 11. (Data for other hole 
types may be added to the library as required.) The tabulated data are 
drawn from experiments reported by Venneman (Refs 28 and 29), Marshall 
(Ref 30), Kaddah (Ref 31), and Dittrich and Graves (Ref 32 and 33)* 
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The discharge coefficient is plotted against the hole pressure- 
loss factor. 

a p h 


P an ' p ft 


I + 


'an 


'an 


(76) 


where AP h = P an - P ft 

= Static pressure in flame tube 

Pressure-Ratio Correction 

It has been shown by Dittrich and Graves (Ref 32) that the effect 
of flame-tube velocity on discharge coefficient can be predicted by: 


where 


C d - ci (0.75 - 0.25 P an /P ft ) 

= the actual discharge coefficient (which is a function 


(77) 


of pressure ratio^ 

Cj = the "corrected" discharge coefficient (which is in- 
dependent of pressure ratio) 

The values given on the library tape are corrected discharge coefficients. 

In Equation 77, p is replaced by P if the hole is a 
scoop or other total-head device. 


Swi rler i'Iow 

The equation used in the program for the flow through a swi rler 
was taken from Reference 3 , k 


A? 


sw 


'ref 


= K 


sw 


ref 2 . 

— sec p 

sw 


sw 


ref 

* ft 2 


m 


sw 


( 78 ) 


m 


J "ref 


where A = area, ft 

m = mass flow, lbm per sec 

J 

q = dynamic head, lbf per ft 
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fi ■ swi r ler angle 

A p d * total-to-static pressure drop, Ibf per ft 
sw ■ swi rler 

ref = reference condition 
ft “ flame tube 

In the program the constant K has been given a value of 1.3, which is 

SW 

suitable for stra i ght-b laded swirlers. For curved-wall swirlers, the 
value of K could be changed to 1.15. 

Swi rler Des i qn 

A comprehensive literature search brought to light virtually 
no work on which such a sound theoretical swi rler design method could be 
based. The method used is therefore based on an examination of swi rler 
designs that have previously been used on 16 combustors, mainly of the 
tubular type. The results are shown in Figure 10. Although the range of 

D r „ does not cover the values normally encountered in annular combustors, 

ft 

the following may be used until more information on swirlers is available: 


Sw i r 1 e r 0 . D . 

- 0.225 D ft 

+ 0.75 in 

Swi rler 1 .D . 

H 

O 

O 

o 

“tl 

ft 

+ 0.25 in 


No. of Blades ■ 8 

Blade angle = 50 degrees 

Rec? rculati on 

In an actual combustor a complex pattern of recirculating flow 
exists, particularly in the primary zone. Representation of recirculation 
is introduced in the computer program by a provision for partial recircula- 
tion of secondary-hole flow upstream into the primary zone. Only the flow 
through the secondary holes is treated in this way. 


77 


The fraction of secondary-hole air flowing into the primary 
zone, m /m , may be specified as input to the computer program. 
This must be set to zero if recirculation is not to be considered. If 
the reci rculat ion fraction is not specified as input, it is generated 
within the program from the following formula: 


= 0.5 sin f 
"sy 


an 


ft 


1/2 


(79) 


This is an empirical expression based chiefly on measurements obtained by 
Rosenthal (Ref 35 ) in a tubular combustor. 

The air entering the primary zone by recirculation from the 
secondary holes is treated in exactly the same way as air entering through 
swirler, dome, and wall holes. Air going downstream is treated as 
ordinary penetration-jet air. 


Primary Zone of Flame Tube 


Pressure 


The pressure in the primary zone, pj , is found from the pressure 
in the diffuser upstream of the dome, , obtained by the diffuser sub- 
program, and the pressure drop across the dome, which is a function of 
the flow through the dome and the characteristics of the swirler and the 
dome holes: 

p p = 


. 2 
"d 


^ ^o/* ref 


A dh C dh + 


1 


/ K 

n/ SW 


2 a 

sec P sw 


, A 
«- sw 


m^ * total flow through dome 

* m + m ,, 
sw dh 

riidh * flow through dome holes 


i?]J 


where 


( 80 ) 
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= discharge coefficient of dome holes 
= coefficient in swirler pressure-drop equation 
= blade angle of swirler 
= cross-sectional area of primary zone 
This equation is derived in Appendix VI. 


Temperature Rise 

The energy equation for combustion in the primary zone takes the 
f ol lowing form: 


m i h • t* 
1 a i r , T 


+ q = (*, +* f „> h T . 


( 81 ) 


where m^ = air mass-flow rate 

m^k = fuel burning rate 
q = heat release rate 

h . = enthalpy of air 

h^j = enthalpy of combustion products 
and subscripts 1 and 2 refer to conditions before and after combustion, 
respectively. Velocity terms have been neglected because mean velocities 
in the primary zone are usually quite small. 

The enthalpies of air and combustion products are calculated 
from the expression: 


■ J* 


c .dT 
P 


(82) 


where the specific heat of air as a function of temperature is correlated by 
c = 0.2419 - 0.8181 x 10" 5 T + 1 7.91 x 10* 9 T 2 - 2.743 x lo' ,2 T 3 


( 83 ) 
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The data of Fielding and Topps (Ref 3&), shown here in Figure 11 , have 
been used to modify Equation 83 for calculating the specific heat of 
combustion products: 


c p = 0.2419 + 0. 1 03 f - (0.8181 - 22. 6f) 10~ 5 T 

+ (17-91 - 29.6f) 10" 9 T 2 - (2.743 - 0. 350 10~ ,2 I 3 
where f = local fuel-air ratio. 

This equation is suitable up to stoichiometric values of the fuel-air 


( 84 ) 


ratio: 


f . = 0.0867 * 

stoi ch 


1 + H/C 


where 


1 + 3H/C 

H/C = fuel hydrogen/carbon ratio by mass 


(85) 


D i ssoc i at i on 

There are two factors in the heat balance for the primary zone 
(or between stations in the flame tube): 

1 . Burni ng of fuel. 

2. Dissociation of products of combustion. 

The equations for c p given above neglect dissociation effects; 
Hodge (Ref 37) suggests that dissociation may be represented by defining 
an effective calorific value of the fuel: 


h r r _ = h -Ah , , 

e “, p p eff, p 


where h p * fuel lower calorific value, per unit mass of fuel 
Ah = energy supplied to cause dissociation of combustion 

products 

In Figure 12, curves of Ah £r are plotted as a function of 

eff, p 

temperature, pressure, and fuel-air ratio, for combustion of a fuel in 
air. The effect of pressure is ignored in the present work and the curves 
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are approximated by: 


Ah 


3 x lo’ 26 ^* 5 


eff ,p 


where 


f = local fuel-air ratio 


Tj = temperature of products of combustion 

If the energy of dissociation is related to the enthalpy of the 
products, then the enthalpy of the products, including the effects of 
dissociation, h 1 ro ^ , is given by: 

' "prod <V * h prod <V + 3 » '°~ 2 S 7 ' 5 f a 


where f = fuel-gas ratio 
9 

_ f 
1 + f 

T b = base temperature for enthalpy 
■ *+59-7 deg R 

Thus dissociation has the effect of lowering the temperature of the 
combustion products. 

Hence for the combustion of fuel in a vitiated gas, when the 
enthalpies do not include the effects of dissociation, the heat release 
rate, q, is adjusted to account for dissociation and recombination. In 


general, for combustion between conditions 1 and 2: 

7.5 


q = m fb h p + 3 x (0 


-26 


*1 T 1 


m„ T_ 
2 2 


7-5 


1 + f 


1 T f. 


(86) 


where 


m 


fb 


= rate of burning of fuel 


For combustion in the primary zone the temperature rise is very large 
and the * ** term in Equation (86) may be ignored. Equation (81) is 
now solved for , the primary-zone temperature. 

The empirical method described above is widely used and gives 
satisfactory results in practical cycle-performance calculations. 
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Density and Velocity 

The density is found from the equation of state: 

p = .JZi 

J 1 R T, 


and the velocity from the continuity equation: 

m„ 


U„ = 


A A 


where = m^ 


2 2 


Hot Loss 


The "hot loss" of total pressure, which is the change in total 
pressure due to heat addition, is given by 


A P hot P 1 " P 2 


(87) 


Ai. „ / jl \ 

"... A ' 


Flow in the Annulus 

The equations of momentum, energy, continuity, and state are 
now set up for a section of one annulus between adjacent calculation stations, 
labeled 1 and 2 in the sketch below. 



p 2’ U 2 ’ m 2 
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The resulting equations are then used to calculate conditions at Station 
2 front the known conditions at Station 1. 


Momentum Equation 

Derivation of the equation for conservation of momentum between 
Stations I and 2 is complicated by the fact that, in the general case, 
the annulus wall between these stations contains a hole through which 
some momentum flows. To deal with this complication the calculation is 
carried out in three steps. 

First, the flow through the hole, m^ , is found from the hole 
type and size, the pressure drop across the hole, and discharge-coefficient 
data. (An expanded discussion of this topic will be found later in this 
sect i on. ) 


Second, it is assumed that the air remaining in the annulus at 
Station l', immediately downstream of the hole, undergoes a sudden 
expansion. The resulting tota 1 -pressure loss is given by the following 
empirical expression, from which P the total pressure at Station 

l 1 , may be calculated: 


P p . ,.8 5 j W . JjP. ,2 

an, 1 an, 1' 2g 


ft. 


. 36m 


an, 1 J 


l/(.5+.242M 2 ' 21 ,) 
an, I' 


(88) 


(This expression has been obtained by fitting the curves of Figure 13, 
which is taken from a review of empirical data on total -pressure losses 
due to sudden expansion in compressible flow, Reference 38.) With the 
total pressure at Station 1‘ known, the other flow conditions at Station l 1 
may be found frort. the equations for the iser cropic flow of a perfect gas 
between Stations 1 and 1‘. The most convenient approach is first to 
calculate the Mach number at Station 1‘, M an j, from the following equation: 
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A 


UhlL 


m 


an, 2 


an, I 1 L Man * U 


■B] 


an, I 


^ an, 1 


r* i 


i + 


2 

M an, 1 ' 


y + 1 

2(r-o 


(89) 


where = stagnation temperature, assumed to be constant across the 

expansion 

2 

u i 

* T . - -22*1- 
an ’ 2C g J 

P O 

The static temperature and static pressure are then determined from: 


an, 1 ' 


, , r- 1 h 2 

1 2 M an,r 


-1 


an, I 1 = 

an, 1 ' 


an, I 


y-j 


Then from the equation of state the density at I' is: 


_ P an, 1 


an, 1' RT 


an, I 1 


From continuity, the velocity is: 


an, I • 


ft 


an, i 1 


P , ,A . 
v an, 1 ' an, l 


Thus, the assumption of a sudden expansion at Station l 1 permits all of 
the necessary flow quantities at that station to be computed. 

The third step is to apply the momentum equation for flow between 
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S tat i ons 1 ' and 2 


d , . A , + m _u ,,/g+ 
an, l ' an, I an, 2 an, I 1 o 


r . 

J „ r an dx 


r * 2 2 

- FAp u / 2g . dx = 
J w an an 


p .A 0 + m _u _ / g 
^an,2 an, 2 an, 2 an, 2 


(90) 


where F - friction factor (Fanning) 

= 0.0035 + 0.264 Re"* 4 * (Ref 39) 

3 n f I 

A^ = wetted wall area per unit length 

A number of simplifications are made in this equation. Since 
the pressure losses due to wall slope and friction are small terms in the 
equation, an iteration loop can be avoided without a significant loss 
in accuracy by writing: 


f 2 F A P u 2 / 2g .dx = F (C. + C D ) p , 
w "an an 3 o A B' "an, 1 

x, . 


' u an,r (x 2 • x 1') /2 9 0 


where C. and C are the inner and outer wall areas per unit length 
A B 

between l 1 and 2. 

The integral pressure term can be evaluated if it is assumed 
that flow between Stations 1' and 2 is incompressible and that the mixing 
rate will be delayed by diffusion. With these assumptions the integral 
pressure term may be integrated as shown in Appendix VII, yielding: 


dx = (p . . - p _) 
VK an , 1 K a n, 2 ' 


2A , A 0 . 

an, 1 an, 2 + p A 


an, 2 an, 2 K an, 1 1 an, 1 


A , + A . 
an,l an, 2 


II 1 
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The final momentum equation thus becomes: 


2 A , A . 
an,l an, 2 , 

A . + A . ' P an , 1 1 
an,l an, 2 


L. 

o ) - F (C + C an >' ' (« - X ) - 

p an,2' f lL A V 2g lx 2 x l' ) 

3 o 


m _ (u - - u ,,)/q 
an, 2 ' an, 2 an, 1 7 


(90 


Energy Equation 

The energy equation for flow in the annulus between Stations I' 


and 2 is: 


2 

f h (T •) + U 2 ./2g j"l m . + P q C A .dx 
|_ an , 1 an, 1 3 o J an, 2 J x H w A 


m 


= f h (T _) + u ./ 2g jl ... 

L an, 2 an, 2 o J an, 2 

where h(T) = static enthalpy of annulus air at temperature T 

q = heat transfer rate from unit area of the flame tube wall 

w 

to the annulus ai r 

The casing temperature is normally assumed to be equal to the compressor 
delivery temperature. Even when a different casing temperature is specified, 
however, the heat-transfer rate between casing and annulus air is assumed 
to be negligibly small. 

Since the variations in annulus air temperature are normally 
small, h(T) may be assumed equal to c (T-T^) , where Cp is a mean specific 
heat. If mean values are also assumed for q and C., the simplified energy 

W M 

equation becomes: 

o 

(c T . + u ,/2g J) rh ? + C. ( x o " x i) 
v p an,l an,r 3 o an, 2 ^w v a 2 1 

= ( c T - + u 2 /2g J) m . 
v p an, 2 an, 2 3 o an, 2 


( 92 ) 
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Continuity Equation 

The continuity equation for flow in the annulus between Stations 


1 and 2 is written as follows 


A an. 1 m an.2 /m an,l “ P an,2 U an,2 A an,2 


0 , u , m . m ,/"< i 

H an,l an, l an,l an ,2 an,l 


Equation of State 


The equation of state for the annulus air is 


an, 1 


an, 2 


y^an.l ^an,l J^a n,2 ^an,2 


Method of Solution of Annulus Equations 

Since the mass flow through the holes between Stations l and 2, 

m, , is known, the downstream flow in the annulus, m 0 , is simply: 
h an, 2 r 1 

• • • • 
m „ = m , , = m , - m, 

an , 2 an , 1 an , 1 h 

There remain four unknowns, u „ , p „,n „ and T „ , which can be 

an , Z an , L / an,*., an, 2 

found by solving Equations 91. 92, 93. and 9^* 

These equations are easily reduced to a quadratic in u 

an , L 


where 


B 2 - 1 - 


B 2 u an,2 + B l 'in, 2 ‘ B o ‘ 0 
^ C p A an , I + A an , 2 


2 Jc 

O r 


an , 1 


an, 2 


an, 1 1 . 
m 


an ,2 


U l I 

an , 1 


A , + A - 
an,l an, 2 

2 g A . 

3 o an,l 


• 1 - 


f ( C A + C B> <*2 - x l'> 


an , l 
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B * u , + 2gJc 
o an,1 op 


6 C, f 

T . + H w A v 2 1' 

an , 1 


1 


c m. - 
p an, 2 


The solution of Equation 95 is 
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where physical reasoning dictates the choice of the positive root. 

W : th 11 . known. Equations 91, 93, and 94 are solved for 
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Introduction to Jet Flow and Mixing 

The flow characteristics of jets entering the flame tube through 

holes and cooling slots is treated by the use of empirical correlations 
for: 

1. The length of the jet (measured along its center line). 

2. The cross-sectional area of the spreading jet as a function 
of jet length. 

3. The jet center-line velocity and temperature as functions 
of jet length. 
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4. The development of the nondimensional ve ocity and tempera- 
ture profiles across the jet with jet length. 

Two sets of correlations are provided, one for wall jets (where the air 
enters the flame tube through a cooling slot in a direction parallel to 
the wall) and one for penetration jets. These two types of jet are 
illustrated in the sketch. 




Wal 1 Jet 


Once the jet properties have been estimated, the rate of mixing 
of the jet air into the main gas stream is estimated according to one of 
three optional mixing models. The same model, chosen by the user of the 
program, is used for both penetration and wall jets. 

The characteristics of penetration and wall jets and of the 
three mixing models are given detailed discussion in subsequent sections. 


Characteristics of Pentration Jets 


The dimensions used in describing jet behavior are shown in 
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the sketch below, where x is the direction of flow of main-stream gas in 
the flame tube, y is jet penetration distance normal to the flame- tube 
wall, s is distance along the jet centerline,^ is transverse width 
of the jet (measured normal to the jet center-line along lines of constant 
s), and the angle between the initial jet trajectory and the flame- 

tube wall at the entry hole. 



Initial-Jet-Anqle Data 

The jet-mixing correlations, and the momentum equations in the 
flame tube, make use of the initial jet angle This is available as a function 
of pressure-loss factor for 47 of the hole types in the library data. Where 
the data do not exist, the following correlation is used: 

sin 2< f=-^~ (97) 

an 

A p h 

where = the asymptotic value of as 1 + — — tends to 

° q an 

inf ini ty 

This approximation gives reasonable results for round flush holes. For 
total-head devices its accuracy is unknown; in particular, the following 
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assumption must be made: 

^ = 90 degrees when 


C d C dn 


Penetration Distance 

The penetration distance y, measured from the wall to the jet 
center-line, is calculated by an empirical correlation of experimental data 

obtained from 18 sources (Refs 40 through 57): 

0.47 0.85 0.32 

v -0.87d h( |f2) (f) sin? (98) 

ft h 

where d, * effective axial length of hole 
h 

■ (Cj)^ x axial length 

Uj Q * initial jet velocity 

* ^h 

P C . x hole area 
y an a 

It is assumed that there is no change in density of the jet as the air passes 
through the hole. Equation 98 is plotted in Figure 14. 


Jet Center-Line Velocity 

An empirical equation for the jet center-line velocity has 
been obtained by fitting carves to the experimental data of Keffer and 


a . 

J 
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5 - 2 <t 


, Fig 

11), shown 

here on Figure 15: 
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J (u J,o' u ft ) (99) 

distance along the jet 

center-1 ine, s', 

used in the above 
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equations, is calculated (Ref 56, page 492) from: 


u . 

- s + (2.93 * 0.279 — Lj ^) d. (100) 

u ft 1 


where s 

i 

j 1 + (^) .dx 

( 101 ) 

d. 

J 

= 

effective initial jet diameter' 



- 

(C,. Hole Area . 4 /Tt) 2 
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Transverse Velocity Profile 

The correlating equations for the jet velocity profile, 

expressed 

in terms of 

the 

transverse coordinate^ , are obtained from Figure 

lb (Ref 

56, Fig 9): 




0<T)<Tl£ 

u. 

J.T| 

* u ft + [o.5 cos (S3-) + 0 . 5 ] (u jft - u ft ) 


T1 >11* 

u . 

J >11 

- u ft + 1.74 e"' 2511 ^ (u Ji4 -u ft ) 

(102) 

where, from 

Figure 12 of Reference 56: Tj ^ = 0.215 s' 




Transverse Temperature Profile 



The transverse temperature profile is assumed to be similar in 
form to the transverse velocity profile. The Prandtl number is proportional 
to the relative spreading rates of velocity and temperature profiles (Ref 58). 
Since the Prandtl number for air in the jets is about 0.7 for typical annulus- 
air temperatures, the dimensionless temperature profile is assumed to be 
roughly 1.4 times as wide as the velocity profile at a given position along 
the jet. For the same reason, the center- line temperature decreases more 
rapidly than the center-line velocity. The expressions relating temperature 
to velocity are then: 
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JlS - u ft 

6.5 ' u . - u 

J-o ft 


) 


(103) 


(104) 


The constants were taken from reviews of jet data by Forstall and Shapiro 
(Ref 58) and Squire (Ref 59 ) - 


Calculation of Jet Properties 

Examination of the data of Jordinson (Ref 57), reproduced here 
in Figure 17, shows that the cross-sectional area of a jet in a plane normal 
to its center-line may be represented by: 


A. = 0.127) 
J >cs ( 


* 1.8 


where VJ * = the value of 7J at the assumed jet boundary. 

The jet area normal to the combustor axis is then given by: 

* 1.8 


A. = 0.1 27) ‘•Vcosfj 


(105) 


where 


. = local jet angle 


= cos 


Other jet properties needed for solution of the flame-tube 
equations may be calculated as follows: 

•w* 


rh. * 
J 

J.'i 

u . 0.12x1 .8*n 

J .*1 

0.8 

.di) 

(106) 
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0.8 

0. 12 x 1 .87j dYJ 

(108) 



The value of if* depends on the mixing correlation employed; this topic 
is discussed after the wall-jet correlations have been presented. 


Characteristics of Wall Jets 

A large number of correlations for wall-jet characteristics are 
available in the literature; most give essentially identical results. The 
methods described have been chosen mainly on the grounds of simplicity. 

The symbols used are defined by the sketch below. 



Max imum Veloc i ty 

The maximum velocity in the wall jet is calculated from equations 
derived by Kruka and Eskinazi (Ref 60 , Eqs 16, 17. and 18): 
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where 
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initial 
/an C d 


jet velocity 


x slot area 
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a 


| .06 

J “ u ft/ u h 



The above equations apply for values of ^ between and u^. Otherwise, 

u. .is set equal to u, . 

J ,i h 


Veloci ty Prof i le 

For purposes of analysis, the velocity profile is divided into 
two sections about its maximum. Near the wall ( < cSj ) the equation of 
Kruka and Eskinazi (Ref 60) is used: 


where "f| 

r 
w 

n 

5 

Since n is so large, it is sufficiently accurate to ignore the precise 
shape of the profile and to assume: 

6 , 

dll - 0.95 6, O'*) 

In the outer region (t|>6j) the equation of Harris ('Ref 61, Eq 21 gi\ s): 
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"Q j = 0.065s‘ (Spalding, Ref 62 Table 1) 


(112) 
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Typical profiles are shown in Figure 18. 

Temperature Profile 

The same assumptions as before (Equations 103 an< * 104) are made 
about the relation between the velocity and temperature profiles. 

Calculation of Jet Properties 

The equations for the properties of wall jets are similar to those 
for penetration jets (Equations 105 to 108). The major difference is in 
the expression for elements of area. 

A. = 

J 

m. = 

J 


h . = 

J 

Jet-Mixing Models 

Three of the mixing models use the jet properties calculated 
above to predict the rate at which the penetration and wall jets mix into 
the gas stream. The quantities calculated are the mass flow, axial- 
momentum flux, and enthalpy flux remaining in the jet as functions of 
distance downstream from its point of origin. The fourth model assumes 
instantaneous mixing. 

For each model the rate of mixing is governed by an entrainment 
constant, C, which is selected by the user, along with the mixing model 
to be used, as input to the computer program. The exact function of the 


C. 1) on the inner wall 
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entrainment constant varies from model to model, as is explained below. 


Mass-Loss Model 

This is a simple model in which a specified fraction of the 
initial mass flow in the jet is assumed to mix with the main stream for 
each unit of distance along the jet. (In this model the above more 
sophisticated calculations of jet properties other than length are ignored.) 


The mass flow remaining in the jet, m^ , is then given by: 

ro . » m (1 - C s ' /d . ) 

j h j 


(H7) 


The entrainment constant, C, may be assigned any value between zero (no 
mixing) and one (instantaneous mixing). 

The other jet-f,ow quantities are calculated directly from the 
residual mass flow, with the jet velocity, density, and temperature assumed 
un i form: 


r?». 


A. 
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p u. cos f. 
*an j ,o ! j 


( 118 ) 


77). ■ m. u. cos 
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(119) 


- m. 
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e p < T an 


T. ) + y U . 
b J >o 


( 120 ) 


Equivalent-Entrainment Model 

This model starts with the actual mass flow in the jet, m. 

J » «Ct 9 

as estimated from Equation 1 06 for penetration jets or 1 1 4 for wall jets. 
In these equations the integration of mass flow is carried across the 

A 

jet to a value of the jet width, 7| , large enough that virtually all 
of the jet is included. The excess of this mass flow over the initial 
flow through the hole, rh ^ - m^, is the mass flow of gases that have 
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been entrained into the jet from the main stream. It is then assumed that 
a fraction C of this entrained mass flow represents the jet mass flow 
that has been mixed into the main stream. The residual mass flow in the 
jet is then given by: 

*j ' *h * C (l Yact * V (l2,) 

An entrainment constant of zero again specifies no mixing; one specifies 
complete nixing when the computed flow of the actual jet reaches twice 
the initial jet flow; infinity (or a large number) implies instantaneous 
mixing. 

Jet velocity, temperature, and density are again assumed uniform 
and the area, momentum, and enthalpy flux are calculated from equations 
118. 119, and 120. 


Profile-Substitution Model 

In this model use is made of the jet temperature profile 
calculated by the methods discussed above. The boundary of the jet is 
defined to be that point on the transverse temperature profile where 
the difference between the jet temperature and the temperature of the gases 
in the main stream has fallen to a specified fraction C of its initial value. 

T. * - T , 

-iJ -±* - C (122) 

j ,o ’ T f t 

Air outside the boundary at is assumed to have mixed with the main 

stream. The value of 7^ defined by Equation 122 is then used to obtain 

m , A ,’>m . , and h. for penetration or wall jets from Equations 105, 106, 107 
J J U J 

and iOd or from Equations 113, 114, 115, and 116, respectively. 

An entrainment constant, C, of one implies instantaneous mixing; 
zero implies no mixing. 
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Aerodynamic Treatment of Porous Walls 

When transpiration cooling is specified, the entire flame** tube 
wall is assumed to be made from material of uniform porosity. The usual 
number of holes for the admission of penetration air are permitted, but 
there should be no film-cooling slots. 

The coolant flow through the wall is governed by the Darcy 
Equation (Ref 63 ): 



t 

w 


2RT ]i G 
w 'an 

^o 


(123) 


where t 


wall thickness 


w 

T * wall temperature 

w 

p. dynamic viscosity 

G = mass flow rate per unit area through the wall 
M « molecular weight of coolant 

The porosity of the wall is characterized by the permeability 
coefficient K which is defined by Equation 123 and has the units of ft 


2 


Calculation Procedure 

When transpiration cooling is specified, the only additional 
input required is the permeability coefficient. The calculation proceeds 
in the following steps: 

1. The air-flow subproqram is run, using an assumed wall 
temperature in Equation 123 to find the transpiration 
cool i ng flow rate. 

2. This air-flow rate is used in the heat- transfer subprogram 
to obtain a first estimate of wall temperature. 



3. This wall temperature distribution is used in the air-flow 
subprogram to obtain a new estimate of the transpiration 
flow rate. 

4. The heat- transfer subprogram is entered once more to obtain 
a new wal 1- temperature distribution. 

It is not possible, without major alterations in the program, to 
continue this iteration until any desired accuracy is achieved. Results 
show, however, that the above procedure gives adequate predictions. 


Flow and Heat Addition in the Flame Tube 

The equations of momentum, energy, continuity, and state will 
now be set up and solved for a section of the flame t ibe between adjacent 
calculation Stations 1 and 2. Contributions to the flow can be visualized 
from the sketch below. 
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Flame tube 
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The flame-tube 
residual jets. Since 
one half, or less, of 


Momentum Equation 

equations include terms for the entering and 
the gas velocity in the flame tube is typical 
the velocity in the annulus, the friction 


iy 
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loss has been ignored* 


The equation expressing the conservation of momentum within 
a control volume bounded by the flame-tube walls and by transverse planes 
at Station 1 and 2 is then 


ft, l ft,1 


+ m 


ft , 1 U ft 


l /9 o 


ZK|j 
entering J 
lets 


residual J 
j ets 
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•X 


p ft 4Ldx 


dX 


“ P ft,2 A f t , 2 + U f t , 2 /g o + aJ1 je ^j,2 


where 


fj 

*fr 


axial component of jet momentum 

m. u. cos 
J J » J 

local angle between jet axis and flame- tube center line 

i 

mass flow of hot gases in flwu tube. (This excludes unmixed 
air in the residual jets but includes the mass of fuel 
burned between stations 1 and 2.) 


The integral pressure term is approximated by: 

x. 


p ft {£ ■ d * 'I < A ft,2 ' A f tl l> (p ft,l + Pft.2 5 
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The final form of the flame-tube momentum equation is then:. 


p ft, l \t,l + "ft, I U ftV 9 o + (A ft ,2 * V.l* + P ft, 2 )/2 + 

Sn Jl0 = P ft ,2 A ft,2 + m ft,2 u ft,2 /9 o + ^ 1 V ij > 2 


entering residual 
jets jets 


al 1 
jets 


Energy Equation 


The equation expressing conservation of energy between adjacent 
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calculation stations in the flame tube is written as follows: 


(j’( T ft.i> + u ft,i / 2 v] *ft. 


, + h- m. + ^ h. ,m. . 

1 . 1 n j , o j , o ^-r* , j , 1 j,1 
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2 i 4. J , 2 j , 2 

al 1 jets J J • 


(125) 


where the enthalpy of hot gases is defined by: 

-TV 


h <V ■ 


c p ( T ) dT 


The specific heat c^ is correlated as a function of temperature by 
Equation 84; the enthalpy of the jets, h^ , is obtained from the jet- 
mixing calculation, Equation 108, 116, or 120; and the heat-release rate, 
q, from Equation 86. 

Heat loss from the hot gases to the flame-tube walls has been 
ignored in Equation 125- Strictly speaking, this is quite unjustified 
since the heat loss reduces the effective value of q, typically be about 
5 per cent in the primary zone. Inclusion of this effect would, however, 
considerably complicate the calculation procedure. 


Continuity Equation 

The equation expressing the conservation of mass in the flame 
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(126) 


where 


fbfk = fuel burning rate 
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The jet areas and mass-flow rates are obtained from the results of jet-mixing 
calculations, Equations 105, 106, 113, II 2 *, 117, 118, and 121. 


Equation of State 

The equation of state may be written as follows: 


P ft.l „ P ft,2 

At, I T ft,l At, 2 T f t ,2 


(127) 


where R * gas constant for air 

It has been assumed in Equation 127 that R ' s independent of temperature. 

This is true to within 0.5 per cent for all conditions encountered in 
practice. It has been further assumed that the gas constant for a mixture 
of combustion products is equal to that for air and independent of the 
fuel-air ratio. The actual gas constant for a mixture having initial conditions 
of 60 deg F and 4 atmospheres varies only 0.2 per cent as the fuel-air 
ratio varies from zero to the stoichiometric value. 


Method of Solution 

The mass flow at Station 2 is known from jet-mixing results: 


,2 ' "ft,! + ^* m j,o + £ m j,l ■^* m j,2 


(128) 


Since Equation 125 cannot be solved explicity for T^ an 

iterative technique is used to arrive at a solution. In this equation the 

velocity terms make a small contribution to the total enthalpy. For the 

first stage in the iterative process, therefore, 2 * s set c< l ua l to 

u -. .. In addition, the specific heat, c , and heat-release rate, 4. are 
tt , 1 r p 

computed for T^ t ^ ** y With these simplifications Equation 125 can 
be solved explicitly for a first estimate of T^ t ^ New values of Cp and 
4 are then computed for the new T ft, 2 and Equation 125 is solved again. 
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This process is repeated until T^ t 2 converges to an acceptable accuracy. 

Using the value of T^ t 2 so obtained, Equations 124, 126, and 127 
are solved to obtain a quadratic equation in^>^ t 

At?2 * T f t*2 ^ A ft,2 + A ft,l* /2 ft, 2 K 1 * A ft,2^ 9 o ^ A ft,2 *H A j,2^“ 0 
wher. K, - ^p ft ) (A ft 2 ♦ A ft>) ) ♦ m ft , « ft /9 0 

+ Dn Jie + I>lj.l -E"Ij,2 

The resulting value of fi ^ 2 '* substituted in the fol lowing version of 
Equation 126 to obtain u^ 2 : 

“ft, 2 " il ft,2 / /ft,2 * A ft,2 *E A j ,J 

This value of 2 ^5 inserted in Equation 125 and a new value 
of T ft 2 is found. This process is repeated until both T ft 2 and u^ ( 2 
have converged to within specified tolerances. 

Structure of the Subprogram 

The air-flow subprogram consists of: 

1. A subroutine which sets up starting conditions and directs 

$ . 

the air-flow calculation through the other subroutines in 
the appropriate sequence* 

2. Subroutines for solving the equations for flow in the annulus 
and flame tube* 

3. Subroutines for calculating the heat addition due to fuel 
burning and the temperature in the primary zone. 

4. A subroutine for calculating the rate of mixing of the 
entering jets* 

Detailed flow charts and program listings for the air-flow subroutines are 
given in Volume II. 
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HEAT-TRANSFER SUBPROGRAM - 
ANALYTICAL METHODS AND PROGRAM DEVELOPMENT 


I ntroduct i on 

The principal operations in the heat- transfer subprogram comprise 
the formation of the various heat flux components at the flame-tube wall 
in terms of the wall temperature, the solution of the heat-balance equa- 
tion to determine the wall temperature, and the evaluation of the wall heat 
fluxes. This section of the report presents the assumptions and relations 
employed in the subprogram. Program options which relate to the heat- 
transfer subprogram are also discussed. 

Heat Balance in Flame Tube 

The heat-flux balance is performed on an elemental length of the 
flame- tube wall; conditions in the circumferential direction are assumed 
to be uniform. The heat-flux components considered in the balance are: 

1. Radiation from the flame to the flame-tube wall, R^ . 

2. Convection in the flame-tube gases at the wall, Cj. 

3. Radiation from the flame-tube wall to the outer casing, R^. 

£. 

4. Convection in the annulus oir at the outside surface of the 
flame-tube wall, C^. 

5. Longitudinal conduction in the flame-tube wall, AK. 

6. Radiation interchange between the flame-tube walls, R^. 

7. Heat transfer between porous wall and transpiring gases, 

Q. , when transpiration cooling is specified. 

Equating the heat gained by the element of wall to the heat 
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lost, the following equation results: 

Rj ♦ Cj +4K * R 2 ♦ C 2 + Rj + G tr (130) 

Each of the components of heat flux depends upon the flame-tube 
wall temperature in the manner shown in the following sections. Substitution 
of the various heat-flux expressions into equation 130 yields a polynomial 
in the wait temperature which is solved numerically. The methods used 
for solving this equation are described in a later section after the 
various heat-flux expressions have been derived. 

Calculation Potions and Assumptions 

Several of the assumptions associated with the heat-transfer 
subroutine apply only under certain program options. The remainder of the 
assumptions are more general and apply in all cases; these latter assump- 
tions will be considered first. 

General Assumptions 

The following is a list of the general assumptions which apply 
under alt program options: 

1. For the purpose of heat-transfer calculations, the combustor 
can be considered as two-dimensional, that is, conditions 
are assumed to be constant In the circumferential direction 
for the annular geometry and in the lateral direction, with 
negligible edge effects, for the rectangular geometry. 

2. Steady-state conditions prevail, 

3. The quantities supplied by the air-flow subprogram, for 
expmp1e f temperature end velocity distributions in the flame 
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tube and annuti, are reasonable approximations to real 
conditions. This assumption is potentially one of the 
largest sources of error in the heat-transfer calculations. 

4. The flame and walls are gray radiators, that is, spectral 
effects can be ignored. In the case of the wall, spectral 
effects are unlikely to be very significant; flame emission 
is discussed again in a later assumption. 

5. Wall emissivities and absorptivi ties are constant. In fact, 
they are functions of temperature, but the range of varia- 
tion in emissivity is small. 

• 6. The effect of the wall-to-bulk temperature ratio can be 

neglected; fluid properties used in the heat-transfer 
relations are evaluated at the bulk fluid temperature. With 
this assumption the rate of heat transfer to the annulus 
air may be overestimated by as much as 10 per cent for Reynolds 
Numbers greater that'. 10,000. 

Calculation Options and Additional Assumptions 

A list of the program options which are concerned with the heat- 
transfer subprogram was presented in an earlier section of the repo. t 
Here the various options are reconsidered in the light of the additional 
assumptions associated with them. 

1. The options to specify film or transpiration cooling of the 
flame-tube walls introduce all the assumptions associated 
with the methods employed for treating these cases, namely 
the method of Spalding (Ref 9) for film cooling and the 
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method of Wheeler (Kef 64) for transpiration cooling. 

2. The use of empirical correlations to estimate the emissivity 
of the flame introduces large uncertainties. One of the 
major assumptions here is that the effect of the fuel type 
on the flame luminosity factor may be characterized by the 
hydrogen/carbon ratio of the fuel, 

3« The one-dimensional radiation option assumes that the tempera- 
ture gradients in the flame-tube are small, a situation which 
is unlikely to be realized in practice. In the case of the 
two-dimensional radiation option, it is assumed that there 
are no sudden temperature changes at the ends of the combustbr. 

4^ The exclusion of longitudinal conduction and radiation inter- 
change between the flame-tube walls from the heat-balance 
equation assumes that temperature gradients along the walls 
are small and that the two walls are at approximately the 
same temperature. These assumptions are, in general, quite 
valid; however, the option does exist to include either or 
both of these effects if they are thought to be significant 
in a particular case. 

5. Unless the outer-casing temperature distribution is specified 
in the program input, it is assumed that the casing tempera- 
ture is close to the compressor-discharge temperature. This 
is not a serious assumption since the radiation from the flame 
tube to the casing is generally small compared wi th the 
convection to the annulus air. 
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Radiation From the Flame to the Flame-Tube Wall - One-D imens i ona 1 Model 

The one-di mens i ona 1- radiation model describes the case where each 
element of wall receives radiation only from elements of flame at the 
same axial position. This model is accurate when axial temperature 
gradients and end effects are small and the length-to-diameter ratio of 
the flam'' tube is large. 

The usual expression for radiant heat transfer between a gray 
flame at temperature T^ and a black container at temperature is: 

- 0 < E ft T n ' a ft T w> 


where 


Stef an-Bol tzmann constant 


= emissivity of flame at temperature 

0L_ t = absorptivity of flame to radiation of temperature 

Since the wall is not black, but has an absorpt i v i ty a which 

w 

win normally be between 0.8 and 1.0, a correction should be made for 

absorption, reflection, re-absorption and so on. McAdams (Ref 26 » p. 90 

suggests a factor (l+QL w )/2 to allow for this effect. 

A further simplification is introduced by the use of an empirical 

formula suggested by Lefebvre and Herbert (Ref 4) ; 

r* t -i l -5 


Then 


R l -T*' +< *w> fift T f ,'' 5 < T ( 


S.5. T 2.5 } 

« l.i * 


(131) 


The emissivity of the flame, is either calculated from one of several 
alternative empirical equations presented below or supplied as special input 
as a function of pressure and temperature. Another way of obtaining 
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flame-emissivi ty data would be to use tables of empirical information 
compiled by Hottel (Ref 26) and others. Although these data are the most 
accurate that are available within their range of applicability, the 
alternative approach of supplying data in functional form has been adopted 
here. This Is equally valid, because: 

1. The functions used fit the empirical data to within one or 
two per cent. 

2. Most of the empirical data apply only at pressures near 1 
atmosphere; the excessive extrapolation required to obtain 
emissivities at normal gas-turbine operating conditions 
renders the data less useful. 

3. The only comprehensive data are for nonluminous flames; 
it is anticipated that luminous-flame correlations will 
prove more useful. 

Nonluminous Flame - Distillate Fuel 

This correlation was produced by Reeves (Ref 65). 

e ft - 1 - exp T- 18.5 p 2 0 b f x )°' 5 T f "' •*] (132) 

where f * local fuel -air ratio 

1^ ■ mean beam path length 
* 3.6 x Volume/Surface area 


HO 


Nonluminous Flame - Residual Fuel 


This correlation (Ref 65) is similar 


ft 


1 - exp -154 P 2 °' 75 O b f x )°’ 5 T 


Lefebvre Correlation 




Lefebvre, in Reference 4, used Reeves's distillate-fuel 
correlation with a "luminosity factor", A : 


£ ft = 1 " exp 


- 18 . 


C A ( y t \ 0 • 5 T ” 1 • 5 
5 Po A (K f J T r«- 


2 #1 v ’b'x' 'ft 
A correlation for A was based on data from two sources. 

A = 7.53 (C/H - 5.5) 0,84 

where C/H = carbon to hydrogen ratio of the fuel, by mass 
Thi;, expression is invalid below carbon-hydrogen ratios of about 5-7 
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0 34) 


(135) 


NREC 1964 Correlation 

By 19^4, further data on the effect of C/H on A were available. 
An examination of these sources showed that a better correlation for A 
i n Equat ion 134 is: 



(136) 


NREC 1966 Correlation 

Schirmer and Qu i gg , in a recent paper (Ref II), examined the 
effect of pressure on flame radiation in a 2-inch diameter tubular combustor. 
Their data have been converted from radiation flux to emissivity at an 
assumed flame temperature of 3500 deg R and are plotted in Figure 19, 
together with the predictions of Equations 134 and 136; agreement is 
poor at C/H ratios above about 7. A new correlation, which fits the data 
reasonably well over the whole range, has therefore been produced. This is: 


Ill 


ft 


1 - exp 


[■ 


4.7 pl‘ 3 AO b f x )°- 5 


. -1.5*] 

'ft J 
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where A " (ffii i 1 ) 

This correlation is compared with the experimental 
in Figure 20. 
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Radiation From the Flame to the Flame-Tube Wail - Two-Dimensional Model 
This section describes the option for calculating the heat 
received by a wall element from elements of flame at all axial locations. 
The radiative heat transfer to an element of the wall from an elemental 
volume of flame is equal to the product: 

Volume of flame element x Emission per unit volume 

x Transmittance x View factor (139) 

The quantities in this product will now be considered individually. 

Volume of Flame Element 

The elements of flame considered have elemental length in the 
axial direction of£X^ t and inner and outer (flame-tube-wall) radii of r^ 

and Tj respectively. Thus, the volume of each flame element is: 

Ax ft - n (r , 2 - f 2 2 ) d4o) 


Emission Per Unit Volume 

Following Hottel (Ref 26), the radiation emitted by unit volume 
of flame at temperature T^ is: 


4<rT ft P G 


{Wt) 


(141) 


Pc'b * 0 


where p^ = partial pressure of radiating gas (l.e., water vapor end carbon 
dioxide) , lbf per sq ft. 
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/< 


de \ 

The quantity/ ft J is chiefly a function of temperature. 

\ dp G' b / P G>L = 0 

For a typical value of p^/ p^ t = 0.05 and for a non-luminous flame it can be 


approximated by (Ref 2i) : 


dp A / 


dp G'b / P G'b = 0 


3.6 

ft 


(142) 


The form of the expressions for flame emissivity (Eqs 132 through 138) indicates 

that for luminous flames, the following equation is approximately true: 

r— “i r . "1 


d £ 


ft 


‘■Vb r gV 0 


= A x 


luminous 


/ de 


ft 


^"Vb / p gV 


non-lumi nous 


Thus a more general form of Equation 142 is: 

'd £ 


ft 


dp G’b 


p G’b = 0 


3.6 A 

T ft 


(143) 


View Factor 

The sketch below illustrates the nomenclature of the axially 
symmetric configuration considered. 


Elemental 
length of 
f 1 ame 


V 


r 3 


Elemental length 
of wall 




Engine center 1 i ne 


H3 


Employing an equation derived by Hamilton and Morgan (Ref 66 ) 
and the principle of superposition, the following equation is obtained 
for the view factor of the flame element radiating to an element of length 


AX in the outer wall: 
w 
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Transmi ttance 


The transmittance accounts for the rad i at ive energy which i s not 

* 

absorbed by the intervening gas; it is a function of the distance the 

» * 

radiation travels through the gas, the gas composition and the gas * 
temperature. Here a curve has been fitted to typical data prepared/ by % 
Hottel (Ref 26) and shown in Figure 21. The expression developed Ijs : 

p , 14,82 # * 

I p. s + l 4. 82 % 

G * 

9 

t ‘ 

where s = distance from element at considered point on the wall to* 

0 

element of the flame, ft * . 

Since s varies over different parts of the flame, an averaging procedure 

\ 

is used over the radial dimension of the flame to determine a single,, 
mean value of the transmittance for each axial location of the flame 
element. The resulting expression for the transmittance to an outer wall 
element is: 
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r mean 20 ^ 


14.82 
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i =1 
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' 82 + P g/ * 2 + [ r 3 - r 2 -(| 0 - 55) (r, - r 2 )] 2 


Partial Pressures 

The partial pressures of carbon dioxide and water vapor that 


are required for Equations 141 and 146 are obtained from stoichiometry 
Then: 

P G = P ft (x C0 2 + x H 2 0^ (147) 

where x = mole fraction of constituent in product mixture. Local air- 
fuel ratios are calculated in the air-flow subprogram in such a way 
that the mixture can never be richer than stoichimetric. 


Total Radiation 

Introducing a wall temperature, T^, and an effective wall 

absorptivity, ( 1 + o< )/2, the radiation absorbed by the wall element 

w 

from the flame element is: 


4 aT 
mean 


(1 


+ ck ) 
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P, (- 44 ) AX f( TT(r 


dp G L 






Hence, the total radiation per unit time at a point on the outer flame- 
tube wall receiving radiation from all positions of the flame, is: 

4t nr<’ + £ 

2 a 1 1 f 1 ame 
pos i t i ons 

The radiative heat flux per unit time and area of flame tube is, therefore, 
given by: 

AX R_ all flame 
W ^ posi tions 

A similar expression applies for radiation to the inner flame-tube wall. 

It should be noted that it is not necessary to specify the length 
of the wall element, AX , since in Equation 143 AX in the denominator 

W 

cancels with AX w appearing in the expression for the view factor, Fp^. 
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Radiation from the Flame-Tube Wall to the Outer Casing 

The radiant heat flux from the flame-tube wall to the outer 

casing is given by the equation: 

R„ - F 0-(T 4 - T 4 ) 

2 wc w c 


(149) 


where T = temperature of casing. This may be supplied as input as 

c 

a function of axial length; if it is not, it is assumed 
equal to the compressor discharge temperature. 

F r - overall interchange factor which accounts for the re-emission 
and absorption between surfaces. 

For the case where only two surfaces are involved and where both surfaces 
have high emi ss i vi t i es , the overall interchange factor is approximately 
given by (Ref 26, p. 76): 


f~ -[ e - - ' 

wc \>* w 
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w 


- — i V s' 

«c ' F 


wc 


(150) 


where 


F = black-surface overall interchange factor 
wc 3 

£ = emissivity of casing (assumed equal to absorptivity) 

c 

A = surface area of flame-tube wall 
w 

A^_ = surface area of casing 

For radiation across an annular space, the black-surface interchange factor 

can be assumed to be equal to unity. With this simplification and some 

rearrangement, Equation 1 50 becomes: 

_ £ w £ c 


wc 


£ c +£ wf - £ c) 


w 


(151) 


A 


Thus, from Equations 149 and 1 5 1 the radiant heat flux from the wall to 
the outer casing is: 




R, 


= 0 


e w e c 

+ C , ( 1 - e ) A /A 
w c w c 


(T 4 - T 4 ) 
w c 


(152) 
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Radiation Interchange latwcqn Flame-Tube Wells 

Under some conditions, radiation interchange between the combustor 
walls can become important. The net rate of heat transfer to a section of 
the fiame-tube wall from all other "visible" parts of the flame-tube wall is 
derived using the nomenclature shown in the sketch: 



Radiation emitted by element which is received by ele- 

ment (neglecting any reflections) is: 

• a * . 

* - 2 L L 

€t 2 TTr 6x F p(T . - T ) 

• w . b b ba wb wa 


( 153 ) 


when 


F, = View factor for radiation from b to a 
oa 

• '' r. - r .» Z 

r = Transmittance between b and a 

14.82 

PrWK - x „) 2 + K - r J 2 + ,1 *' 82 


(I5 1 *) 


Radiation interchange per unit time and area of inner wall is given by: 
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r * 
a 


ix b < r b- O 

2 [ (r b ' r a )2 + (x a ‘ *b )2 ] 3/2 


( 155 ) 


Thus the net rate of heat transfer per unit area to the point on the 
inner wall under consideration, from all parts of the outer wall is: 


A similar 
wall to a 


r = <n ‘ 

3 w 


y ^ 


,> 2 r(T wb 4 - t . 6 , 


wa 


b-wal! r a [(r b - r ^ 2 + (x, - x fa ) 2 J 


all 
pos i t i ons 


213/2 


expression applies for radiation from all parts of the inner 
point on :he outer wall. 


Convection at the I nner Surface of the Flame-Tube Wall 

The correlation employed for the internal convective heat transfer 
at the flame-tube wall is one suggested by Humble, Lowdermi Ik and Desmon 
(Ref 67) for subsonic flow through smooth tubes at high surface and fluid 
temperatures : 


c 1 = 0.023 


D 


Re°- 8 Pr °- 4 (T 


ft 


V' 


ft 


(156) 


where k = thermal conductivity of gas 

d| = hydraulic diameter of the flame tube, equal to twice the 
distance between the flame-tube walls for both annulus and 
rectangular geometries. 

The effect of annulus diameter ratio is small and is ignored in this 
treatment. Reynolds and Prandtl numbers and thermal conductivity are 
based on properties of a gas mixture taking account of the amount of fuel 
present. The methods used to evaluate these properties are described in 


a later section. 
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Entry-Length Effect 

When starting conditions on the wall are uniform, there is a 

developing boundary layer and the expression for the Nusselt number is 
usually multiplied by a factor for the entry-length effect. Typically 
(Ref 26} , this is: 


1 + 


-Jt 


-0.7 


Because of the intensely turbulent conditions in the combustor, the 
relationship for fully developed flow (Equation 1 56 ) gives a better 
representation and has been used in the program. 


Internal Convection with Film Cooling 

The contribution of internal convection to the wall heat bal- 
ance is computed from Equation 156 only when film cooling and transpiration 
cooling are absent. For the case in which film cooling is employed, 
Equation 1 56 is modified by replacing the hot-gas temperature, T^ t> by an 
adiabatic wall temperature, T . , thus: 

dCl 


C. « 0.023 Jr~ Re 0 * 8 Pr 0 * 4 (T . - T ) 
1 ad w' 


057 ) 


The adiabatic wall temperature is a quantity measured experimentally 
and has been correlated by Spalding (Ref 9 ) to geometric and flow variables 
by means of a film-cooling effectiveness, X , and a downstream-distance 
parameter, X: 


, - T ft ~ T ad 
= T ft - T cf 


u,„ 0.8 -0.2 

- 0.91 (TT^-TT") (Re cf ) + 
U ,. y cf 


cf 


cf 
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ft 
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where T ^ * inlet temperature of cooling film 

y c f * slot width 

x = downstream distance of wall element from cooling slot. 

Figure 22 shows an envelope enclosing 90 per cent of the data 
points employed by Spalding; these data may be represented by the three 
straight lines: 

1*1 

X 0.22 

or z - (y^) (160) 

X 

or T ■ ~ whichever is the smallest 

A 

where X * 7.0 

o 

These relations correlate well the nine sets of experimental data on which 
they were based, but their validity for practical combustion systems, in 
which the slot may be obstructed by structural elements and imperfections 
left after fabrication, is uncertain. This uncertainty is supported by 
the data published by Sturgess (Ref 10} which are also plotted on Figure 
22. These data for, so-called, “dirty" slots may be represented by the 
same equation (Eq 160) that was used to represent the data for “clean" slots; 
however, the value of X Q is reduced from 7.0 to 3.5. 

Equation 160 is employed in the computer program and X^ is 
supplied as program input by the user; its value can be determined in one 
of two ways: 

1. If experimental data are avai lable, these are plotted in 

the form T — X, and a value of X^ that fits the data best 

o 

is chosen and supplied as input to the program. 

2. if no experimental data are available, the recommended values 
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for X are: 3*5 for "dirty 1 ' slots; 7.0 for "clean" slots, 

o 

The procedure employed by the program to determine the 
convective heat flux is as follow: For a particular wall element 

considered, a value for X is computed. This value, together with the 
value of X Q supplied as input, are used to calculate the cooling effective- 
ness X from Equation 160. Equation i$8 is then employed to give the 
adiabatic wall temperature, T which is substituted into Equation 157 
to give the convective heat flux. 


internal Convection wi th Transpiration CooKno 

The transpiration-cooling flow rate is obtained from the air-flow 
subprogram, it is assumed that the coolant leaves the wall at the wall 
temperature. x 


The heat-transfer rate is calculated by a method suggested by 
Wheeler (Ref 64). With this method, the convective heat-transfer coefficient 
for the transpiration-cooled wall is given by: 

Nu - 0.002 Re ft 

The Nussait and Reynolds numbers are based on the bulk properties of the 
gas in the flame tube. 

Hence, the convective heat flux from the hot gas to the wall is: 

C, - 0.002 R. ft (T ft - T w ) (161) 


It should be emphasized, however, that no independent assessment of the 
validity of this method for the present application has been made. It 
is expected that* If a serious study of transpiration cooling is to be 
undertaken with this program, more sophisticated and powerful methods 
such as that of Spa Id I ng (Ref 68) , or actual teat data from pract I ca I 


combustors, will be used. 
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Convection at the Outer Surface of the Flame-Tube Wall 

External convection from the flame-tube walls to the annulus 
air is governed by an equation similar to Equation 1 56 : 

C 2 = 0.023 4 - Re°; 8 Pr 0 - 4 <T„ - Tj (162) 

an 

where 0 is equal to twice the distance between the flame-tube and 

an 

outer casing for both annular rectangular geometries. 

In the case of transpiration cooling, there is an additional component 
Q representing heat transferred to the coolant as it passes through the 
wall: 


Q 


tr 


G C. 


< T w- 


T an> 


( 163 ) 


Gas Properties 

The gas properties used in the above relations are obtained 
from empirical worrelations derived from a survey of all available 
literature. Values for air are used in the equations for heat transfer on 
the annulus side of the wall. In the case of gases in the flame tube, the 
molal gas composition is found from stoichiometry, and the properties of the 
constituents are combined according to formulae recommended by Perry (Ref 69) 
to obtain the properties of the mixture. 

The effect of pressure has been ignored. 

The equations used for the specific heat, the dynamic viscosity 


and the thermal conductivity can all be expressed in the form: 


Z = a + bT + C T 2 + dT 3 


064 ) 


where 


Z represents the gas property 
The constants take different values for each property and for 
each gas considered as shown in Table 1. Table I also shows the temperature 
range covered by the data employed to determine the values of the constants 
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for each gas; the percentage accuracy of the equations is also indicated. 


Properties of Gas Mixtures 

The specific heat of a mixture is obtained from: 

C * C . 

p * P» 

m. * mass fraction of ith component 
C p . = specific heat of ith component 

The viscosity of a mixture is obtained from: 


where 


where 


H * 
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X. 

I 


2 x i < M i) 

mole fraction of ith component 


M. = molecular weight of ith component 


( 165 ) 


( 166 ) 


The thermal conductivity of a mixture is obtained from: 

1/3 


k - 


S'x. k. (M.) 

i i i 
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Longitudinal Conduction Along Flame-Tube Wall 

The longitudinal conduction effect can be estimated on the 

basis of average temperature gradients. For the inner flame-tube wall 

the rate of heat flow to the section of wall (n) under consideration 

from an adjacent one (n + 1) is given by: 

k t C A (T . - T ) 
w w A w,n+l w,n' 

c * 

s x . I - x 

n+l n 

where k * thermal conduct ivi ty of wall material 
w 

t = wall thickness 
w 
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Per unit area, the heat gain is: 

k t (T . 
K _ w w w.n+1 
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Combining this with a corresponding quantity for heat flow from the 


n-1 element, the net heat received becomes: 


AK - 


2k t r T 
w w w 

i . i - x , 
n+l n-1 


,n+l 


- T 


w.n 


n+1 


T - T 7] 
- _w,_n w,n-l 

x - x - 
n n- 1 J 


(169) 


Solution of the Heat-Balance Equation 

When all of the heat-flux terms are considered, the heat-balance 
equation for an element of the flame- tube wall is: 

Rj + Cj +AK = 1*2 + ^2 + ^3 + ^tr 

Both the longitudinal conduction term,AK,and the radiation interchange 
term, R^,are functions of the wall temperature distribution (rather than 
the local temperature of the wall element) and, therefore, an iterative 
method of solution of liquation 130 is required when these terms are 
considered. Often, however, the user will not wish to consider A K and 
Rj>in which case a non-iterative solution of Equation 130 is possible. 


Non-Iterative Solution 


In this case Equat on 130 becomes: 


R, + C, = R 2 + c 2 + d tr 


(170) 


Q ^ is zero when there is no transpiration cooling. 

Substituting for R ] from Equation 131 or 146; for C ] from Equation 156, 

157, or 161; for R 2 from Equation 149; for C 2 from Equation 162; and 
for Q from Equation 163; Equation 170 becomes: 

°1 T w + °2 + D 3 T w * D 4 (,7,) 

where Dj, D 2> D^, and are known, their values depending upon 
the program option. Newton's approximation is used to solve Equation 171 • 

For use in Newton's approximation, a first estimate of T w is 
obtained from one of the following equations: 
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The chosen T is that which leads to tfee smallest value of £ when Equation 
w 

171 is written in the form: 


D. T 4 + D 0 T 2,5 + D, T - D. 

1 w 2 w 3 w 4 
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With this initial estimate of T Newton's approximation is 

W » 

used to find a value of that satifies Equation 171* 


I terative Solution 

The iterative heat-balance solution is only called for after a 
noniterative calculation; an approximate wal 1- temperature distribution will 
therefore always be available. The iterative heat balance is repeated 
until the maximum temperature change between any two iteration cycles is 
less than the specified tolerance. 


End Effects 

The boundary conditions assumed at the ends of the combustor can 
have an important effect on the overall wall temperature distribution. 

Since the highest temperature gradients frequently occur away from the ends, 
it has been assumed that the wall temperatures at the ends are as given 
by the noniterative calculation. 
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Structure of the S ubprogram 

The heat-transfer subprogram consists of: 

1* A subroutine for solving the noniterative heat balance* 

2. A subroutine for solving the iterative heat balance* 

3* A group of library subroutines which provide gas 
properties, correlations for wall cooling and flame 
emissivity, and the solution of the fourth-power heat- 
balance equation. 

Detailed flow charts and program listings for the heat-transfer subprogram 


are given in Volume II. 
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L IMITATIONS OF THE COMPUTER PROGRAM 

The computer program represents a significant advance in the 
analysis of annular combustors, in that it draws together reasonable represent* 
ations of the pertinent flow, burning, and heat-transfer phenomena Into an 
integrated computing tool. Because of limitations In the models and in the 
experimental data available, and the need to keep the program to a manageable 
size, however, it has been necessary to make a number of assumptions and 
compromises that limit the accuracy of the results* Most of these have been 
discussed in connection with the analytical methods* To emphasize the 
importance of understanding these limitations in making intelligent use of 
the program, they are summarized here. 

General 

The analysis is one-dimensional : circumferential variation of flow 

quantities such as pressure and temperature is ignored throughout the program*, 
radial variation of these quantities is considered in the streamtube method 
of diffuser analysis, but not elsewhere* It Is well known that the flow in 
an actual combustor departs significantly from one-dimensionality, because 
of the presence of swirlers, fuel atomizers, holes, swirl at the compressor 
exit, and so forth. At present, however, neither the theoretical methods nor 
the experimental data are available upon which to base a more detailed model. 
This is, then, a limitation of the program as a whole. Its effect is difficult 
to assess. It should be kept In mind when using the program in connection 
with experimental data, however, that unusually large departure from one- 
dimensionality In the experiment may cause considerable disagreement between 
experiment and analysts* 
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Diffuser Subprogram 

The analytical methods used in the diffuser program contain several 
important limitations: 

1. Flow in the diffuser is assumed to be unaffected by downstream 
conditions. If the flow is well matched at the snout, this 
assumption is expected to be valid. Otherwise, serious errors 
will undoubtedly result. Snout design should, therefore, receive 
careful attention. 

2. No provision is made for including splitter vanes in the diffuser 
analysis. 

3. The program does not attempt to treat flow after separation. 
Air-Flow Subprogram 

The more important limitations of the air-flow subprogram are as 

follows: 

1. The treatment of the primary zone as a stirred reactor limits 
the information obtained. For example, if the detailed burning 
rate and flame temperature were known, the flame-tube wall 
temperature could be calculated as a function of position. Present 
understanding of the flow, mixing, and combustion processes in 

the primary zone is Inadequate to support a more detailed analysis. 

2. The program does not calculate the fuel-burning rate directly, 

but receives it as input and only serves to keep it below stoichio- 
metric at any location. The results are limited by the judgment 
of the user and by his access to experimental data. For a more 
accurate estimate of flow conditions, particularly at low pressures 
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where reaction rates decrease, it should be possible to devise 
a fuel-burning rate correlation. 

3. The jet-mixing models suffer from a lack of experimental data 
upon which to base the choice of models and entrainment constants* 
in addition, the program does not relate the mixing rate of wall 
iets to the film-cooling analysis carried out in the heat-transfer 
section. 

4. The swi rler-design method contained In the program is crude. 

When more experimental data become available, a more refined 
method should be incorporated. 

5* IT the use of transpiration cooling is to be explored in detail, 
an assessment of calculation methods should be made. 

Heat-Transfer Subprogram 

There are several significant limitations to the heat- transfer 

program: 

1. The film-cooling model is idealistic; there is a serious lack 

of experimental data upon which to base choice of the characteristic 

distance X used In the correlation; in eddition, as mentioned 
o 

above, the film-cooling and wal 1 -jet-mixing models are not related. 

2. Data is lacking on the emissfvlty of luminous flames. 

3. There is a lack of data on transpiration cooling, with which to 
judge the reliability of the theory used. 

4. No heat balance Is performed to calculate the temperatures of the 
inner and outer casing walls* 


129 


Conclusion 

Before making extensive use of this program, It would be useful 
to make a more exhaustive survey of the Influence of program options than 
has been possible within the limitations of time and budget In the present 
project* For example, Jet-mixing rates, and consequently the flame- temperature 
distribution, are quite sensitive to the choice of entrainment function and 
constant* By experiment, a 'Vnost reasonable" choice may be found* 
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TEST CASES 

In this section, the results of a number of test cases 
illustrating aspects of the program are presented. 

Diffuser Test Cases 

I ntroduct ion 

The diffuser test cases are mainly concerned with checking the 
validity of the analytical "st reamtube" method. The following cases were 
speci f ied: 

1. Comparison with experimental data on two-dimensional 
diffusers, compiled by Cochran and Kline in Reference 70. 

2. Comparison with data on annular diffusers, compiled by 
Abramovich and Vasil 'yev in Reference 71 • 

In addition, the results of a number of test cases used during program 
development are presented sn this section. 

Choice of Separation Shape Factor 

Before the test cases were run, it was necessary to select a 

value for H g , the value of the shape factor at which separation occurs. 

In conical diffusers, H is typically 3-0, but for airfoils at moderate 

sep 

incidence, a value in the range of 1.8-2. 2 is more usual. 

The value of H * 1.9 is used for comparison of the 
sep 

streamtube calculation method and the empi ricai -data method. The 
latter is based in the case of two-dimensional diffusers on the 
experimental results of Reneau, Johnston, and Kline (Ref 7) • A com- 
parison between streamtube and empi ricai -data methods is shown in 
Figure 9 for diffusers with area ratios in the range 1.2 to 5*0 and 
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nondimensional lengths of 1.0 to 20.0. The only assumption made in 
obtaining these results was that the shape factor at entry is 1.4. 

For annular diffusers the empirical -data method is based on 
data of Sovran and Klomp (Ref 6). The experimental results are based 
on diffusers with inlet hub-to-shroud ratios of 0.55 end 0.70 and various 
values of 0j and 0 Q (see sketch below) and inlet blockage. 



In Figure 23 the experimental results of Sovran and Klomp are 
compared with the predicted performance using the stream tube method. The 
following assumptions were made! 

1. Shape factor at Inlet ■ 1.4 

2. Separation shape factor ■ 1 .9 
3* Blockage at inlet * 3 per cent 
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4. Hub* to- shroud ratio ■ 0.55 

5. 9. - 0 degrees 

The Inlet blockage assumed in the calculations was somewhat higher than 
the 2 per cent assumed in Reference 6 for the following reasons: 

1. The flat-plate Reynolds number at entry appears from the 
geometry of the system to be of the order of 10^ giving a 
blockage of 4 per cent. 

2. The maximum diffuser effectiveness recorded was 85 per cent 
which indicates a high inlet blockage. 

The difference between Figures 9 and 23 Is not due only to the different 
values of inlet blockage, but is also due to differences in the variation 
of area ratio with length. For a two-dimensional diffuser: 


A ■ A. + (A 2 - A.) . x 

I 

<1 ( A/ A, - ') . *2 . , 

d (*/i) A l 

For a simple annular diffuser: 



(173) 



« 
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Equations 173 and 174 reduce to the same form for r^ * r^, l.e. 
or for large values of r . Thus the performances of annular and two- 
dimensional diffusers show the most difference at large area ratios. 


Two-Dimensional Diffuser Geometries 
The test cases for two-dimensional diffusers are taken from 
Cochran and Kline (Ref 70) and the earlier work on similar apparatus by 
Moore and Kline (Ref 72). The geometry of the diffuser section is shown 
below: 



This geometry is not entirely suitable to test the streamtube method for 
the following reasons: 

1. The boundary layers at inlet are formed on highly curved 
surfaces, and not, as assumed, oh a flat plate, 

2. The static*pres$ure profile at inlet is distorted due to 





the curvature of the flow. In the streamtube method the 
static pressure is assumed to be uniform across each cross- 
sect i on. 

These difficulties notwithstanding, it is felt that the results provide 
an adequate basis for assessment of the streamtube method. It is to be 
noted that these results are also included in the compilation and cor- 
relation of experimental data presented in Reference 7* 

For the purposes of obtaining a concise comparison of the 
diffuser effectiveness obtained from the streamtube method with experi- 
mental results over an adequate range, a family of diffusers was selected 
fr om the data of Reference 70. These diffusers all have a nondimens ional 
length of approximately 8, and cover a range of area ratios from 2 to 
6.5. The geometric data are summarized below. 


1 nc 1 uded 
Ang 1 e 
degrees 

Nondimensional 

length 

L/Wj 

Area 
Rat i o 
a 2 /Aj 

7.0 

8.04 

1 -99 

14.0 

8.08 

2.96 

16.8 

8.10 

3.40 

21 .0 

8.12 

394 

28.0 

8.09 

4.94 

311 

8.06 

5.44 

35-0 

8.03 

5.94 

38.2 

8.00 

6.46 


For these configurations, measured values o* the shape factor at inlet were 
in the range 1.43 to 1.79 and measured values of the blockage were in the 
range 0.004 to 0.006. To circumvent the difficulty associated with obtaining 
the experimental velocity profile at the inlet to the diffuser, the streamtube 
method was used assuming a uniform inlet velocity profile; to account for the 
fact that the actual velocity at the wall was some 20 per cent greater than that 
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In midstream, the test cases were run with a blockage of 0.0075, and a cor- 
responding shape factor of 1.47 (since the deceleration rate along the wall 
will be greater for the nonuniform velocity profile, a larger blockage should 
be used wi th a uniform velocity profile to simulate this effect). 

For obtaining a comparison of the results of the streamtube method 
for the prediction of separation with experimental results, the data of 
Reference 72 were used. Again, the streamtube method was used with values 
of the inlet blockage and shape factor of 0.0075 and 1.47, respectively. 

Comparison of Results for Diffuser 
Effectiveness 

The effectiveness as computed by the streamtube method is shown 
with the experimental results in Figure 24. The agreement is reasonable, 
considering both the approximations made with regard to inlet conditions and 
the fact that the experimental accuracy was reported as +5 per cent. 

Results - Line of First Stall 

In Figure 25, lines are drawn dividing geometries in which 
separation occurs from those in which the flow does not separate. A 
comparison is made between the results of Moore and Kline (Ref 72), in 
which a flow visualization technique was used, and the results obtained 
using the streamtube method. Agreement is quite good. 

Annular-Diffuser Geometries 

The test cases for annular diffusers were taken from Abramovich 
and Vasil'yev (Ref 71 ) . The geometry considered is shown below: 
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LLLLLLL III 




Sol id cy li ndr i cal 
1 


i nsert 


77777 


The two test cases analyzed were: 



77777777777 


Area Ratio 

Range of 8 

V A . 

degrees 

2.0 

5-40 

38 

5-40 


The following assumptions were made: 

1. Shape factor at inlet * 1.4 

2. Shape factor at separation « 1.9 

Annular-Diffuser Results 

In Figure 26 a comparison is made of the results of the stream- 
tube method with the experimental results given by Abramovich and Vasil'yev 
(Ref 71)* As the blockage is not quoted, theoretical curves are given for 

1 per cent, 2 per cent, and 3 per cent blockage. A blockage of 1.5 per cent 

gives good agreement between the theoretical and experimental results. It is 
to be noted that the inlet blockage has a significant effect on the computed 
effectiveness, as is expected from the consideration that the effectiveness 
is determined by the displacement thickness at the exit (until separation 
occurs) and that this thickness is roughly proportional to the inlet displace- 


ment thickness. 
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Unfortunately, no data on the occurrence of separation in annular 
diffusers is available for compar i son wf th analytical results. 

Shear Flow in a Diffuser 

The streamtube method provides a means of predicting the change in 
velocity profile through the diffuser. The main assumptions made in the 
method are that the profile at entry is in equilibrium, and that mixing 
effects can be neglected. 

To test the method, the experimental resjlts of Horlock and Lewis 
(Ref 73) are analyzed. In these experiments strai ght-wal led diffusers and 
nozzles were used and velocity profiles at inlet and exit were measured. 

Two geometries were considered: 


Type 

Area 

Total Included 


Ratio 

Angle, degrees 

Nozzle 

0.7 

10 

D i f fuser 

1.47 

15 


The results of these tests are shown in Figures 27 and 28; the agreement 
between the streamtube method and the experimental results is good. 

Air-Flow Test Cases 

Introduction 

The air-flow test cases were designed to check the pressure- 
drop and air-flow-distribution calculations. For comparison, the 
experimental and theoretical work of Dittrich (Ref 74) and Graves and 
Grobman (Ref 1) was used. 
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The tests used for comparison were carried out on a straight- 
walled tubular combustor without a diffuser. Seven test cases were run, 
varying three main parameters: 

1 . Reference Mach number 

2. Overall temperature ratio 

3. Ratio of total hole area to reference area 

Input Conditions 

Since this program only deals with annular combustors and 
combustors of rectangular cross section, it was impossible to match the 
geometry used in the original tests. Instead, a hypothetical combustor 
of rectangular cross section was set up. This was 12 inches wide and 
had all cross-sectional areas equal to those in the NASA combustor. The 
dimensions of the hypothetical combustor are shown in the table below. 


Geometry employed in Test Case 

NASA Combustor (Tubular) 

Flame-tube height 2.37 in 

Flame- tube diam 6.0 in 

Casing height 3-95 in 

Casing diam 7-75 in 

Overall length 52.0 in 

Overall length 52.0 in 

No. of holes per . 

No. of holes , 

row on each wall ^ 

per row 

Area per hole 1.48 sq in 

Area per hole l .48 sq in 


input quantities specified by NASA are listed below: 
Compressor-di scharge total temperature 80 deg F 

Compressor- di scharge total pressure 15 psia 

Fraction of air entering secondary holes that 

recirculates upstream 0 


Rate of jet mixing 


instantaneous 
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The total pressure was assumed to be 15 psia at the beginning of each 
annulus . 


Three hole arrangements were used, as indicated below: 


Axial location, inches 

6 

10 

14 

18 

22 

26 

30 



34 



Run 1 

D 

D 

■ 

y/ 

B 

■ 

B 

B 

Run 1 1 

D 

D 

D 

>/ 

D 

D 

H 

a 

Run III 

D 


■ 


D 



B 


In Run I, input conditions were varied from case to case as 


fol lows : 



Overal 1 

Fuel-air Ratio 

Ai r Mass Flow Rate 
Ibm per sec 

Case 1 

0.0 

2.46 

Case 2 

0.0 

1.23 

Case 3 

0.0 

4.92 

Case 4 

0.00724 

2.46 

Case 5 

0.0152 

2.46 


The three fuel-air ratios correspond, according to Figures 21 and 23 of 
Reference 36, to temperature ratios (combustor outlet/inlet) of 1, 2, and 
3. The three mass flows correspond to reference velocities of 50, 100, and 
150 ft per sec. 

Results 

In Figure 29, the variation of pressure- loss factor with reference 
Mach number as calculated by the present method is compared with the results 
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of Reference, 74. The agreement with the experimental data is very good, 
considering the difference in geometric configurations. However, the 
sharper increase in pressure drop at higher Mach numbers, apparent in the 
experimental data, is not so marked in the curve representing the present 
predictions. 

Figure 30 shows the effect of changing the ratio of total 
o 1 c* area to reference area. Here the agreement with experimental results 
is exceptionally good, particularly at high values of the area ratio. 

The effect of heat addition on the pressure-loss factor is 
shown in Figure 31; again the agreement with the experimental results 
Ts good. 

Figure 32 is a plot of the fraction cf total air flow in the 
flame tube against the ratio of cumulative open hole area in the flame- 
tube wall to total hole area. The results obtained here differ from 
those of Graves and Grobman (Ref 1) at high values of the ratio of 
total hole area to reference area, A^^./A re ^. No experimental results 
are available to check the validity of the two methods. 

Gone 1 us i on 

In view of the differences between the geometrical configuration 
of the experimental combustor in Reference 74 and the configuration 
employed for the present computations, the comparison between the calculated 
and the measured pressure losses is particularly good. The only effect 
which the computations apparently do not accurately predict i s the effect 
of high Mach numbers on the pressure-loss factor. However, the Mach 
number at which the divergence becomes appreciable is outside the normal 
range of reference Mach numbers for gas-turbine combustion chambers. 
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Heat-Transfer Test Cases 

The predictions of the heat- transfer subprogram were compared 
with experimental data presented by Milford and Spiers in Reference 75. 

Experimental Conditions 
The apparatus used is shown in the sketch. 



Three sets of results were obtained: 

1. Convective cooling (on annulus side) only 
2-. Fi lm cool ing only 
3. Convective and film cooling 
Some of the test conditions were as follows: 

Hot-gas temperature 3510 deg R 

Hot -gas Mach number 0.25 


Convective Cooling 
(Cases 1, 2, 3, and 4) 


'Cooling-gas temperature 1440 deg R 
- Cooling-gas Mach numbers 0, 0.11, 0.14 


and 0.20 
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Fi lm Cool ing 
(Cases 5, 6 $ and 7) 

Convective and 


Cooling-gas temperature 1440 deg R 
Slot-gas velocities 260, 570, and 
800 ft per sec 

j'Cool ing-gas temperature 1530 deg R 

r t i \i£± f m imka r A 1C 


Fi 1m Cool ing 
(Cases 8, 9, 


10, and 11) 


Slot- gas velocities 0, 260, 580, and 
825 ft per sec 


Assumptions 

The information provided in Reference 75 was insufficient to 
specify the problem. A thorough search of the literature revealed no other 
published experimental results in which heat-transfer measurements had been 
made with combined film and convective cooling. Thus, other than abandon 
the heat-transfer test case, the only course left open was to estimate or 
assume values for the information missing from Reference 75« The information 
required to complete the problem specifications and the assumed values are 


as follows: 

Distance across flame tube 4.5 in 

Annulus height 0.45 in 

Slot height 0.036 in 

Static pressure Atmospheric 

Fuel-air ratio 0.067 

Fuel hydrogen-carbon ratio 0.17 


The casing temperature was calculated assuming natural con- 
vection and radiation to a room- temperature background. (The actual 
heat-transfer coefficients were obtained from McAdams, Ref 26, Table 7-2.) 
Heat transfer from the casing to the annulus air was ignored. The casing 
temperature distributions thus obtained were as follows: 



143 


Case 

Casing Temperature, deg R 

Distance 

2.0 

Downstream of 
4.0 

>1ot, in 
5.5 

1 

1997 

1872 

1786 

2 

1808 

1780 

1750 

3 

1770 

1752 

1720 

4 

1720 

1707 

1698 

5 

1800 

1762 

1723 

6 

1610 

1661 

1662 

7 

1580 

1639 

1647 

8 

1762 

1757 

1720 

9 

1566 

1617 

1609 

10 

1461 

1542 

1550 

11 

1448 

1518 

1537 


The flame-tube wall absorpt i vi ties and emissivities were taken 

as O. 85 ; the casing emissivity was taken as 0,8. The flame was assumed 

to be luminous (1966 NREC correlation. Equations 1 37 and 1 38 ) . and a 

value of X * 7.0, the value for "clean" slots, was assumed in Spalding's 
o 

film-cooling correlation (Equation loO). 

Discussion of Results 

The results are shown in Figures 33. 34, and 35. Figure 33 
shows Cases 1 to 4 - those with convective cooling only. Agreement is 
reasonably good, the predicted results showing a somewhat greater effect 
of annulus convection than the experimental ones. The inconsistency may 
be due to: 

1. Incorrect assumptions about the dimensions of the test 
combustor. 
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2. Incorrect assumptions about the heat loss from the casing. 

The casing temperature does, in this case, exert an important 
influence on the flame-tube wall temperature, since radiation 
from the flame- tube wall to the casing becomes significant in 
the absence of film coolinq. 

The steep slope of the '’radiation-only" curve was not well 
predicted. It is interesting to note that Milford and Spiers' own 
calculations also failed to predict this trend. This is clearly due to 
two-dimensional effects which become increasingly important as the 
convecti ve-heat -transfer components are reduced. Two-dimensional 
effects were not included in this calculation. 

Figure 34 shows cases 5 to 7 - those with film cooling only. 

At the lowest-slot gas velocity, the correlation predicts poorer cooling 
near the slot than was obtained in the experiment. At higher velocities, 
the effect of film cooling was found to be greater than in the experiment. 
Ayain this is probably due to incorrect assumptions about the dimensions 
of the slot. 

Figure 35 shows the results for combined convective and film 
.tooling. As expected, the predicted effect of film cooling at low slot- 
gas velocities is again lower than the experimental one. 

Conclusion 

The most striking differences between the calculated and the 
experimental curves are ?n their slopes. There appears to be a dominant 
two-dimensional heat-transfer mode that is causing the experimentally 
measured wall temperatures to fall off rapidly with distance downstream. 

This effect is most apparent in the radiation-only results (experimental 
curves A of Figures 33 and 34), but it must also be present in the other 
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experimental results. Since the experimental conditions of Reference 75 
are poorly documented, it is not possible to simulate this two-dimensional 
effect in any meaningful way. 

From the results of the present computer calculations it is 
concluded that the film- and convective-cooling models used are producing 
results that agree with the main trends of this experiment, ignoring 
undocumented two-dimensional effects. 

Overal 1 Test Cases 

Operation of the complete computer program has been checked by 
running five overall test cases having identical geometry and inlet con- 
ditions (Table II). The program options specified by NASA for these test 
cases are given In Table III. The results of the test cases are discussed 
in the following paragraphs. 

The pressure-loss characteristics of the diffuser are summarized 
on Figure 36. Subsequent figures are grouped in pairs by individual test 
case. For Test Case 1 the cumulative fraction of inlet air present in 
the flame tube and the fractions of inlet air contained in the residual 
jets, as functions of axial position, are shown on Figure 37. the axial 
distributions of the flame temperature, the flame-tube wall temperatures, 
and the cumulative fraction of fuel burned on Figure 38. Figures 39 
through 46 are similarly ordered. On Figure 47 are shown the axial dis- 
tributions of total and static pressures in the flame tube and the annuli 
for Test Case 3, as a typical example of pressure distributions. 

The complete input and output data for Test Case 3 may be found 
in Appendices II and III of Volume II. 
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Diffuser Performance 

Figure 36 portrays the pressure-loss characteristics of the 
diffuser as obtained from the five diffuser calculation options employed 
in the test cases. The range of flow splits obtained for each option 
is that which occurred during the air-flow iterations in the process of 
obtaining a solution (the narrow ranges shown for Cases 1 and 2 are a 
result of accurate initial estimates of the flow splits). The large 
difference between Case 1 and Cases 2 and 3. each of which employs the 
streamtube method in the diffuser section between the compressor discharge 
and the snout, is due to the fact that when the streamtube method is used 
in the annuli, the calculated blockage decreases considerably between 
Stations 2 and 3. and hence small mixing losses result. In Case 2, 
where separation occurs between Stations 1 and 2, the blockage at Station 
3 is assumed to be equivalent to that at Station 2, and hence rather 
large mixing losses result. In Case 3» the mixing loss again includes 
the loss associated with the blockage at Station 2. The generally 
lower pressure losses in the inner annulus indicated by Cases 1, 2, 
and 3 as compared to Cases 4 and 5 are due to the distorted velocity 
profile at the compressor discharge which results in a higher local 
total pressure in the inner annulus; only when the streamtube method is 
used between Stations 1 and 2 is this distortion in total pressure re- 
flected in the annulus calculations. In the preset: t examples, which 
result in a badly mismatched flow at the snout, the use of the stream- 
tube method between Stations 1 and 2, and the mixing method between 
Stations 2 and 4 (as in Case 3). is recommended. 

The pressure losses in the diffusing passages influence the mass- 
flow split at the snout. Table IV summarizes the fractional flows through 
the snout and the inner and outer annuli for the two test cases. The 
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dependence of flow split on pressure drop may be seen by comparing entries 
in Table IV with corresponding points on Figure $ 6 . 

The flow split at the snout in turn influences the temperature 
level in the flame tube. As the amount of fuel available for burning in 
the primary zone does not vary from case to case, and as the total 
temperature of the entering air is substantially constant, the effect of 
a decrease in flow through the dome into the primary zone is to raise 
the temperature of the combustion products. This trend is evident in 
Table IV, where the primary-zone temperatures are listed. The primary- 
zone temperature influences the general level of temperature in the flame 
tube for some distance downstream. 

The temperature and flow distributions for the individual 
test cases are discussed in subsequent sections. Reference should be made 
to Table III for the program options used in running these cases. 

Test Case 1 

The choice of entrainment constants for this case has a 
pronounced effect on the results. ^ As shown on Figure 37, the penetration 
jets dissipate before the first calculation point downstream of their 
origin; this is to be expected for the entrainment constant (1) used here. 
The wall jets, however, with an entrainment constant (0) that implies no 
mixing, remain unchanged from their origin to the end of the flame tube. 
The virtually instantaneous mixing of the penetration jets results in a 
sharply falling flame temperature immediately downstream of each hole 
row (Fig 38 ) • As the wall jets do not mix at all, the amount of gas 
available in the flame tube is reduced by the contents of the wall jets 
and the resulting flame- temperature distribution is considerably higher 
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than it would have been, had the wall jets mixed with the main stream in 
the flame tube. 

In the program option used for this test case, film cooling is 
ignored, (it should be recalled that there is no interconnection in the 
program between wall jets and film cooling.) Oespite the lack of film 
cooling, the general level of the flame-tube wall temperature is relatively 
low. The explanation is, of course, the assumption of a nonluminous flame. 
The temperatures of the flame-tube walls follow the general variation of the 
flame temperature reasonably closely, as would be expected with a one-dimen- 
sional calculation of radiation from the flame to the walls, up to an axial 
position of about 24 inches. Downstream of this position the sharply 
dropping radiative beat transfer from the flame is roughly balanced by 
a drop in convective heat transfer in the annuli caused by sharply reduced 
ai r veloci t ies in the annuli downstream of the dilution holes. Near 
the end of the combustor the air velocities in the annuli fall nearly to 
zero, forcing the wall temperatures to rise sharply. 

Test Case 2 

Although the change in the wow split at the snout for this case 
causes a drop of about 140 deg F in the flame temperature at the end of 
the primary zone, the most pronounced changes in the temperature distri- 
butions for this test case, when compared with those of Test Case 1, are 
attributable to the change in the entrainmbnt constant for wall jets 
and the inclusion of film cooling. 

As in Test Case 1, the penetration jets effectively mix 
instantaneously (Fig 39). The wall jets, however, dissipate within 
about lj inches of their origins and their cool air is mixed into the 


149 


flame-tube gases. As a result, the flame temperature drops sharply down- 
stream of the cooling slots. Because of the lower primary-zone temperature 
and the more rapid jet mixing, the flame- temperature profile and the peak 
temperature are much lower than in Test Case 1. 

The inclusion of the film-cooling effect and the relatively low 
flame temperature combine to hold the flame-tube wall temperatures (Fig 40) 
at a general level about 400 deg F lower than in Test Case i. The wall- 
temperature profile is quite jagged in this case: essentially step changes 

in temperature at the cooling-slot posi tions . are followed by rapidly 
increasing temperature as the cooling film dissipates. 

Test Case 3 

The most significant change in program option in this case is 
from nonlumi nous to luminous flames. The radiation heat transfer for Case 3 
is greater by about a factor of three than that for Case 1. In addition, the 
behavior of the jets is substantially altered by the change in jet-mixing model 
(Fig 41). The penetration jets now remain distinct from two to four inches 
downstream of their origins; for the model used here (equivalent entrainment) 
the entrainment constant (1) implies complete jet dissipation when the mass flow 
in the "actual" jet has grown to twice the mass flow at the hole of origin. 

The wall jets mix very slowly and are extant at the end of the flame tube. 

Because of the relatively slow mixing of the jets, the flame 
temperature (Fig 42) is relatively high throughout the combustor, compared 
to that in Case 2, where mixing was rapid. The high flame temperature and 
luminosity combine to increase the flame- tube wall temperatures to a general 
level even higher than those for Case 1. The predominance of radiation 
as a heat- transfer mechanism is apparent from the manner in which the wall 
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temperatures decrease in step with the flame temperature in the downstream 
end of the combustor. 


Test Case 4 

The luminous- flame correlation used in this case gives luminosities 
considerably lower than were found in Case 3 (about 0.25 vs about 0.4). This 

accounts for the generally lower level of wall temperatures in this case 

(Fig 44), despite the exclusion of film cooling. The combination of no film 

• • • 

cooling and two-dimensional radiation causes the wall temperature profile 
to be relatively flat, compared to the profiles of previous cases. 

Cases 3 and 4 use the same mixing model, with changes in the 
entrainment constants for penetration and wall jets too small to affect jet 
behavior significantly (compare Figs 41 and 43) • For this reason the 
flame- temperature profiles of these two cases are quite similar. 

The calculated effect of radiation interchange between the flame- 
tube walls is negligibly small. 

Test Case 5 

With the entrainment constants and the mixing model used, the 
penetration jets mix in a distance too short, for a residual jet to 
be observed at the next calculation station downstream of the jet origin. 

The wall jets survive for roughly the distance to the next penetration-hole 
row (Fig 45). As a result, the flame- temperature profile is relatively low, 
with periodic bumps attributable to the influence of lingering wall jets. 

The luminosity correlation and relatively low flame temperatures 
result in low wall temperatures (Fig 46). The typically jagged profiles for 
wall cooling are tempered somewhat by two-dimensional radiation. 
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It should be noted that the effects of including radiation inter- 
change between the flame- tube walls and conduction along them are small and 
cannot be distinguished on these plots because of the stronger affects of 
changes in luminosity and jet-mixing correlations. 

Cone 1 us ions 

The results of the test cases indicate the sensitivity of 
calculated conditions in the combustor to the models used. The diffuser options, 
mixing and f lame- lumi nos i ty correlations, and the presence or absence of film 
cooling, can influence flow and temperature distributions dramatically. 

Choice of reasonable correlations is mostly a matter of judgment. A 
poor choice can result in ridiculous values. For example, a poor choice 
of entrainment constant in the prof I le- subs ti tut ion model can result in 
wall jets that do not disappear, or worse yet, grow with downstream distance. 

It has not been possible, within the scope of the present program, 
to investigate fully the influence of various arbitrary program Inputs, such 
as entrainment constants, on program results. Such an investigation should 
be carried out by the user, preferably in conjunction with appropriate 
experimental data, to select reasonable values of these constants. 



CONCLUSIONS 

The main conclusions of this work are: 

1, A computer program for analysing the aerodynamic and heat- 
transfer characteristics of annular and rectangular gas- 
turbine combustors has been completed, 

2, All computer- prog ram options have been checked and are found 
to operate satisfactorily, 

3, The diffuser, air-flow, and heat-transfer sections of the 
computer program have been tested individually by comparing 
their predictions of appropriate quantities with experimental 
data available in the open literature. The correlation between 
calculated and experimental values is judged to be reasonably 
good, 

4, The complete computer program has been tested for satisfactory 
operation by running a series of five hypothetical overall test 
cases, using geometry, inlet conditions and program options 
specified by NASA, The main trends of the results, as inferred 
by comparing changes in pressures, flows, and temperature distri- 
butions from case to case, are reasonable. It is therefore 
concluded that the program as a whole is operating properly, 

5, The wide range of program options available makes this program 

a useful tool for predicting the effect on combustor performance 
of changes in combustor geometry, inlet conditions, and fuel 
distribution. 
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6. The results of the stream tube diffuser calculation are 
surprisingly good in view of the gross assumptions made* 

If the equations are examined further and the potentialities 
and limitations of the method explored, this technique could 
provide a powerful new tool for more general application. 
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TABLE I - VALUES OF CONSTANTS IN GAS -PROPERTY EQUATION 

The gas-property equation (Eq 164) is: Z ■ a + bT ♦ cT^ + dT^ where Z 
is the gas property. Values of the four constants ere given in the table 
for various gases for the gas properties indicated. 


Specific Heat 


Gas 

a 

b x 10 4 

QO 

O 

X 

o 

d x 10 n 

Range 
Deg R 

Per Cent 
Accuracy 

Air 

0.2419 

-G r08l8l 

1791 

-0.2743 

180-3600 

1.9 

Carbon Dioxide 

0.1397 

1 .349 

-2.646 


400-2500 

1 .1 

N i trogen 

0.2540 

-0.2067 

2.658 

-0.4718 

200-3600 

1.9 

Oxygen 

0.2126 

0.1564 

0.9044 

-0.2294 

700-2700 

1.1 

Water Vapor 

0.5632 

-2.346 

18.51 

-3.351 

700-2700 

1.1 


V i scos i ty 


Air 

-0.2853 

3.268 

-0.8253 

1.239 

200-1800 


Carbon Dioxide 

-1.180 

1.582 

0.7728 


400-1100 


N i trogen 

-0.01937 

3.194 

-0.7820 

0.8990 

200-2200 

0.5 

Oxygen 

-0.4704 

3.270 

-0.5915 


200-1100 


Water Vapor 

-2.316 

2.280 

0.1756 


700-1500 



Thermal Conductivity 


Air 

3.057 

8.607 

-2.279 

2.908 

700-3200 

1.8 

Carbon Dioxide 

-1.203 

7.786 

-1.796 

2.320 

400-3000 

0.5 

Nitrogen 

2.142 

8.838 

-2.766 

4.295 

400-3600 

2.5 

Oxygen 

4.770 

9.024 

-1.994 

2.368 

400-3600 

3.0 

Water Vapor 

6.200 

2.442 

1 .926 

-4.428 

700-2700 

1.6 
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TABLE II - COMBUSTOR GEOMETRY AND INLET FLOW 
CONDITIONS USED FOR OVERALL TEST CASES 


Configuration of Combustor Walls 


Axial 
Pos i t i on 
From 

Compressor 
Discharge 
1 nches 

1 nner 
Casing 

Inner 

Snout 

Inner Dome 
or Flame 
Tube Wall 

Outer Dome 
or Flame 
Tube Wail 

Outer 

Snout 

B 

Geometric 

Input 

Point 

Number 

0 . 

32.500 

-0. 

-0. 

-0. 

-0. 

WfgM 

1 

0.500 

32.300 

-0. 

•0. 

-0. 

-0. 


2 


32.000 

-0. 

-0. 

- 0 . 

-0. 

36.200 

3 

I 1.500 

31.800 

-0. 

- 0 . 

-0. 

"0. 

36.400 

4 


31.500 

-0. 

-0. 

-0. 

-0. 

36.700 

5 


31.300 

-0. 

-0. 

-0; 

-0. 

36.900 

6 

2.800 

31.100 

33.100 

- 0 . 

-0. 

35.100 

37.100 

7 

3.000 

31.000 

33.000 

- 0 . 

-0. 

35.200 

37.200 

8 

4.000 

30.500 

32.500 

-0. 

-0. 

35.700 

37.700 

9 


30.100 

32.100 

- 0 . 

-0. 

36.100 

38.100 

10 

6.000 

29.600 

31.600 

- 0 . 

-0. 

36.600 

38. 600 

11 


29.100 

31.100 

34.100 

34.100 

37.100 

39.100 

12 


28.600 

30.600 

31.100 

37.100 

37.600 

39.600 

13 

8.800 

28.200 

30.200 


38.000 

38.000 

40.000 

14 


28.200 

30.2C0 

30.200 

38.000 

38.000 

40.000 

15 

9.100 

28.200 

29.900 


38.300 

38.300 

40.000 

16 

14.000 

28.200 

29.900 

29.900 

38.300 

38.300 

40.000 

17 

15.000 

28*200 

29.900 

29.900 

38.300 

38.300 

40.000 

18 

15.100 

28.200 

29.600 

29.600 ' 

38.600 

38.600 

40.000 

19 

20.000 

28.200 

29.600 

29.600 

38.600 

38.600 

40.000 

20 

21.000 

28.200 

29.6Q0 

29.600 

38.600 

38.600 

40.000 

21 

21.100 

28.200 

29.300 

29.300 

38.900 

38.900 

40.000 

22 

26.000 

28.200 

29.300 

29.300 

38.900 

38.900 

40.000 

23 

27.000 

28.200 

29.300 

29.300 

38.900 

38.900 

40.000 

24 

27.100 

28.200 

29.000 

29.000 

39.200 

39.200 

40.000 

25 

29.000 

28.200 

29.000 

29.000 

39.200 

39.200 

40.000 

26 

30.000 

28.200 

29.000 

29.000 

38.200 

38.200 

40.000 

27 

31.000 

28.200 

29.000 

29.000 

37.200 

37.200 

40.000 

28 

32.800 

28.200 

29.000 

29.000 

37.200 

37.200 

40.000 

29 


Swlrler Design (Specif lad as Input) 

Number of Swirlers • 24 

Number of Blades • 8 

Blade Stagger Angle * 45.00 Degrees 

Inner Diameter *1.00 Inches 

Outer Diameter • 1.50 Inches 

Area Per Swirler • 0.98 Square Inches (ignoring Blockage 

Due to Vehes) 
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TABLE II. CONTINUED - COMBUSTOR GEOMETRY AND INLET FLOW 
CONDITIONS USED FOR OVERALL TEST CASES 


Details of Holes and Cooling Slots 


Hole 

Row 

Number 

Axi Jl 
Position 
of Hole 
Center- 
Line 
1 nches 

Inner 

or 

Outer 

Wall 

Number 
of Holes 
i n 
Row 

Description of Hole 
or Cool ing Slot 

1 

9.000 

■a 

1 

\ Cooling slot: stepped louvre 0.095 in 

2 

9.000 


1 

i high, with wiggle strip 

3 


wtSM 

50 

\ Round holes 0.75 in diameter with scoop 

4 

12.000 

Outer 


/ 0.75 in wide and 0.636 in high 

5 

13.500 


o 

Dummy hole row (no holes) 

6 

^MlWi yiyn 

Inner 


1 vCooling slot: stepped louvre 0.095 in 

7 

15.000 

Outer 


I high, with wiggle strip 

8 

16.500 



Dummy hole row (no holes) 

9 


I nner 

■ ' l ■ 

1 Flush holes 2.0 in long and 1.0 in wide 

10 


Outer 

jjTv ' iy 

! 

11 



0 

Dummy hole row (no holes) 

12 

21.000 

Inner 

1 

) Cooling slot: stepped louvre 0.095 in 

13 


Outer 

1 

1 high, with wiggle strip 

14 

22.500 


0 

Dummy hole row (no holes) 

15 

24.000 

Inner 

75 

1 Flush holes 2.0 in long and 1.0 in wide 

16 

24.000 

Outer 

75 

) 

17 

27.000 

Inner 

1 

1 Cooling slot: stepped louvre 0.095 In 

18 

27.000 

Outer 

1 

1 high, with wiggle strip 

19 

32.700 


0 

Dummy hole row (no holes) 


Inlet Flow Conditions 


Total Temperature at Compressor Discharge ■ 
Total Pressure at Compressor Discharge * 
Air Flow Rate at Compressor Discharge ■ 
Overall Fuel-Air Ratio • 


1150.000 Deg F 

90.000 psia 

96.000 Ibm per sec 
0.018 
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Program Option 


Test Case Number 


a 

a 

a 

b 



Calculation method in first part of 
di f fuser: 

a. Streamtube 

b. empirical Data 

Calculation method in second part 
of diffuser: 

a. Streamtube and Mixing 

b. Empirical Data and Mixing 

c. Sudden Expansion and Mixing 

The following combinations of dif- 
fuser calculations are permitted: 

First Part a 

Second Part a 

Entrainment function for jet mixing: 

1. Direct Mass Loss 

2. Equivalent Entrainment 

3. Profile Substitution 

Entrainment constant for penetration jets 

Entrainment constant for wall jets 

Luminosity correlations for flame 
emissivity: 

1. Reeves Distillate Fuel 
(nonl umi nous) 

2. Reeves Residual Fuel 
(nonl umi nous) 

3. Lefebvre (luminous) 

4. NREC 1964 (luminous) 

5. NREC 1966 (luminous) 

One- or two-dimensional radiation? 
Radiation interchange between walls? 
Cooled walls? 

Longitudinal wall conduction? 

Heat transfer to annulus air? 
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TABLE IV - FLOW SPLITS AND PRIMARY-ZONE TEMPERATURES FOR OVERALL TEST CASES 


Test Case No. 


Fraction of Inlet Air 
Entering Snout or Dome 

.04847 

.06174 

.06230 

.05759 

.06313 

Fraction of Inlet Air 
Entering Inner Annulus 

.43117 

.42020 

.42994 

! 

.40121 

.41369 

Fraction of Inlet Ai r 
Entering Outer Annulus 

.52036 

.51806 

" 1 
.50775 

.56120 

.52317 

Temperature of Gases in 
Primary Zone, deg F 

3890. 1 

3730. 1 

3758.8 

3800.0 

3725.7 
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FIGURE 1 - OVERALL FLOW CHART FOR ENTIRE PROGRAM 
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FIGURE 2 


OVERALL FLOW CHART FOR DIFFUSER SUBPROGRAM 
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DISJET 
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After 
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Conditions at Secondary 
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Determine conditions in 
the flame tube at the 
secondary holes 


Solve annulus Eqs 
for m . u .p p , T 


PRTEMP 

Calculate the tem- 
perature in the 
primary zone 


— L 

Flame Tube 


EQUFT 


JETMIX 

Move from calcula- 
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Solve flame tube 
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— < 
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jets 
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HEATAD 
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burning fuel 


NOTE: Fortran titles of subroutines 

are defined in Volume II. 
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FIGURE 4 - OVERALL FLOW CHART FOR HEAT-TRANSFER 
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Nond i mens i ona 1 Length, L/Wj 

— Experimental results of Reneau, Johnston, and Kline 
(Ref 7 ) for 0.015 inlet blockage 

— Streamtube method with 

Initial shape factor = 1.4 

Shape factor at separation ** 1.9 

Inlet blockage = 0.015 

FIGURE 9 - COMPARISON OF THE EXPERIMENTAL RESULTS OF 
RENEAU, JOHNSTON. AND KLINE FOR TWO-DIMENSIONAL DIFFUSERS 

WITH THE STREAMTUBE METHOD 
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FIGURE 10 ■ SWIRLER DATA FOR 16 COMBUSTOR DESIGNS 




















Pressure atm 


Temperature T 9 deg F 


From Hodge (Ref 37) 


FIGURE 12 - DISSOCIATION CORRECTION TO EFFECTIVE FUEL CALORIFIC VALUE 
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From Reference 38 


FIGURE 1 3 - PRESSURE LOSS IN A SUDDEN EXPANSION IN COMPRESSIBLE FLOW 
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Distance along Jet Center-Line s'/d. 


From Keffer and Baines (Ref 56 ) 


FIGURE 15 - DECAY OF CENTER-LIN E VELOCITY IN 

PENETRATION JET 









Dimensionless Di stance Across Prof i le ^ /K)x 

2 

From Keffer and Baines (Ref 56) 


FIGURE 16 - L ATERAL DISTRIBUTION OF VE LOCITY 
IN PENETRATION JET 
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Numbers on curves indicate 
Figure ir. Reference 57 from 
which data were taken. 
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Jet Wi dth T|, ft 
Data from Jordinson (Ref 57) 


FIGURE 17 - CORRELATION FOR CROSS-SECTIONAL 


AREA OF PENETRATION JETS 








Local Jet Velocity, ft per sec 
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FIGURE 18 - TYPICAL VELOCITY PROFILES IN A WALL JET 
























Pressure p^ t , atmospheres 

Data from Schirmer arid Quigg (Ref 11) assuming 
T ft « 3500 deg R 

Values predicted from NREC 1964 correlation 
(Equations 134 and 136) 


FIGURE 19 - EFFECT OF PRESSURE AND FUEL 
COMPOSITION ON EM1SSIVITY OF LUMINOUS FLAMES J 











Pressure p^ t , atmospheres 

Data from Schirmer and Quigg (Ref 11) assuming * 3500 deg R 

Values predicted from NREC 1966 correlation (Equations $37 
and 138) 

FIGURE 20 ~ A NEW CORRELATION FOR THE 
EMISSiVlTY OF LUMINOUS FLAMES 
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FIGURE 21 - EFFECT OF PRESSURE AND 
GAS THICKNESS ON RADIATION TRANSM l TTANCE 
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FIGURE 22 - CORRELATIONS OF FILM-COOLING DATA 
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Diffuser Effectiveness £ 


Nondimensionat Length, L^R 


Experimental resuits of Sovran and Klomp (Ref 6) 


Sfreamtube method with: 
Hub/Shrouc ratio at inlet 
Initial blockage 
Initial shape factor 
Shape factor at separation 


FIGURE 23- COMPARISON OF THE EXPERIMENTAL RESULTS OF SOVRAN 
AND KLOMP FOR ANNULAR DIFFUSERS WITH THE STREAMTUBE METHOD 
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Total Included Angle, 20 

A— -A Experimental results of Cochran and Kline (Ref 70) 

Streamtube method with: 

Inlet blockage » 0.0075 

Inlet shape factor ■ 1.47 


FIGURE. 24 - COMPARISON OF EXPERIMENTAL AND THEORETICAL 
EFFECTIVENESS FOR TWO-DIMENSIONAL DIFFUSERS 
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Nondimensional Length, L/W. 

* ■ 

Experimental results of Moore and Kline (Ref 72 ) 

> 

Streamtube method 

FIGURE 25 - POSITION OF FIRST STALL 
FOR TWO-DIMENSIONAL DIFFUSERS 
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Total Included Angle, 0, deg 

Experimental results of Abramovich and Vasil'yev 
(Ref 70 for area ratio of 2.0 

— Streamtube method with: 

Inlet shape factor * 1*^ 

Shape factor at separation * 1.9 


FIGURE 26 - COMPARISON OF EXPERIMENTAL AND 


| CAL EFFECTIVENESS FOR ANNULAR DIFFUSERS 
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Total Included Angle, 9, deg 


Experimental results of Abramovich and Vasil 'yev 
(Ref 71) for area ratio of 3.8 


Streamtube method with: 

Initial shape factor * 1.4 

Shape factor at separation ■ 1.9 


4 


FIGURE 26 CONTINUED 
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FIGURE 27 - COMPARISON OF EXPERIMENTAL AND THEORETICAL 
DEVELOPMENT OF A SHEAR FLOW IN A TWO-DIMENSIONAL DIFFUSER 
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FIGURE 29 - VARIATION OF PRESSURE-LOSS FACTOR 
WITH REFERENCE MACH NUMBER 
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FIGURE 30 - VARIATION OF PRESSURE-LOSS 
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FIGURE 31 - EFFECT OF HEAT RELEASE ON 
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FIGURE 33 - WALL TEMPERATURES. CONVECTIVE COOLING ONLY 
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FIGURE 37 - NASA TEST CASE NO. 1 — ACCUMULATED FRACTION OF INLET AIR IN FLAME TUBE 
; AND RESIDUAL AIR IN JETS AS FUNCTIONS OF AXIAL POSITION IN COMBUSTOR 
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FIGURE 39 ~ NASA TEST CASE NO. 2 -- ACCUMULATED FRACTION OF INLET AIR IN FLAME TUBE 


AND RESIDUAL AIR IN JETS AS FUNCTIONS OF AXIAL POSITION IN COMBUSTOR 
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FIGURE 41 - NASA TEST CASE NO. 3 — ACCUMULATED FRACTION OF INLET AIR IN FLAME TUBE 


AND RESIDUAL AIR IN JETS AS FUNCTIONS OF AXIAL POSITION IN COMBUSTOR 
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SURE 42 - MSA TEST CASE NO. 3 — FLAKE TEMPERATURE. FLAME-TUBE WALL 
TEMPERATURES. AND FRACTION OF FUEL BURNED AS FUNCTIONS 
OF AXIAL POSITION IN COMBUSTOR 



Cumulative Fraction of Fuel Burned 






Fraction of Inlet Air In Flame Tube and In Residual Jets 



X Inner Wall Holes 

O Outer Wall Holes 

Air In Flame Tube as 

Fraction of Inlet Air 
(Includes Residual Jets) 

Air in Residual Jets 


8 12 16 20 2k 

Distance Downstream of Compressor Delivery, in 


FIGURE 43 ~ NASA TEST CASE NO. 4 — ACCUMULATED FRACTION OF INLET AIR IN FLAME TUBE 
AND RESIDUAL AIR IN JETS AS FUNCTIONS OF AXIAL POSITION IN COMBUSTOR 
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FIGURE 44- NASA TEST CASE NO. 4 -- FLAME TEMPERATURE. FLAME-TUBE WALL . 
TEMPERATURES. AND FRACTION OF FUEL BURNED AS FUNCTIONS 
OF AXIAL POSITION IN COMBUSTOR 
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FIGURE 45 - NASA TEST CASE NO. 5 — ACCUMULATED FRACTION OF INLET AIR IN FLAME TUBE 
AND RESIDUAL AIR IN JETS AS FUNCTIONS OF AXIAL POSITION IN COMBUSTOR 
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FIGURE 47 - NASA TEST CASE NO. 3 -- TOTAL AND STATIC 
PRESSURES AS FUNCTIONS OF AXIAL POSITION 
IN THE COMBUSTOR 







NOMENCLATURE 


Description Units 

Area - sq ft 

Area occupied by boundary- layer sq ft 

displacement thickness 

Hole area sq ft 

Cross-sectional area occupied by a sq ft 

specified mass flow at diffuser 
station 2 

Wetted wall area per unit length ft 

Critical area, for which Mach number sq ft 


equals 1*0 

Area Ratio A2/AJ 

Entrainment constant 

Rate of heat transfer by convection 
from the hot gases (or the cooling 
film) to the flame-tube wall 

Rate of heat transfer by convection 
from the wall to the annulus air 

Inner wall area per unit length 

Outer wall area per unit length 

Discharge coefficient 

"Corrected" discharge coefficient 

Specific heat 

Pressure-recovery coefficient 

Ideal pressure- recovery coefficient 

Ideal pressure- recovery coefficient 
in the presence of mixing 

Hydraulic diameter of annulus ft 


Btu per Sq ft sec 

Btu per sq ft sec 

ft 

ft 

Btu per Ibm deg F 


216 


Symbol 

Description 

Units 

D ft 

Wall-to-wall width of flame tube 

ft 

°i 

ft 

Hydraulic diameter of flame tube 

(** 2Df t ) 

ft 

. « 

D ref 

Combustor reference diameter 

ft 

d h 

Effective axial length of hole 

ft 

d j 

Effective initial jet diameter 

ft 

l-E 

Diffuser blockage 

- 

*“ E bi 

Diffuser boundary- layer blockage 

m 

f 12 

View factor of receiving surface 
2 from radiating object 1 

- 

F 

Fanning friction factor 

. m 

f 

Fuel -air ratio 

m 

G 

Mass flow rate per unit area 

ibm per sq ft sec 

9o 

Constant in Newton 1 s law 

ft Ibm per lbf sec 2 

h 

Enthalpy 

Btu per Ibm 

H 

Boundary- layer shape factor 

m 

H, 

Shape factor at diffuser inlet 

- 

d sep 

Shape factor at which separation occurs 

• 

H/C 

Fuel hydrogen-carbon ratio 

m 

^eff^p 

Effective fuel lower calorific value 

Btu per Ibm 

H p 

Fuel lower calorific Value 

Btu per Ibm 

J 

Mechanical equivalent of heat 

ft lbf per Btu 

K d 

Number of dynamic heads lost in flow 
from diffuser statiop 2 or 2 1 to dome 

** 

AK 

Net rate of heat transfer by conduction 
i nto unit: area of the f 1 ame*tube wail 

Btu per sq ft sec 


from adjassrtt wall elements 



2*7 


Prandtl number 
Static pressure 

Partial pressure of radiating gas 

Heat transferred .to transpiration 
coolant as it passes through the v 

-2 

Dynamic head* .puj^EgQ 

Rate of heat add i t ion due to fuel 
burning V 


Heat transfer .rate from the flame* 
wait to the annulus air 

Gas. constant 


Units 


Gas thermal conductivity Btu per ft sec deg F 

Thermal conductivity of wall material Btu per ft sec deg F 

Length ft 

Mean beam path length ft 

Mach number 

Molecular weight of coolant 

Axial component of jet momentum flux ft Ibm per sec 

Mass fraction of component In mixture 

Mass flow rate Ibm per sec 

Mass flow rate through hole Ibm per sec 

Number of streamtubes in diffuser 
analysis 

Mussel t number « 

Total pressure Ibf per sq ft 

Loss of total pressure due to combustion Ibf per sq ft 


ft Ibm per sec 2 


Ibm per sec 
Ibm per sec 


Ibf per sq ft 
Ibf per sq ft 
Btu per sq ft sec 


Ibf per sq ft 
Btu per sec 


Btu per sq ft sec 


ft Ibf per 
deg R 
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Symbol 

R| 

r 2 

r 3 


liS 


r 

Re 

^ e cf 




s 

s' 


T 

T o 

T b 

Tc 

Tcf 

Tft 

Tw 

TA 

■V 

U 

U ft 


Description Units 

Rate of heat transfer by radiation from Btu per sq ft sec 
the hot gases to the flame-tube wall 


Rate of heat transfer by radiation from Btu per sq ft .sec 
the wall to the casing 


Net rate of heat transfer by radiation Btu per sq ft sec 

from the flame-tube wall to all parts 
of the opposite wall 

Diffuser inlet width, for annular diffusers ft 
Radius ft 


Reynolds number 

Reynolds number of cooling film, based 
on slot height 

Reynolds number of gas in flame tube, 
based on flame-tube hydraulic diameter 

Distance along jet center- line 

Distance along jet center-line from 
virtual origin 

Static temperature 

Stagnation temperature 

Adiabatic-wall temperature 

Base temperature for enthalpy 

Casing temperature 

Temperature of cooling air 

Static temperature of hot gases 

Wall temperature 

Total cross-sectional area of streamtubes 

Wall thickness 

Main-stream velocity 

Bulk velocity of gas in flame tube 


ft 

ft 

deg R 
deg R 
deg R 
deg R 
deg R 
deg R 
deg R 
deg R 
sq ft 
ft 

ft per sec 
ft per sec 



Description 


Units 


Symbol 


u 


W 

X 


*o 

x 

X 

X 

Ax ft 


4*w 

Y 

y 

Yfif 

y 

ttft 

rt w 

K 

/S 

r 

5 

* 


Mean initial velocity of cooling film 
Velocity 

Area-average velocity 

Mass-averaged velocity 

Diffuser width, for rectangular diffusers 

Nondimensional distance downstream of 
cooling slot 

Parameter in film-cooling correlation 
Mole fraction 
Downstream distance 
Distance downstream of cooling slot 

Thickness of slab of hot gas 

Axial length of strip of wall 

Nondimensional distance across diffuser 

Distance across diffuser or flame tube, 
measured from the wall 

Height of film-cooling slot 

Streamline slope 

Flame absorptivity 

Wall absorptivity 

Permeability coefficient 

Profile parameter 

Ratio of specific heats 

Boundary- layer displacement thickness 

Flame emissivity 

Transverse jet -coordinate 


ft per sec 
ft per sec 
ft per sec 
ft per sec 
ft 


ft 

ft 

ft 

ft 

Ml 

ft 

ft 

radians 

a* 

sq ft 

m 

ft 

m -, ... . 

ft 
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Symbol 


Y|* 

6 

A 

^diff 


(A 

V 

<P 




Description 

Value of l)at which nondimens ional 
velocity profile parameter equals 0.5 

Value of *Jat the assumed jet boundary 

Boundary- layer momentum thickness 

Luminosity factor 

Diffuser total -pressure- loss coefficient 
Dynamic viscosity 
Kinematic viscosity 
Initial jet angle 

Diffuser effectiveness in the presence 
of mixing 

Diffuser effectiveness between Stations 
1 and 2 

Density 

Stef an^Bol tzmann constant 
Fi Im-cool ing effectiveness 
Transmittance 


Units 

ft 

ft 

ft 

m 

Ibm per ft sec 
$q ft per sec 
radians 

m 

m 

Ibm per cu ft 

Btu per sq ft-(deg R) 
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APPENDIX I 

IDEAL PRESSURE- RECOVERY COEFFICIENT IN A DIFFUSER 
WITH NONUN I FORM INLET VELOCITY PROFILE 


Governing Equations 

Consider the incompressible flow in a diffusing passage; it is 
assumed that the velocity profile is nonuniform at the inlet, and mixes 
to a uniform profile at the outlet as indicated in the sketch* 



u. 


The equations of continuity and momentum for an elemental control 
volume of length dx can be written as 


d' 0® “/) » o 


0 - 1 ) 


dp 

dx 


-44- (/>jj u A) 
g A dx am 
o 


0 - 2 ) 


Equation 1-1 can be integrated to yield 


u, 


m 


jou £ m h 0 - 3 ) 

where m is the mass flow In the passage* Combining Equations 1-3 and 1-2 
yields . ‘ 


dx />.S 0 ‘ 


mL d_ (A \ 

A dx V A r 



This is the equation which expresses the static pressure as a function 
of the passage area and the uniformity of the velocity profile* it is 
to be noted that for a uniform profile j3 m 1, while for a nonuniform profile 
JB>1. 


Solution fo 


mum P 




To obtain a solution for the maximum possible pressure recovery, 
it is first necessary to determine the form of f5 (x) which yields the 
maximum value of the function 


"f Z d * 

lx. ***’■■ 


(1-5) 


This expression can be integrated by parts to obtain 

X. m « /•■S-l * \ <rv mm dA (1-8) 

• ^ \ a) r ■ 

The maximum value of X will accordingly be obtained when the last term in 
Equation i-6 is a minimum; this Is obviously achieved by having reach its 
minimum value instantaneously at Station 1* Hence, the maximum possible 
pressure recovery is achieved by instantaneous mixing at Station 1, i*e«, 


then 


fi m at x * Xj 


The solution of Equation 1*4 for 


*2 *2 

m m 

P , - “ i - 


pressure recovery 


1 1 1 


‘■2 2A, 2A 2 
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or* from Equation 1-3 * 



+ 1 - jjp (1-5) 

m 



APPENDIX 11 


BOUNDARY- LAYER MOMENTUM THICKNESS 

Introduction 

In the present streamtube method, a relationship is required 
between the boundary- layer momentum thickness along the walls of the 
diffuser and the local flow variables* The situation considered is 
indicated in the following sketch: 


Boundary Layer 


The method used herein is a modification Of standard integral 
treatments of the turbulent boundary layer (see Ref Chapter XXII) to 
account for compressibility effects. The integral momentum equation for 
compressible flow is (see Equation 1-23 of Kutateladze and Leont’ev, Ref 

« 

where /> _ » density in free stream 

8 “ f^u 0 - u 3 

v Q 

If it is assumed that H Is constant and that the shear stress at 


dd 

dx 


£ (H 

u vn 


* 2 ) 


dU 

dx 


A 


dx 
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where n and (X are constants and V is assumed to be constant, then 


Equation 2-1 may be integrated to gives 


h 

u. e 'n 


1 


(-^) * 


■ (l+*) -b 


U 


x 


+ a 


b n +l 
U b fi n dx 


(2-3) 


mJ, 1 

where . a ■ (k 
n 

b ■ (h+2) 


i 

n 


C| 88 constant of integration 


Taking values of n « 6 and c{ * 0.0065 (Ref 76), and assuming a mean value 
of H « 11/7, then 1-3 becomes: 


U 0 '4> 

0 


x ' V 


•) 


1 


'x u x L 


C| + 0.0076 


dX 


V»x. 


(2-4) 


0 . 


6 


^ ,V 25/6 . 0. 


6/7 


(~ZZ) 

'U ' 
x 


(f) 


..0076 V ,/6 1 U 4 D 7/6 , 

+ 7175T 25 7 r\ u / dx 

X .X 


(2-5) 


Equation 2-5 is used in the computer program to calculate 0; a constant 

1/6 1/6 
ue of V * *=0.23 (sq ft per sec) f is assumed. 
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APPENDIX Ml 

NEW ESTIMATE OP BOUNDARY- LAYER 
DISPLACEMENT THICKNESS 

* 

This appendix describes the technique used to obtain convergence 
in the iteration on bounda ry- 1 ayer thickness* 





Consider any section where the present value of displacement 

thickness is <5 and the predicted value is £ . The correct value lies 

o 

between 6 and 6 since, as 6 is increased, the flow velocities increase 
o o 

and 6 is reduced. The best estimate, 6 » to use as a starting value for 

the next iteration is such that, as £ is increased to § , S is reduced 

o 

to <5 \ Let 

A<S = <5' -6„ 

O O 


At any section, from Equation 2-4: 

® a “24/7 

u 


AB 24 Au 
6 " ” 7 u 

From the continuity equation for two-dimensional flow: 

Au A<£q| * A(f 0 2 
u ~ EA 


(3*1) 


(3-2) 
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where EA is the flow area and subscripts 1 and 2 refer to inner and 
outer wall s, respectively. 


By definition: 


S = He 


6 1 = (H +6H) (0 +A0) 

S' - 6 [' * f * Y 


(3-3) 


From Equation 36; 


AH 

AH 

H 


- 70. (H o - 1 .05) A (~) 

= 70. (H ° - ’-°5) _ (A 0. -a 6,.|) 


H 


"(x.-x. .) 

v i i-r 


(3-4) 


Combining Equations 3-l» 3-2, 3-3* and 3-4: 


<T - 6 


1 - 240. 


OL - 1.05) 


H 


+ 24/7 


x.-x. . / \ 

« i-l / v 


> 

| 4ioj *A£o2 


70(H q - 1.05) AS,., 


H 


X.-X. . 
i « = ! 


(3-5) 


i.e., for walls 1 and 2 


6, [i -f,l + ?, <* 0 , +<y o2 ) = (f, ♦ Dal +t P> 4 (3-6) 


1 "^2 +l fz (<5 o? +tf o2 } 


f 2 .<Sj + (f 2+ 1)4 (3-7) 


where Ac 5 q1 = cJ 1 - <S ( 

Cl) ^ JL 
T EA 


(H - 1 
o 


a 


H 


Wl 
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¥ - 70 


(H q - 1.05) 
H 


A0 


1-1 


x r x i-i 


Solution of Equations and 3*7 provides a new estimate of the boundary- 
layer displacement thickness; experience has indicated, however, that this 
estimate will not in general yield convergence and hence the actual estimate 
used Is a weighted average of <5 Q and as indicated in the text (Equation 


37 ). 
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APPENDIX IV 

ESTIMATES OF BOUNDARY* LAYER SHAPE FACTORS 


As inputs to the streamtube method of diffuser analysis, boundary 
layer displacement thicknesses and shape factors on the two walls are 
normally provided. If, however, momentum thicknesses or shape factors 
are not available, they may be easily calculated, assuming the boundary 
layer has built up as on a flat plate. 

The equation for the shape factor at inlet, Hj, is derived as 

f 

follows. Hama (Ref 78) » from experiments with boundary layers on various 
types of flat surfaces, produced the following correlation: 

H 1 (4-1) 

1 - 0.78 R ' * 

where R * Mjp 
x V 

x * length of flat plate 
U ■ free stream velocity 
V ■ kinematic viscosity 
For flat plates Equations 2-1 and 2-3 give: 


ft 

“5T “ 7ig ,/n 


(4-2) 


where R * M . u • £ 

wnere - e V V H 
I ntegrat i ng Equat i on 4- 2 : 

n*l 

“iff Rg n » 0t . R x (4-3) 

Using values of ft * 0.0065 and n « 6 (Ref 76 ), Equations 4-1 and 4-3 give 
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1-0.51 


1/12 



This equation may be solved to obtain » if the- inlet displacement 
thickness (or boundary* layer blockage) is known* 


(4-4) 






APPENDIX V 


DEVELOPMENT OF MIXING EQUATION FOR 
j Nfcffi'PRfe'SsYBLE FLOW * 

Governing Eguptlon*. 

The equations governing the Incompressible flow In a diffusing 
passage have been derived In Appendix I; they are: 


J*v 

i S 


m 


(5-0 


£1° da . . I 4- (A\ 

.2 dx A dx 'A' 
m 


(5-2) 


Mixing-Process Assumption 

To obtain a solution to Equations £-1 and £-2, it is necessary 
to postulate the form of the mixing process, which is characterized by y@ (x). 

a 

it is assumed here that mixing will occur in such a way that 



constant, K>yS 


(5-3) 


which is merely based on the plausible hypothesis that the mixing rate will 
be delayed by diffusion* The value of K may be obtained by applying Equation 
£-3 at the inlet and exit: 


K » 


5 


AR ^2 
1 

' AR : 


( 5 - 4 ) 


Solution for Pressure-Recovery Coefficient 

Substitution of Equation £-3 iflto Equation £-2 yields: 
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K 



which can be integrated between inTet and exit: 







Combining Equations 5"i* and 5-5 provides the desired result: 



\ 



J 2 + AR* ^1 " AR ^2* 


If the flow mixes to a uniform profile ( p ^ » I), then: 



(5-5) 


(5-6) 


(5-7) 
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APPENDIX VI 

PRESSURE DROP ACROSS THE DONE 

The total -to-statlc pressure drop across the dome results from 
the combined mass-flow rate through the swirler and through holes in the dome: 


& p d - f fa d ) 

where m^ = total flow through the dome (or snout) 

= m .. + m 
dh sw 

The mass-flow rate through the dome holes is given by the discharge equation; 

™dh = C dh A dh/ 2 ®o P ref 

The mass-flow rate through the swirler may be found from the swirler pressure- 
drop equation: 


A p 


^ ref 


sw _ 


A 


- K (- 


ref 


sw 


\ Z 

sw 


sec 2 j3 


A 


ref. 


A 


The reference dynamic pressure is given by: 


Vef 


_ /Vef U ref _ 


ft 


.2 

m ref 


.2 

m 

sw 

.2 

m ref 


(6-2) 


2 g. 


2 9 o J* ref A ref 


This expression is substituted for q in Equation 6-2 
solved for the swirler mass-flow rate: 


and the result is 



(6-3) 
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The total flow through' the dome may now be found by adding Equations 6-1 
and 6-3: 

m d . “ C dh A dh / 2a o ref AP d 

Since the pressure drop across the dome holes and the swirler must be equal, 
this equation may be solved for this pressure drop: 
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APPENDIX VII 

INTEGRATION OF INTEGRAL PRESSURE TERM 


The Integral pressure term in the annulus momentum equation 


(Equation 90) may be integrated by parts to yield; 


n dA . 

p c£ ox 


(pfl) 


A dp 


(7-1) 


An appropriate relation between A and p for mixing has been derived in 
Appendix V; 


(7-2) 


When Equation 7*2 is substituted into Equation 7*1 » the latter may be 
integrated; 


f dx.(pA) + ^4 2 -*J 7 > 


(7*3) 


Equation 7-2 may be integrated directly: 


■^ 9 o . / \ 

.2 ^ p 2 p r 

m 


K / l 


2 a. 2 
A 2 


(7*4) 


The desired result is obtained by eliminating K between Equations 7*3 and 

rZ * • ' . 


dA j 

p -r- dx « 


|* ;+ 2{p|‘* Pg)/^ * X J 

• 12 * ; -1 


(7-5) 


